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1. INTRODUCTION 


Although nonparametric procedures are powerful 
Leeks to the analyst, they are currently underused and 
Gfven avoided by potential users. Perhaps one reason 
Moe woes is the ditficulty in generating the exact 
G@werePOubLLONS of the test statistics, even for moderate 
sample Sizes. Consequently, tables eu these 
G@eertbutions are only available for very small sample 
“ges ana normal theory based approximations must then 
pe usc . 

ine pmiroose Of this thesis is to make a variety of 
nonparametric procedures quick, easy and accurate to 
apply using menu driven computer programs in APL.! 
These programs use enumeration, recursion, or 
combinatorial formulas to generate the exact null 
Gia rTroution of the various nonparametric test 
Stem lStics. fulgial ss allows hypotmesis testing and 


confidence iT mcer va estimation to be based on exact 


Gisetcoibutions without the use of tables. For larger 
sample sizes. the normal , a MG ee aeSeripubtions are 
TAPL was chosen because it is an interactive 


language that is especially powerful at performing 
Sruecibations dealing with rank order statistics and 
Veeuors arithmetic. Menus are not included in the 
workspace designed for the mainframe. 


used to approximate the distrivutions (61 ) eee 


statistics with three Jdeei male face  —accurae . 


section II addresses workspace design issues, © 
include: workspace requirements and assumptions 
regarding Sires use section III discusses the methods 


used to assess the accuracy of different asymptotic 
ADP rOxi mations, and the sample "“SiZeés"requlired tor “aa 
approximation to yield three decimal place accuracy. 
Section IV gives background information and discusses 
programming methodology for nonparametric tests based 
on single and paired sample data. [nee SCCELOM rey 


nonparametric tests for two or more independen (sean. 


are considered. section VI discusses nonparametric 
tests for association; and, Section VII déalsa wig 
nonparametric simple lincanwrecress1 on - Ssectiem VEE 


recommends other nonparametric tests that may be added 
to the workspace and areas for further work= 

To show application of nonparametric statistical 
methods to Department of Defense problems, several 


military examples are worked in Appendix C. 


LO 


PwWVORKSEACImDESLENS LSSUES 


This section presents a brief overview of the 
design Seems teerattons —~Used jim developing tne APL 


workspace for both the mainframe and microcomputer. 


ae BQOUIPMENT AND SOFTWARE REQUIREMENTS 

The mECreocompuccr  MUsSst be an IBM PC or  sAT 
Semeeaaltole, equipped with Si2 kilobytes of RAM and the 
Peewee bus; PC system software, release 3.9 or later, and 
IBM’s DOS, version 2.00 or later.! The 8087 math 
Peeocessor chip is not required two run this software, 


but will increase the computational speed. 


B. KNOWLEDGE LEVEL OF THE USER 


iiacome user iS expected to nave had Some exposure to 


APL and a working knowledge OL nonparametric 
SeeolSsctics. Pitiiari ty owlth MmlLeCcreoceomputers or the 
Naval Postgraduate School mainframe computer is 
assumed. 


PPneeArPL system software requires 144 kilobytes of 
RAM while the NONPAR workspace requires an additional 
190 kilobytes. 


i 


Ge SHLDSCTLION OPES es 

The nonparametric tests chosen for this workspace 
are some of the more widely known, and are considered 
basic material for any nonparametric Statist te] coun. 
More information about the tests can be found in any of 


the textbooks that are referred to in this document. 


O 
Zz 


ENU DIiSPhAvs 

The microcomputer’s workspace is designed around 
the use of menus. This was accomplished using the 
software vyvackage PC TOOLS from STSC. These menus are 


designed to guide the user through the selection of the 


tests without an excessive amounz Of Wer onp es The 
main menu displays (the (choices) availa) io the 
workspace, while the test menus give the background 
information and options available for each test. Help 


menus to provide additional information about the 


tests are also available. 


ie ORGANIZATION OF WORKSPACE DOCUMENTATION 

separate documentation is included Heels the 
microcomputer’s and mainframe computer’s workspaces 
(see Appendices A and Ee respectively). These 
appendices explain the organizatiton anew ore. 6 0) en 
the workspaces. Appendix C, which provides example 


problems for each nonparametric test, is applicable to 


both workspaces. 


Led 


DEE GENERAL SAMPLE SIZe CONSIDERATIONS AND ASYMPTOTIC 


APPROXIMATIONS 

Hue cots thesis’, the term alpha value is used in 
feeeeeeral sense, and refers to the probability of 
Meeeees ts a true null hypothesis. Mots Anishe nl Jere DiS 


me@ers tO she probability that a test statistic will 


1 
* 


exceed (or not exceed in the lower-tailed test) the 


Semeilbea Value, when tne hypothesis being tested is 
true. 

Fr OF selected values, the exact cumulative 
Met. OULLON =Unctions (C.D.F.) of the test statistics 


are comparecd with those obtained from normal based 
Ss9MocvoOctic approximations. The results of the compari- 
sons are used as a basis for assessing the accuracy of 
wre approximations. In those cases where more than one 
Sevy@M@oLOtic approximation has been suggested in the 
me beracvure, the accuracy of each approximation is 
compared over arrange of desired C.D.F. values’ and 
sample sizes. PeotietiecmrtecastlEs, the most consistently 
accurate approximation, and the sample size for which 
eac approximation ase Ae Postma bpirese decimal 
place accuracy is determined. 

Once the accuracy comparisons were completed for a 


specific nonparametric test, microcomputer capabilities 


Ve 


were “Cconstacredn. In some cases, generation of the 
exact distribution up to the €esired samptc ss 42Z-ueree. 
too long or was NOt possivlec on tee wee When this 
Occuresa, the mainframe computer was used to generate 
the required distrioutions with the results stored in 
numerical matrices Lior GQuicls recall DY the 


Nonparametric spr Ol ams. 
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IV .THSTS FOR LOCATION BASED ON SINGLE AND PATIRED- 
SAMPLE DATA 


The Tests assume that the data consists Qn a 


Ft) 


Merwe set Ot independent observations Xj; or paired 


Meret ons (X5 77), i-1,2,...,N, from a continuous 


rey 


Seser tOutTLOn. or the 


i oc ods al Le seCasSe.. 


) 


M@emonwteht hypotheses are concerned with the median of 
Seem cee na —2ne median of tne differences. Xi- Y, 
measoeCuLVel7T. meee ss CONS VectoC@eare tne Ordinary 


Sign Test and the Wilcoxon Signed-Rank Test. 


A. ORDINARY SIGN TEST 

The Siew test can be used to test various 
hypothesis about the population median (or the median 
or the population at differences). Confidence 
Maem vea lem tOr nese parameters can also be constructed. 
cme lia! OpLren,  nonparamnetrire contirdence intervals 
mere ene quantiles of a™ continuous distribution are 
oem ine Cd . 

7. Comoutaulon Of Sider Test Statistic 

Reme Sime lsc—-sample data, mre wGese Statviscic A 

SSD ugea a5 ute Nuniter of = observations Aj, greater 
Pimretoewaypotiesized median Mg. For the paired-sample 
Certo ims ttemmumMbper Of differences X;- Yj that exceed 


Mg. All observations Xyj (or X4- Yy) that are equal to 


SB) 


Mg are ignored and tte  sampte size decreased 
aceor dines, . As long as the number Cfeeucn 32. > eo 
small relative to the size of the sample, the test 
results are not greatiy affected. Gibbons [Ref. 2:pp. 
198]. | 
2 The Null and Asymptotic Distrimimommorwe< 

The WUutl Gs e@awoution “ore: is binomial 
wollte D = gfe In Table 1, the exact values of the 
C.D.F are compared with tne Correspondin2 weber etme ne 
values using a normal approximation With angie eeu. 


COnNtCINULTY “CORmCce ton. 


TABLES 1. ©.D.). COMPARTSONS FOR TERR Vslae se 


PRUBLX £ ki; PUR SeMPLS esi ee EOURL T1024, 
TEST SiAt, vebun | 3 | 6 (| , 4 8 | oe 10 a afl 
EXACT C.D. 8, | 00332 | .01133 1 ROgtSc tr 2 o7o19 1oniads 2 eee oe eee 
ERROR; NORMAL | .OO0117 | .00417 1 7Ot134 | 702456 1204339 1 ec2s2 eens 
ERROR; NORM. W/CO 1.00063 1°,.00104 1°.00114 17.00073 | .00001 | .90043 1 .00023 
re ree 


PROBLX £ k]; FOR SAMPLE SIZE EQUAL 10 2s. 


TES te o7n teehee onal 6 | ? 3 | : 10 ne 42 

! ! 1 ! ! 
EAnG1 ea | .00732 | .02164 | .05398 | .11476 | .24213 | .34502 | .50000 
TREOR) eO eeL | »V0deo.! soem ! MOL eRe | 1.93400 | oUsatZ ! ert. | i aes 
e2oGR; HUeM. W/O | .0008a | 00LiT g0u92 | Ce ieee | 200044 | 00000 


rrr eter ee peor Ne 
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AS can be seen, for sample sizes greater than 
i eetiom@mal Approximation withecontinuity “correction 
Hs eeecurate to at least three decimal places. 
wr Hypothesis Testing 
P-values At euecOnpuwead for three basic 
mppeOotnaeses comparing the median of the population or 
Cmem Mecianeor the population of differences M with a 
oypothesized median Mg. P-values are taken from the 
cumulative Gqdistrissutiene@eof the binomial for the 
Lowes CeSts of hypothesis. 
a. One-sided Tests 
(7) HO: M =Mg Versus Hi: M< Mg. The 
Pevelue Equals Pr{K < kK], where k is the -computed 
maeue Of the test statistic. 
(2) H®: M = Mg Versus H1: M > Mg. The 
P-value equals Pr[K > k]. 
b. Two-sided Test. 
Ged? Selene VersuceeHl: M 4 Mg. The 
P-value equals twice the smaller value of a(1) or 
a(2), but does not exceed the value one. | 
For sample sizes greater than 25, a normal 
exe ON With CcCOntinuizy correction is used. 
fe. Ponsa LoeneGecs ANvervaly aStamation 
Confidence iio yas Ol a Enc pepulat Lon 
median are based on the ordered observations in the 


Sample. For paired-sample data, confidence bounds are 


le? 


obtained from the ordered differences of the pairs of 


data. A 100(1- a )% confidence interval is determined 
in the following manner. Let k be the number such that 
Prfk < kil < (Sa2er Then, the (k+1)th and (N-k)th 


order statistics constitute the end @geamcs (ot Fine 
confidence interval. Gibbons [{Ref. 1: pp.104]. 

HOw Compuming confidenese  “2ngervars when 
sample size N is greater than 2955 a normal 
approximation with continuity correction is used. 

Also included” under thisS tést is ante ems 
generate nonparametric confidence intervals for any 
specified quantile given a sample size N from a 
contimeews. distribu; tede The end points of the 


intervals aré sample order srageecres. 


B. WILCOXON SIGNED-RANK TEST 

The signed-rank test requires the added assumpc ie) 
that the underlying distribution is (Ss ime rn ic- eS 
test uses the ranks of the differences Xj- Mg (or Xj4- 
Yi- Mg) together with the signs of these differences to 
determine the test Statistic. Cont taence wim te mie eee ers 
the median can also beweenstrietea- 

1 Computation of the Testeotatiwc 
For single-samplée data, thew geo es Cott su em, 


1s computed as follows. 
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| 1 aif X4y- Mg > @ 


{| ® if Xi- Mg < @ 
N 

and let ry = rank(|Xi- Mg!/). Then, W= ) Ziri. 
j=] 


For paired sample data, Wis calculated in 


the same manner, except the differences to be ranked 
are the paired-differences minus the hypothesized 
median. Zero differences are ignored and the sample 


size is decreased accordingly. When ties occur between 
ranks, the average value of the ranks involved are 
asSigned to the tied positions. It has been shown that 
a moderate number of ties and Zero differences has 
lgttie effect on the test results.! 
ae, These Sande AswmipmoeticaDistribudden of W 
Pacweraerenull disperipburion of Weis given by: 
Pr{wWw = w] = uy(w)/2N, w= 0,1,2,..-.,N(N+1)/2, where 
uyn(w) is the number of ways to assign plus and minus 


Signs to the first N integers such that the sum of the 


positive integers equals w. It can be shown (see 
Gibbons [Ref. 1:pp.- 112]) that un(w), for successive 
values of N, can be computed using the recursive 
relationship: 


uy Cw) = UN—1(w-N ) a ne Cw) 


IlFor more information on the effects that zeros and 
tied ranks have on the Wilcoxon Signed-Rank Test, see 
Pratt and Gibbons [Ref. 3]. 


19 


Exact C.D.F. values were compared with those 
obtained using the following asymtotic approximations: 
student’s T with (N-1) degrees of freedom (Le 
Student’s T with continuity correction (TC), normal 
(Z), normal with a continuity “comrect veda 4e). (aaa 
average of T and Z as suggested by Iman [Ref. 4], and 
the average of TC and ZC. 

As can be seen in Table 2 below, the average 
of TC and ZC gives the most comsiStemwmmy accurate 
results with three decimal place AaAcecuUracye Woes 
sample size exceeds 9. 


TAS aes C.D.F. COMPARISONS FOR THE WILCOXON 
SIGNED-RANK TEST 


PROBLW { w); FOR SAMPLE SIZE EQUAL TO 9. 
TEST STAT. VALUE | 3 | Sea 6 | 8 | 9 | 2s 
EXACT Cy Deen. | ,00977 | .04953 | .02734 | .04383 | .06445 | .42500 | .47969- 
ERROR; NORMAL [7.00067 | .00046 | .00204 | .00594 | .00958 | .041824 | .02272 
ERROR; NORM. W/CO 1 7.00243 17.00247 17.00167 | .00023 | .00269 | .00693 | .003806 
$$ eee 


ERPORs sf Dist { .00318 | .00608 | .00673 | .90701 | .00817 | .00826 | .00318 
ERRORS hau ce 1 .00O355 | .00276 1 2.00233 1! .00012 | T7O0G10) | s.00ts7 al ee 7ce 
LZKUR; AVE T/Z { .00225 | .00327 | .00433 | .O0e45 | .00988 1 .01325 | .04549 
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PASBILE 2. CContinueds 


Peootues WiyerOk onMroe olzk EQUAL TO 10. 


, T i ? rt => | ce esse 
Bor sini. VALUE | oT) (nn 8 | 20° | 2 | Lo | 1? 


~Z 


oor te 

E2202; NORMAL Lets j 2OUOds | wOCGI7 1.004575 1 -.0Ud37 I .0t490 | 2Otece 
’ : ! 

Peer; tenn. WCC 1 1.00270 17.002T9 17.0058 | .00043 | .00229 | .00558 | .00710 








SPROR; TDIST 1 .00399 1 .COS12 | .00525 | .00665 1 .00664 1 1.00664 1 .00684 
ERROR; TWCC —«s|-: «00249 1 .00243 1.00492 | .00451 1 7.00052 17.00376 17.0057! 
ERROR: QUE T/Z | .00442 1 .00267 | .00342 1.00574 | .00749 1 04075 1.01284 
ERROR; AVE TE/ZC | 7.00040 | .00012 1700008 1 .00097 | .00089 1 .00094 | .99070 


} ' t 


o bypOLNeSsis. Testing 
P~values are computed for three basic 
hypotheses comparing the median of the population or 
the median of the population of differences M with a 
hypothesized median Mg as shown below. 
a. One~sided Tests 
(1) H@®: M = Mg Versus H1: M < Mg. The 
P-value equals Pr[{W < w], where wis the computed 
Veue Of the test statistic W. 
G25) HOe M -9My Versus ™ H1: M > Mg. The 
P-value equals Pr{W > wl. 


ay Two-sided Test 


(1) H@O: M Mg Versus H1: M # Mg. The 
P-value equals twice the smaller value of a(1) or 


a(2), but not exceeding the value one. 


ou 


For sample sizes greater than 9, an average 
of the normal and student’s T approximations, cach with 
continuity correction, is Used. Computation = oie a. 
value for each@ort the altermative feet eee. 

a. Hee < Nig 


Let Pgc = Pr[ Z < (w +.5 - fy;)/ oy] and 


let Ppc = Poe [Eee 
Now2 [lw - owl =] »912 |e? 
N - 1 
where 4 is standard normaly T(N_-1))24se42 Seteen tae 
distri om Gerd with (N-1) degrees of freedom, uy = 
N(N+1)/4 and dwS = (N(N+1)(2N+1)/24)). Then, the 


P-value for the test is (Pzc + (1 - Poc))/2 if w is 
less than Uy and (Pac + Prc)/2, otherwise. The above 
formulas are obtained from those given by Iman “~Rer oa 
after inclusion of a GOntinurty  conprecruaen 

b. Hil: M SMe 


The P-value equals ((1 - Pac) + Prc)/2 if 


w is less than Uy and ((1 - Pac) + (1 - Pre))/2, 
otherwise. The computatwion es Pye ance e 2s camila. 
to the above except” "he 0s Seneca. SOmernul Va 


The P-value equals twice the smaller 


value of aor b above, but not exceeding the value one. 


ZZ 


4. Contidence Interval Estimation 

er Single-sample data, the eonridence 
interval for the population median is based on the 
ordered averages of all pairs of observations (Xj+ 
Meee such that i < j. A 100(1-a@)% confidence 
interval is determined in the following manner. Let w 
be the number such that Pr({W < w] <« (a@/2). Then, the 
(w+1)th and (m-w)th order statistics, where m = 
NCN+1)/2 or the total number of paired-averages, 
Sememeeoe Cie end points of the confidence interval. 
A confidence interval for paired-sample data is 
computed in tne same manner, except the end points are 
taken from the paired-averages of the differences Xj- 
Yy;. Gibbons [Ref. 1:pp. 114-118]. 

For computing confidence intervals when 
Sample sizes are greater than 9, a normal approximation 


MeaeecOMNCINuLlty correction 1s used. 


ae 


V.TESTS BASED VON TWOSOR MORE SAM vings 


The. tests sassume Ghat the davameon sis: On, 
independent random samples from two or more continuous 
Gistreaout i onmse The general null hypothesis is that the 
samples are drawn from identical populations. The 


Mann-Whitney and Kruskal-Wallis tests are considered. 


ae MANN-WHITNEY TEST 


The Mann-Whitney test is based on the distribution 


of the test statistic U, which can be used to compare 
the equality of tne population medians Or variances 
for two samples.'! The Mann-whitney test with a 


modified ranking scheme can be used to test (iar 
equality of variances if the population means or 
medians are assumed to be equal (Conover Len D> DDir 
229-230]). If the medians differ by a knaown —amounteeone 
data can be adjusted before applying the test. A 
confidence interval for the difference in the medians 
of the two populations can also be estimated. 
Ve Computation of the Test Statistic 
For the comparison of population medians, the 


test statistic U is computed from the combined ordered 


IThe test statistic U and the method used to 
compute it are taken from Gibbons [Ref. 1:pp. 140-141]. 
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arrangement of observations X4 and Yj, ss Pe een gts 
ae ee ks Ml Het Urq — fank (x, ) in the combined 


N 
Sm@g@ered sample and Ry = ere: ine Tae, 
= 
U = Ry - M(M+1)/2. 


For testing the equality of variances, ipals 


computation of U is similar except for the method of 


assigning ranks to the ordered sample. This method 
ranks the smallest value 1, largest value 2, second 
eamsest Value 35, SeCOme Stat LesG eo tuc 4, and sono 


Eero Ss, “Until the middle of the combined ordered 
sample is reached. For either test, tied ranks for the 
combined sample are assigned the average value of the 
ranks involved. Paiodctateo we mmincecimOr tes Nas little 
Etm@ect On the test results 

Be The Null and Asymptotic Distribution of U 

The exact Meet Sin Ou om) Of U is 
determined using a recursion algorithm due to Harding 
Pet need |. 

Exact Crepe F. values were compared with 
approximate values obtained from the following 
asymtotic distributions: student’s T with (n-2) degrees 
Sreurecdom swhere nm =N +M, the total number of 
observations in both samples (T), student’s T with 
pomemdity cormection (TC), normal (Z), normal with 


continuity correction (ZC), the average of T and Z 
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(Iman [Ref. 7]), and the averacetopeGeoncee ee The 


results for various sample sizes are yztven An fenton 


To C.D.F. COMPARISONS FOR THE MANN-WHittiia 
TEST 


PROBCU { ul; PUR TSRMPEE GiZzs N EQUAL TG> 2 nib eeOua comes. 








TEST STAT. VALUE | $4 4 47 1 48 § 24 1 2004 27 1 29 
SXAGT Cadet | 00938 | 01999 | .02543 | .04696 | .06796 1 112304 | .17005 
TRROR; NORMAL «7.00028 | .00100 | .00163 | .00444 | .00682 | .04243 | .01514 
TRROR; NORM. WCC 1700146 17.00444 1 7.00083 | 00028 | .00130 1 .00354 | .00436 
ERROR; TOIST | 00237 | .00937 | .00372 1 00484 1 00525 | .00876 | .00773 
TRROR; TUCO  —«t-«. 00449 1.00408 | .00035 1 .00007 1 7.00078 1 7.00259 1 7.00328 


‘ 


EO 2 esa eee f .OO10S | .00229 | .00279 | .00453 | .00603 | .00957 | mits 


, 


° 

* 

e 
7 ~}— 
— 


L2R0R; AVE TC/ZC $7, 00014 9) = VOCUS er CCCG ta mmc ecud | .00025 | .00048 | .900S4 


csr csi ft fe fp le 


PROBCU £ ul; FOR SAMPLE SIZES N EQUAL TO 7 AND M EQUAL TO 12. 


en eee a, cc ee a 
TES? Sipe plus i4 | fie i 2 ai o4 of | 30 


SS a a a 


EXACT C.D.F. | .00853 | .01792 | .o27s25 se Csise Ae Ori | raiase feo 


See ee 


ERROR; NORMAL (7.00045 | .00062 | .00187 | .00359 | .00701 | .01097 | .01437 


—— a ee 


ERROR; NORM. W/CC 17.00452 17.00128 17,.00079 17.00003 | .00154 | .00322 | .004339 


ERROR; T DIST | .00199 | .00292 | .00356 | .00422 | .00527 | .00643 |! .00796 
ERROR; T W/CC Pe OOUga a ,00092 | ,00068 | ,00026 17.00066 1 7.00178 1 °.00252 
LsrUR AYE 1/2 | .OO0077 | .O0O177 15.0027 se O0s [i ee or mec 1 | Ot1se 
Sae0R: AVE Te/To 1.90029 Wo. J00da"s . 0000s Tt 000 I] a0 = ge ee eee 


' } } f 
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TABLE 3. (Continued ) 





Decoy osha. VALUE | oe | 15 3 al 20 | os | oi | 30 

es | Gets 2 Pees | Ute Osmiee ome balntoto | 070ml «izes | .17335 
SS SE 
2205; ORHAL eee cOU rie! . Ceese Ooo 1. 00e49 | «09639 | .Of2t) | .901564 


t ! 


at wan p i 


Pee eet. “3/CC 1.90190 | .00TSO 1 °.00087 | .cOOO7 | .00193 | .00444 | .00573 


ean0s: F OST PM SOI OO coo toes! 200890 | .00Sa4+ 1 1.00774 1 .003S5 
—_— 00 oov>q$q5498 0 9 99 )3ssSs£(n4g 
Peron? T W/CC re COs logs sl 1002s O00 er 1) pO000? | .000235.1 4.90051 
a 
Mees; Ave. t/t eC Uae Jo ee er Uae seie.cOss? | .006is | .00992 | .O12S0 


i 
: gs 
4 2 
oy o 
¢ 
we =e 
> 
4 
Cc 
© 
é 
< 


meee; HYE TCC Sse. 00 Cemee JO0cegk .UCUls | .COO37 | .00209 | .00264 


@ 
t ' t l ! ! ' 


Beer) Sui: POReSenPle Si Zs 4 SQUAD TOS AND M EQUAL TO 27. 


a 





meteor. VALUE | fo 10 | a co | 2 | 23 | 26 


—————————————E————E————EeE——————————————————— 
BACT C.D.F. | .00764 | .01650 | .02512 | .04296 1 .07734 1 .12660 | .17433 
a Se ener ee 
BR2O2; NORMAL 1-,00265 17.00099 | .00072 | .00389 | .00874 1 .04342 | .01680 


1 $ } t 


ER20R; NORM. U/CC 1 7.00363 | ~.00254 1°.00133 | .00089 | .00405 [| .00663 | .003831 


rr ne 
Oreo 


- 
q 


BeeCR; T DIST [~.00124 | .00034 | .00173 | .00414 | .00741 | .01026 | 

Enaon; + W/CC f~.00219 1°.00129 1°.00042 {| .00097 | .00250 | .00323 | .00330 
I ; 

ERROR; AVE T/Z [~,004133 17.00034 | .00122 | .00402 | .00803 | .04184 | .91466 


, : 1 7 1 1 


eee eee | 002g | Cogis2 | .000g7 | .00093 | .00327 | .00496 | .00619 
Nee ee eee eee ee ee eS el 

Nome enmeese Seen from themtaprves. C2eweeer er ace 
“eco emda tC 2£ives Ghe most consistently accurate 
Gesultcs. For sample sizes NxM > 8@, nearly three 


Geecumeam place accuracy is obtained in all cases. 


ey 


Da Hypothesis Testing 
P-values are computed for three basic 
hypotheses comparing the medians or variances of the 
two populations as shown below. 
a. ONe@w cho cemlecaus 
(1) H@: My = My Versus Hi: My < My or 
H@: Vy = Vy Versus Hi: Vy > Vvo [heer ool ue secular 


Pr[{u uj], where u is the observed value of the test 


JA 


statistic. 
(2) H®@: My = My Versus Hi: My > My or 
H@: Vy = Vy Versus Hi: Vy < Vy-o ihe P-vyeiee sequaike 
Pe > ule 
De bVO-o dean tee G 
(1) H®@: My = My Versus Hi: My # My or 
H@: Vy = Vy Versus Hi: Vy # Vy. The P-value equals 
twice the smaller value of a(1) or a(2), but not 
exceeding the value one. 
For sample sizes NxM greater than 8@, the 
average of the normal and student’s T approximations, 
each with continuity correction, is uscea. Computations 


of the P-value for each alternative hypothesis are: 
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cope: EvieeGeion or Wy >” Vaz 





cover gee- meee Cae oe iy) fo, | and 
Bemopo = Pri T(n-2).s ole eae 
(N+M-1) a," [ju-wyl-.312 |? 
— a2 ~~  NeMG2 
peewee 4 is standard normal, T(n-2) has a student’s T 
distribution with (n-2) degrees of freedom, Uy = NxM/2 
and 4° = (N(M)(N+M4+1))/12. Then the P-value for the 
test is (Pzco + (i - Ppc))/2 for u less than uy, and 
(Pro + Prc)/2, otherwise. The above formulas are 
Gerearned from those given by Iman [Ref. 7] after 


Piewuslon ot the continuity correction. 
Db. Hi: My > My or Vy < Vy 


The P-value equals ((1 - Pac) + Prc)/2 if 


u is less than wy and ((1 - Pac) + (1 =- Poo))/2, 
otherwise. The computation of Pvc and Pyro is similar 
co @iewe a0Ove Except the sign of the continuity 


Corprection is changed. 
c. H1: My # My or Vy # Vy 
Laer P-value. egal saatwice. the smaller 
Value of aor b above, but not exceeding the value one. 
4. pom Jdeaecs A kerva som maulon 
Confidence intervals for the difference in 
medians, (My —- My), are based on the ordered 


arrangement of the differences (Yj- Xi), j = 1,2,...,M; 


Zo 


i= (1.2,.44a for »ali 2s gaacee 5 A 1@0(1- @)¢ 


confidence interval is determined in the following 


manner. Let ube the) nNUmber stich that Sree 
( Qj aoe Then, the (u+1)th and (m-u)th order 
Statistics: where m = NxM or the total number of 


possible differences, constitute the end points of the 
confidence interval. 

For computing confidence intervals when 
sample sizes NxM are greater than 8@, a normal 


approximation with continuity “Correction is Mecae- 


ioe KRUSKAL-WALLIS TEST 

The Kruskal-Wallis test is a nonparametric analog 
of the one-way classification analysis of variance test 
al equality of several population medians. Gibbons 
[Ref. 1:pp. 199]. 

[ Computation of the Test Statistic 

Calculations of the test statistic HS cemie, 

around the ordered arrangement of the combined samples 


from which the sum of ranks for each sample is derived. 


Let Xij> j=1,2,.6.,Ny and 1-15 25 oy Ky ree ele nan 
random samples from k populations. Let rij = rank(Xjij), 
Ni K 
Ry = Des Ei eae aati Then, 
j= i= 


k 
ee C2 /08061))] >. (RY2/nqz)}—- 3(N41) 
i= | 
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Iz ties occur in the combined sample, they 
are resolved by assigning the average value of the 
ranks involved. A correction based on the number of 
observations tied at a given rank and the number of 
memes involved, is included in the calculations. A 
Semplete description of the correction factor is given 
in Gibbons [Ref. 2:pp. 178-179]. 

te. moe Nita easymptOuLC distribution of H 

Tne mnuli distribution of H is generated by 
enumeration. Sach. DesssmomeDermutaction of ranks is 
listed for the combined sample, and the corresponding H 
[mee computed. The frequency distribution of H.,is the 
total mumber of occurrences of each distinct H value. 
The H values are arranged in increasing order while 
maintaining the frequency pairings. The null 


Piet ouUutlOn is Obtained by dividing the cumulative 


frequencies by nj!no!...ny!/N!. 

Due to computer limitations, generation of 
the exact distribution of H was only possible for k = S 
populations with n = 4 observations in each, and 4 
populations with 3 observations in each. Most of the 


distributions wers generated on the mainframe computer 
and saved in matrices for quick recall by the Kruskal- 
Wallis test program. 

Exact C.D.F. values were compared with the 


corresponding approximate values using the following 


Bib 


distributions: chi-square with (k-1) degrees of freedom 
(C), F distribution with (k-1) and (N-k) degrees of 
freedom (F), and F with (k-1) and (N-k-1) degrees of 
freedom (F1). The chi-square distribution uses’ the 


Kruskal-Wallis) .heestets soe while Cie eee waa cl F' 1 


distributions use a modified H SiieuerSaeC.,. cette 
((N-k)H)/(k-1)((N-1)-H); see Iman and Davenport [Ref. 
Bl. As can be seen in Table 4, F1 gives the most 
consistently accurate estimates. 


TABLE 4. C.D.F. COMPARISONS” PUR The 23 
WALLIS TEST 


PROBCH 2h1: FOR 1 CaAOUP OF 3 SAMPLES CONSISTING OF %, 4, AND 3 OBS. 


TEST STAT. VALUE | 7.08439 | Sandan | 6.1828 0Meeesoes Neo sad. || Goce eee ae 
skacr ¢.0.2, [7.90870 [192908 | .02961 | .owaes [107810 | :a2sia [i778 
SRROR; CaISQUARE |~ .01840 |~.01982 | .01585 [9.012205 001lsu |) (007 o meer o 221 
zeaon; F Dist | 00304 | 00736 | .o0aa6 | ‘01113 | ‘01820 | -o16e6 | .009a0 
ERROR: F W/"1 DF | .00084 | .00426 | .00403 | .00544 | .01138 | .00895 | .00172 


PROBCH 25]; FOR A GROUP OF 3 SAMPLES CONSISTING OF 4, 4%, AND 4 OBS. 


TEST STAT. VALUE | 7.6538 | 6.9615 | 6.5000 | 5.6923 | 4.9615 | 4.2692 | 3.5769 
oe |_,00762 | .01939 | .o2996 | .ougés | .08000 | .12190 | .17299 
ERROR; CHISQUARE |" 01416 |" .01139 |" .00882 |" .oo941 |" .00368 | .00361 | .00577 
ERROR; F DIST 1.90290 | .0083a | .01204 | .01100 | .02272 | .02225 | .00263 
ERROR: FP W/'1 DP | .00149 | .00600 | .00890 | .00647 | .00708 | .00592 | .00384 


eae eee tae O88 BF FF esBZFerweweweswr_eswesweesreewere B@erewe_ere_esw_eireswetresweeswedereesweswee@wetresrereeeswereeereae Beet BF@e ewes ewe st@ew wBFewteeeae & & = = 


PROBCH 2h]; FOR A GROUP OF 4 SAMPLES CONSISTING OF 3, 3, 2, AND 2 OBS. 


PEST STAD. VaLUBola7 536% | 7-181 gui\ 00 oe ee 
EXACT C2PeF ean b 201000 | 20292) | 02828 Tee ae 
ERROR! CHISQUARE ||" .OUM16 Ve OUT L271 eo re eee oe eee oe ee 
EREOR; 2 Ohor | .9028% | .90260 | .00730 | .00880 | .010S54% | .01123 | .9C0902 
$2202; P d/71 5F [7100172 7.30udu |7,302"7 17.30307 17130378 17,30870 |7.a08ua 


sme eee Zee se ewe Fes 2 2 2 Se Se 2s 2s ew eee ee ewes ewesFr ew eweweee es Fee Fees Fees ese se es ese ese eee ee mew eewwtrewe Be Be Se | = = 


PROBCH 2hI; FOR A GROUP OF 4 SAMPLES CONSISTING OF 3, 3, 3, AND 2 O&S. 


2S wer eres eee es Oeste eseee@ereese sew st ete @eeewewew ew eww eFZ BFF ee ZTeFZ BT eeZTZ wseweww ws esFT aeseBVererw ew wFZ BF BBW weweBZFet wee F ws F282 8 2 | = = 


TES? STAT. VALOE | 8.0152 | 7.926% | 7.1549) | Oh ee eee 
OO {00962 | 02832 | 102074 Oot See ee 
ERROR? CHISQUAREN | .03609 7203564 5.02 eee eee 
BRROR; F DIST ___{_-00215 | .00uga | coud | 00820 tate esse eee 
ERROR: F W/'1 DF |~.00133 |]~.00019 |~.00269 | .00030 | .00098 |~.00230 | .00091 


A final accuracy comparison between the C and Ft 
approximations was conducted by computer simulation for 
SePpopulations with 8 observations each. Neale all any, 
39,9990 permutations of the 49 ranks were randomly 
generated (no tie ranks allowed) and the H statistic 
@eeculated for each permutation. Pen pane  emogricall 7 
Getermined percentiles Hy for selected values of op 
petween | and -18 were compared with the 
SaeonextMations given by the C and Fl distributions. 
Steer esulLts are shown in Table 5. Pecan be Seen thar 
the F1 approximation compares well with the simulated 
rPesuits, giving three decimal place accuracy, while 
tee G approximation is less accurate. 


beeen 86D C.D.F. COMPARISONS FOR THE KRUSKAL- 
WALLIS TEST USING COMPUTER SIMULATION 


PROBCH 2hijl: BASED ON 10000 GENERATED H'S FOR 5 SAMPLES OF 8 OBS. EACH. 
Tes? STAT. VALUE | 12.229 | 11.065 | 10.24a | 9.212 | 9.129 || 7.030 | 6.232 
C.0.F. VALUE {_,01000 | 02000 | .03000 | .05000 | .o8000 | .13000 | .18000 
er Azz) 00a74 Focuses 17 006%6 | 00601 | .o0ed6 | .co¥a2 | .902K6 
Memon: F ¥V/ 1 DF | .00081 | .00231 | .00259 | .00350 | .00161 | .00109 | .00103 


PROBCH 24]; BASED ON 20000 GENERATED #'S FOR 5 SAMPLES OF 8 OBS. EACH. 


= eweweee ws @ eee ewes ewewewewese FFT Bs swe ewe wZ FT FT FT FT FT FT FT BTSs TZ TZ ZT FI wss ZZ FFT SF FTF WOFseZ FF FFs ese FEZ FSF wTse swewrwws22 77 BF 2 


a ee oe ea es ele 2 89) ee 228 | 7088 | 82220 
eee ee. ee ee eee eo Ooo Seem C8000 4 213000 1 218009 
et eek ees eo a 2 O00 716 |esgosee | oours 1 290393% 
eee 47s OF | 1.00126 | .0902272 | .00166 | .00235 | .00161 | .00130 | .00199 


=e eeesweenwes S82 SB Beers ewese wBewswesesw es BO ewWsw ese BTEZewZwwZTwZT Ss OFTZ FTF TZ FZ T= ZT BOS FBO TSI ST FT Ss SF SF TE TFTFZ TF FT FZ Ze TST FTF FT Ze ws FF | @ 


Peoe@ee 225 34520 IN 20000 CENERATED d'*S FOR 3 SAMPLES 9F 3 OBS. Z2Acd. 


=aeae ea eaeeeegeaeqqeg@eeqge#ess =sS eaqeee e222 22 @2ewe Fe esse eaeeeseweae es FF eee we eZee sees eT esse ses eweweseess Fe es ee Ss ees SB we FTF FTF Fess es ses 2 es |] | 


aa eRe ee a 828 | 72072 _ | 6.265 
eee _ ee eee ome ecee| C0007 | 2-13000 | 18000 
ee ee) eee | 7 seu. | -G08o2 [7.00677 Wngoses 1700273 | .00020 
ERROR; ? W/"1 DF | .00121 | .00143 | .00111 | .00273 | .00301 | .co34u | .00142 
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De HyYDOUnes Ise becom 


P-values for the test H@: the population 


medians are all equal versus H1: at least two 
population medians are not equal, are computed as: 
PrtH > hj], where h is the value of the observed test 
Stabs iace 


For three or more populations with at least 4 


observations in each, the Fl approximation is used. 
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Pence OR PASC OC LATION EN. PATRED-SAMPLES 


wee SGesus cCescribed herein assume that the data 
S0Nsists of independent pairs of observations (X4,Y¥j) 
tieemm a Olvariate distribution. The general null 
hypothesis is that of no association between X and Y. 


Kendall’s B and Spearman’s R are considered. 


Ree KENDALL’S 8B 
le Computation of the Test Statistic 
mae SSUestessceclStic 1S compmured by comparing 
fecme observation (X;,¥;) with all other observations 
Sia 5 ) in the sample. If the changes in X and Y are 
of the same sign, sgn(X; =e) = sen(Y; ~ Yj), the pair 


(X4,Y;,) and Ce a) ia) "Concordan: ' mcdrea +1 is scored. 


nEeeeweec SipNnS are different, the pair is "discordant" 
and a -1 is scored. Any ties between either the X’s or 
mre ¥’s scores a zero for that pair. The sum of all 


scores divided by the total number of distinguishable 
pairs , (N(N-1))/2, gives B. If zeros are scored, the 
demominator is reduced by a correction factor which is 
pesca On ithe number of observations tied at a given 
rank and the number of ranks involved in each of the X 
and Y samples. econ aero aesen i piusOn Of the 


Gorrection for ties is given in Gibbons [Ref. 2:pp. 


2 


289]. The value of B ranges between 1, indicating 
perfect concordance, and -1, for perfrecy discoreamec- 
Gibbons [Ref. 1:pp. 209-225]. 
oe The Null and Asymptotic Distr2bution of 
. The null distribution of B°2sederived fam 


the following recursive formula given in Gibbons [Ref. 


1? pp 422 cle 

u(N+1,P) = u€(N,P) + u(N,P-1) + u(N,P-2) + ...+ u(N,P-N) 
where u(N,P denotes the number of P concordant 
pairings of N ranks. This formula is used to generate 


the frequency with which the possible values of P 
OCCUR Division “Dy eh. results in the probability 
distributvene |on oF. Since, B= (4P/(N(N-1)))-1, =the 
null distribution of B is easily determined. 

Exact C.D.F. values were compared with those 
obtained using a normal approximation, with and without 
a continuity correction factor (CC = 6/N(N2-1), pro- 
posed by Pittman [Ref. 11] for the Spearman’s R test). 
The results for various sample sizes are provided in 
Table 6. As can be seen, for sample sizes greater than 
|e. a normal approximation With conus Gy weer ween es 


provices three decimal place accuracy. 


oo 


feeseec. ) C.D... COMPARISONS FOR XENDALL’S B 


PWORLR 2 bisa sUR Shinbee slit ZQURL TO 13. 


Meeeeoiess YRLUS 1 0.5125 | O.8515 | 0.410391°90.3590 1 9.333371 0.2564 | 9.2308 





SAGTOC Der, | meO7t2 | .Ols2e4 |. 02366 | motego) | .Jotaeml 112593 1.13399 
———_—-_—_‘_--_--—-—- anreeecc — —=~$@osmmamssSssSsSsSsS3)8s SS oD MU 
E2208; NORMAL | .vO0Ols | .0OtsD 1.00313 | 700620 | .008038 | .Of473 >) .01703 


! ' ' ' t 4 
1 4 l 


meer seorn. WeC | 001s | "6007S | .00239 pmpo04s3? | .00658 | .01223 | .01445 


' ' t : ] ! 


Presta 2 bl; FOR SeHPLE SiZr fQUAL TO i+. 








Bemeeomsnee value | 0.4725 | 0.4296 | 0.4066 | 9.3626 | 9.2967 | 9.2527 | 9.2083 
’ 1 \ t ' ' ‘ 

‘Cs ee [eRtOeca nO nendenoeces | 03493 i 2O0cds3g | .tioso | «tase! 

eee eet one aa 

B220R; NORMAL P.2090gs | .00tHO | .90223 | .00431 1 .90889 | .O1235 | .OL6aT 
‘ t } ’ ’ ' ! 

Meeeeeenean. #700 [| .0CO0S | .0009S | .00152 | .00S45 | .00742 | .010R7 | .O2374 





De. Hypothesis Testing 


Cairo OnmerosicG 8Ol MG aSsoclaulOn De tween 
X and Y are computed for three types of alternative 
hypotheses. Because the distribution of B is symmetric, 
all probabilities can be taken from the upper tail 
using the absolute value of b, the observed value of 
mae Gest statistic. Linear interpolation is used when 
Bb lies between tabulated values. ieee aluies are 
Pemempeda as -OlLows. 
a. One-Sided Alternatives 
The one-sided alternative tested depends 


on the sign of b. A positive b will automatically test 


for direct association or concordance, wi le 4222s 
b will test for indirect association er dic cerdae 
The P-value equals Pr[{B > |bl]. 

b. Two-Sided Alternative 

The P-value equals twice the probability 

computed for the one-sided hypothesis. 

For sample sizes greater than 12, a normal 
approximation with conQinui vy .correct lone sae a. The 


approximate P-values are then: 


1 = Pr a < Geli CC) = 4,7 2 5 ees 
is standard normal, CC is the continuity coOrpeciueun 
pee = (0) and T= = (4N + 19)/SN(N-1), for the one- 


Sided test and twice this P=value for the two-sided 


test- 


Be SPEARMAN’S R 
The Spearman’s R Test requires the added 


assumption that the underlying bivariate ~distrfeutien 


is “GOntinueus: The test measures the degree of 
correspondence between rankings, instead of the actual 
variate values, and can be used aS a measure of 


association between X and Y. Gibbons [Ref. 1:pp. 226]. 
Computation of the Test Stage 
The test statistic  K iS Witeeneurea in Vrnie 
following manner. Let ry = rank(X;j) and sj = rank(Yj) 


and © Dj = Fj = ss.) lien: 
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N 
yn 


Rae | =i 
N(N@ - 1) 
where N is the size of the sample. Tei = ee Oe en 1) ak 
or Y, they are resolved by assigning the average value 
em coe ranks involved. MeCOmnecu. On fact Olgas Dasedwaon 


tne number of opdservations tied at a given rank and the 
number Ot iano et 7Ol Ved, omc Luded in nie 
ea culayvions. A complete description Of the 
MepeeecctlOon caccor is Ziven in Gibbons [Ref. 2:pp. 279]. 
The value of R ranges between 1, iteGmecating perfect 
merectu association, AgiGe i. for perfect indirect 
Mesociation. Gibbons (Ref. 1:pp. 226-235]. 
a. The Null and Asymptotic Distribution of R 

The null distribution of R for a given sample 
size N is generated by enumeration. The method, as 
presented in Kendall [Ref. 9], involves generation of 
an N by N array of all possible squared differences 
between any two paired ranks of X and “YY. All WN! 
permutations of N ranks are used to index values from 
the array. The sum of these indexed values for each 
permutation gives rise to N! sum of squared differences 
Piewmmat co bien Converted to the R statistic. The 
PoedicticvweGistribution Of KR is the total number of 


Securtele@cs Ot GCach distinct value of R divided by N!. 


Sg 


Due to mainframe computer memory limitations 
in the APL environment, generation of the distribution 
of R was limited to sample sizeemon)/ er less.  eeUeime 
tables, provided by Gibbons [Ref. 2:pp. 417-418] to 
supplement computer computations, a numerical matrix, 
called PMATSP, was created to store the cumulative 
distributions of R for sample sizes ess eee ‘Aes 
Wau i xX allows LTOr “quero. recall Ou, cumulative 
probabilities by the Spearman’s R Test program. 

Exact C.D.F. values were compared with those 


obtained using a student’s T approximation with (N-2) 


degrees of freedom (see Glasser and Winter [Ref. ota 
and a normal approximation. botn normal ang i 
approximations were computed with and without a 


continutity correction factor, CC = 6/N(N@-1) (Pittman 
[Ref. 11]). From the results presented in Table 7, the 
most consistently accurate approximation is given by 
the Tf distribution with ageconwrecemrom 
re Hypothesis Testing 

P-values for tests Of no associat on. bDetwecm 
X and Y can be computed for three types of alternative 
hypotheses. Because the disStrisution Of fies ssyumecrl tees 


all probabilities are taken from the upper tail using 


lThe memory capacity of the mainframe computer in 
the APL environment is limited to 2.5 megabytes. 
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TABLS 7. C.D.F. COMPARISONS FOR SPEARMAN’S R 


PomoeR 2) sl roOR sanPlessrze SQUAL Td 3. 


i ee SK Me a  ——  — —<—-- 1 - SCC 
fee estAr. “LUE 91 0.7833 | 0.7167 | 0.6507 | 9.6000 1 90.5333 | 0.4333 | 0.33500 


§ ' { 


ee. 


a= 
Or 
& 


poce C.D... f .00¢eLl | .01943 | .02944 | .04840 | .O7S75 3 12243 


t ' 
‘ 





‘QO 


ERACR; NORMAL terousi2 | ~oueo 1 .JO0es | .00Ss6 | .00805 | .cts77 1 .9137 


= 
- 
— 
- 
— 


EreORs PHORM U7eC | 6.00853 1°.00443 1°.00135 | .00123 1 .00493 | .01029 | .01296 








bee; + O1St ees | 90352 ! 200451 | .09459 | 00415 | ,00238 | .00138 
ERROR; T WCC | .00153 [| .00203 | .00251 | .00175 | .00042 1°.00212 | °.00479 


t i I t ! ’ { 


ceo i ieeenr 2 Silt EOURL 70m. 


’ ' 


Meme tacoevalJE | O.¢255 | O.ofar | 9.5254 i 0.3626 | 0.4909 | 9.4061 | 9.3333 
— A ee 
Beace@ Cod.?. | 300670 | s7ot2ee| .O2i2e | .04dit | 907744 | .42374 | 117437 
TERR; HORMAL ie em es 0-20) *.0C0e ln 0087. ) (0700.1, 9iaisa| .0Lsre2 


beRUR; WORM. W/CC 1.00457 1°.00327 1 °.00219 | .00095 | .00451 | .00857 | .O011423 


ee 
BESOR; T DIST Doo eta ooees) | .CO326 1.00323 | 1.00253 1 .00159 | .00107 


, } u 


E220R; T W/Ce } .00444 | .00435 | .00134 | .001414 1°.00035 1°.00230 1 °.00361 
a ne rn 


the absolute value of r, the observed value of the test 
Statistic. The P-values are computed as follows. 
a. One-Sided Alternatives 


The one-sided alternative tested depends 


See ote Si2n of r. PaO sieGmwewmr Siete test for direct 
association, ie ee zeautvower tests for indirect 


association. The P-value equals Pr[R > |r|]. 
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b. Two-Sided Alternative 
The P-value equals twice the probability 
computed for the one-sided hypothesis. 
Mor Sample sizes  werecater ~erean TGs an 


approximation based Om the student’s SR distrisurieca 


Wwiteny (N-2) degrees of freedom aime COmt i nme 
COFTrectl ong. s Useca- The P-values are: 
1 - Prf TCN2=2) < GC —- COR) ec oe 


where T(yN-2) denotes the 2 dis oni bu al orm eee licece 


degrees of freedon, CC is the continuity Correction: 
Ur = @, and ore = (1 -(|[r|]-CC)2)/(N-2), for the one- 
steed reste and twice this P-value for the two-sided 


test. Gibbons [Ref. pe. Zieie 
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VII.NONPARAMBTRIC SIMPLE LINEAR REGRESSION! 


Nonparametric Linear Regression assumes that the 
data consists of independent pairs of observations from 
PolLvariLatve Gistribution and that the a Gomessi on OL ~¥ 
on Mele) ba mear.. The program estimates linear 
regression parameters based on the data samples. eG 
then allows the user to input X values to predict the Y 
values. Hypothesis testing and confidence interval 
estimation for the slope of the regression equation is 
Gee cred. it tne estimated slope lies outside the 
conridence interval, an alternate regression equation 
is offered with an opportunity to input X values’ to 


predict the corresponding Y values. 


a. COMPUTATION OF THE ESTIMATED REGRESSION EQUATION 
The least squares method is used to estimate A and 
Pmt wle PTesxresslon equation Yj = A + BX; + ej 
fer-1;2,...N), where ej; (unobservable errors) are 
assumed to be independent and identically distributed. 


A and B are computed from the following equations: 


'Except for program design considerations, the 
information and concepts provided in the section are 
paraphrased from Conover [Ref. 12:pp. 263-271]. 
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B= =i i || joao 
N Z 
N >> X42 oa ( ra) 
| ie) 
N N 
wii Baa 
og 1=1 1S 
N 
Be HYPOTHESIS THSTING 


P<values for testing hypotheses about the slope of 
the regression equation are based on the sSpearman’s 


rank correlation coerriciléemit k Setween Tie eae 


Yy -—- Boag, where Bo is the WypeeeestZe0 Steve he 
appropriate one-sidéa test Gi ay potrtesis, "hv ene 
versus Hl: B< Bg Or Hi: B > Bf, °1S  aucemavre oe 


chosen based on the sign of the computed test statistic 
r (positive r tests, H1: B > Bg; negative r tests, H1: 
B <98o))= The P-value is computed as: Pr[R > [rl]. 
. P-values for two-sided tests ene. 5 —- epee hels oe 
Bg, are also presented. 

For sample sizes N greater than 190, P-values are 
approximated using a T distribution with (N-2) degrees 


or freedom and Continul vymeorrece ore 


C. CONFIDENCE INTERVAL ESTIMATION 
100(1- a )% confidence bounds “for ) tite slope 


parameter B are determined as follows. The n possible 
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Peemess e506 C= CY GH YG)/(X4-X5), are computed for all 
ees Of data (X%,,¥%;) and (%7,.Y;) such that i< j and 
Xj # Xj and rearranged in increasing order to give g(1) 
mene mee SK) Chet w be the (1- a/2) percentile 
meee dlstrloution or Kendall’s statistic with sample 
size n.! Let d be the largest integer less than or 
equal to (n-w)/2 and u the smallest integer greater 
than or equal to (nt+w)/2 + 1. Then S(d) and S(U) are 
the desired lower and upper confidence vounds , 
Pea pee pe vely . 

sce: sample sizes larger than 13, a normal 
me GeaematblOn Wwicn continuity correction is used to 
Pee lase the coniidence intervals. 

If the slope of the estimated regression equation 
does not lie within the computed confidence interval, 
moe soroOcram automatically calculates a new regression 
equation where the slope is the median of the two-point 
Slopes Sjij and the intercept is the difference of the 


medians of the X and Y samples, My- My.2 


'Kendall’s statistic is defined here as Ng - Ng, 
water ée Wels the rne@hiber of concordant pairs of 
Se@cerYautens and Nq is the number of discordant pairs. 
penover Wrer. i2:pp. 256]. 


2This procedure is recommended by Conover’ [REF. 
neem. 250 |. 
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Vili ARRAS SOR PUR BER Orme. 


To create amore versatile and powerful software 
package, the NONPAR workspace could be expanded to 


include some or all of the following nonparametric 


tests: tests for randomness based on runs, Chisquare 
and Kolmogorov-Smirnov(K-S) Goodness-of-fit tests, 


Chisquare and K-S general two sample distribution 
CESS. Chisquare test for independence, and the 


Friedman test for associawten 
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Pere aN DIT A 


DOCUMENTATION FOR THE MICROCOMPUTER WORKSPACE 


ie General Information 

This appendix describes the organization and 
operation of the IBM-PC (or compatible) version of the 
workspace. Appendix C continues from where ageless 
appendix leaves off, to walk the user through each test 
EyemvOorKing practical examples. 

Before proceeding any further, the user’ should 
Deazem tO section II (Workspace Design Issues) for 
2zeneral information about workspace requirements) and 
assumptions regarding its use. 

To get started, enter the APL environment in the 
usual manner and load the NONPAR workspace. 

Zot Workspace Menus 

This workspace is designed around the use of 

menus. They guide the user through the selection 


process of choosing a nonparametric test and a test 


ODEON. Three types of menus are used; the main menu, 
test menus, and help menus. 
a. The Main Menu 


Wat ian moments of loading the NONPAR 
workspace, the main menu will appear. fase cee led 


Nonparametric Statistical Tests. MPesoeienlum= presents 
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general information abcut the workspace. Ets primary 
purpose is to list the choices of nonparametric tesus 
available and provide an option which allows the user 


to exit the main menu into APL CeomeGpwdaca Jace ee 


workspace or return to DOS. Each test choice is listed 
with some information about the test’s area Oa 
appiilecartioen-. To make a selection from the menu, move 


the cursor (using the cursor keys) tom@hazolighiee tie 
desired choice, and press enter. As a reminder to the 
user. a footnote at the bottom of the screen describe. 
the procedure for entering a choice. Once a test has 
oeem selected from the main nem a SUD-menme 
appropriate’ to the test "appears: To exit from any menu 
back to the main menu, press the Escape key. 
b. Test Menus 


The title of the test menu is the name of the 


nonparametric test chosen. The text portion of the 
menu gives a general overview of the test, to include, 
the method used to compute the test statistic, and a 
description of the various “opttens, Sthat ‘may be 
exerciaca. The third section consists of the list oF 
test options available. These Options include 


returning to the main Menu Or CHOSS tie ae pe heme 
Test menus may have options listed in single or 
mubtiple-paged formats. The comment in the final block 


of the menu lets the user know if a certain menu is 
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multiple-paged oor not. To make a selection from a 


maitiple-paged menu, use the page-up or page-down key 


PepeiCecemtiomune GCosired Opt Lon... Proceed with the scroll 
mevys to Dighliignt the choice, ao. Presa en being Once a 
test option is entered, tie Ser hom pitom py wc —vememer 1D UG 
the data required to run the test. When the option for 
more inrtormation is selected, the help menu is 
displayed. 

on Help Menus 


The title of the help menu usually begins 
mete ole words “More Information About..." followed by 
the titie or the nonparametric test. Pace cexu DOrgion 
of the menu explains the test and its options in 
greater detail. No choices are offered in the menu. 


To return to the test menu, press any Key. 


3) IL 


AP PEE oS 


DOCUMENTATION FOR THE MAINFRAME COMPUTER WORKSPACE 


Te General Information 
This appendix describes the organization and 
operation of the mainframe computer workspace. To load 


a copy of the NONPAR workspace from the APL library, 
enter the APL environment and type: )LOAD 9 NONPAR. 
Within a few moments the variables LIST and DESCRIBE 
are displayed on the screen. These variables provide a 
description of the workspace. 
oe The NONPAR Worksvace 

The NONPAR workspace consists of seven 
programs which call several subprograms during their 
execution. The Seaeue syntax for each test and its 
corresponding nonparametric test name is given in the 
following formaw: 

SYNTAX: Nonparametric Test and Application. 


a. SIGN: Ordinary Sign Test for bocation in single 
and Paired-sample Data. 


b. WILCOX: Wilcoxon Signed-rank Test for Location 
in Single and Paired-sample Data. 


c. MANNWHIT: Mann-Whitney Test for Equal Medians or 
Variances in Two Independent Samples. 


dad. KRUSKAL: Kruskal-Wallis Test for Equal Medians in 
K Independent Samples. 


Dr 


e. KENDALL: Kendall’s B: Measure of Association for 
Paired-sample Data. 


rt) 


SPEARMAN : spearman’s R; Measure of Association 
Between Rankings of Paired Data. 


coe NPSUR’: Nonparametric Simple Linear Regression; 
Least Squares. 


The list presented above can be displayed at any 
Pete wove, typing: List. 
Germ Sacn cesc program, tunere exists a HOW variable 


that gives a full description of the test and the 


meoerous, OOtLTIONnS that may be exercised. To display any 
See ene HOW variables, Pict onuet wae Fest program ’s 


name with the suffix HOW appended (i.e. SIGNHOW). 

Be ctesteis fun by Entering the orogram’s name. Tne 
user is immediately prompted to input data. Enter 
numerical data separated by spaces or as a variable to 
which the numbers have been previously assigned. 
Several of the tests require a considerable amount of 


prompting before all the necessary data has been 


entered. 
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APPENDIX C 


WORKSPACE FAMILIARIZATION THROUGH PRACTICAL BXAMPLES 


Ig General Information 

This appendix applies to both the Mainirame and 
microcomputer workspaces. Its purpose is to acquaint 
the user with the organization of the programs and the 
type oL prompts to be expected. 

Bxtensive error checking has been included in the 


programs to ensure that the data is of the proper form. 


snould a program become suspended, clear the state 
indicator BbYyeventerans: JRESET, check over the data for 
errors, and restart the program. 33@ kilobytes of 


computer memory are needed to load APL and the NONPAR 


workspace; to avoid filling up the remaining workspace 
area, the ~wWser should mMinimrze daca srorace ewe. 
NONPAR workspace. To exit a program at any time, press 


the Control and Escape keys, simultaneously. 
2% Practical Examples 
a. Sign Test 
(1) Descrintion@et Probl ena A Sinelet.) mine 
is manufactured to nave a median explosive weight of 
not less than 16 ounces. The explosive weights of 15 


mines, randomly selected from the production line, were 
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recorded as follows: Wesel e ome 1>.8 Voeo 16 16.1 
meso oro oe 1G.1 15.7°15.8 15.9°975.8. 

(a) Is the manufacturing process packing 
enough explosives in the mines”? 

(db) What range of values can be expected 


for the median of the explosive weights 90% of the 


comeme 
C25 moot taione To see if the manufacturing 
process is meeting the specifications, we test the 
[perme sis Hoe WM = 16 versus Hi: M < T6. 
(3) Workspace Decision Process. 
(a) Microcomputer: Choose the Sign Test 
prom the main menu, and clic ODtLOn, Single Sample; 


eso Hi: IM —-eMoeversus Hi: M < Mo, from the test menu. 
lp to the Program Interaction section below. 

(bob) Mainframe: Enter SIGN at the keyboard 
and receive the prompt: 

DID YOU ENTER THIS PROGRAM FOR THE SOLE 
PURPOSE OF GENERATING CONFIDENCE INTERVALS FOR A 
SPECIFIED SAMPLE SIZE AND QUANTILE? (Y/N). 

Enter N (If Y is entered, the user will 
go directly to this last option of the test). The next 
Drompct ins: 

THE NULL HYPOTHESIS STATES - THE 
POPULATION MEDIAN (M) IS EQUAL TO THE HYPOTHESIZED 


MEDIAN (Mo); H®: M = Mo. WHICH ALTERNATIVE DO YOU WISH 


2 PS. 


TO) Tee 7 ENTER: 1 FOR Hi: M < Mo; 25FORJHi=) iwc. 
SPOR ee hie ee vier 
Enter 1. The Mesa see emet is: 
ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 2 
FOR PAIRED-SAMPLE PROBLEM. 
Bia cig 
(4) Programminteraction. [hes pregpiuss.- 

ENTER THE DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED ). 

Enter the data separated bY Spaces Or as a 
variable to which the data has been previously 
assignecd. Tie. WexG Digenpt 1s: 

ENTER THE HYPOTHESIZED MEDIAN. 

Enter 16. The following is dislayed. 


COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF: 


lo 

THE TOTAL NUMBER OF POSTIVE SIGNS IS: 3. 

THE P-VALUBE FOR H@: M = 16 Versus Hi: 
Mos 62S: -9461. 

Consider a significance level of .@5. 
Since the P-value of .@461 is less than .05, we reject 


H@: M = 16 in favor of Hi: M < 16 amd conclu@e tia tme ce 
Manufacturing process is not packing enough explosives 
im the, Sanclasim; mine: The next preempt ise 


WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE MEDIAN? (Y/N). 
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Enter Y (If N is entered, the progam asks 
'M—mGOnridence intervals for a quantile are desired). 
mae Next prompt is: 

ENTER THE DESIRED CONFIDENCE COBFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 90. The following is displayed. 

A 90% CONFIDENCE INTERVAL FOR THE MEDIAN OF 
Pees OorUbATION IS: ( 15.8 < MEDIAN < 16 ). 

Them@mext prompt is: 

WOULD YOU LIKE CONFIDENCE INTERVALS FOR A 
SPECIFIED QUANTILE? (Y/N). 

To see the form of the results, we generate 
confidence intervals for the 30th quantile. sample 


size is automaticly set at the number of data points 


entered eariler. Enter Y (If N is entered, the 
mainframe program ends; Or; the Sign test menu 
reappears). The next prompt is: 


ENTER DESIRED QUANTILE; FOR EXAMPLE: ENTER 


BoweEOR THE 20TH QUANTILE. 


Enter 50. The following is displayed. 


ORDER STATISTICS i: COBFFICIENTS 
3 8 FO 2p "oO 
2 9 .94949¢ 
1 1@ ! .99160@ 


evo Gao ke ee GlivViSs ~€ONFIDENCE COBF- 
FICIENTS FOR VARIOUS INTERVALS WITH ORDER STATISTICS AS 


mo POINTS FOR THE S0TH QUANTILE. 
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The mainframe program ends. The menu- 
driven microcomputer program pauses for input from the 
keyboard by prompting: 

PRESS ENTER WHEN READY. 

Press Enter and the Sign test menu 
reappears. | 

b. Wilcoxon Signed-rank Test 
(1) Description of Problem 2. A special train- 
ing program is being considered to replace sine regula: 
training that Radio Telepnone opera .e eee ee. Ce In 
order to evaluate the erfectiveness of the new training 
program, proficiency tests were £1Ven Surine sneer ne 
week of regular training. Twenty-rour trainees were 


chosen at random and grouped into twelve pairs based on 


proficiency test scores. One member of  e€acha spank 
received specialized training while the other member 
received regular Strains Upon graduation, the 


proficiency tests were given again with the following 
resullgs: 

Specially Trained Group (X): 60 505 55. 
43 59 64 49 61 54 47 7@ 


Regularly Trained Group (Y): 40 46 60 53 


re 
WO 


~] 


45 59 40 535 


UI 
UI 
A 


W) 


(a) Does the special trainine program 


ensure higher scores? 
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(b) By what range of values can the 
scores of the two groups be expected to differ 95% of 
ric time ? 

(Ap ecottnton. Towstest the hypothesis that the 
special training program raises profficiency scores, we 
test H@: M(X-Y) = @ versus Hi: M(X-Y) > @. 

(3) Workspace Decision Process 

(a) Microcomputer: Choose the Wilcoxon 
Signed-rank Test from the main menu, and the option, 
Paired-sample; Test H@: M = Mo versus Hi: M > Mo, from 
the test menu. Skip to the Program Interaction section 
ine Ow . 

(b) Mainframe: Enter WILCOX at the key- 
Doard and receive the prompts: 

THE NULL HYPOTHESIS STATES - THE 
POPULATION MEDIAN (M) IS EQUAL TO THE HYPOTHESIZED 
MEDIAN (Mo); H®: M = Mo. WHICH ALTERNATIVE DO YOU WISH 
Oe This Tt? Pitbke: t POR Hl: M < Mo; Zoned > Mos 
> FOR H1: M # Mo. 

Paver 2.) ile next prompt is: 

ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 2 
FOR PAIRED-SAMPLE PROBLEM. 

Batter “2% 

(4) Program Interaction. The prompt is: 

ENTER xX DATA (MORE THAN TWO OBSERVATIONS 


ARE REQUIRED). 


a9 


Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Enter the ¥Y data. The next prompt is: 

ENTER THE HYPOTHESIZED MEDIAN ' FOR THE 
DIFFERENCES OF THE PAIRED DATA. 

Enter @. The following is dislayed. 

COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF: 
ee 

THE TOTAL SUM OF POSITIVE RANKS IS: 60.5. 

THE P-VALUE FOR COMPARING THE MEDIAN OF THE 
POPULATION OF DIFFERENCES TO THE HYPOTHESIZED MEDIAN, 
H®@: M(X-Y) = ® Versus H1: M(X-Y) > @, IS: -@505. 

Consider a significance level o} EE Lo 
Since the P-value of .@505 is greater than .@5, we do 
not reject the null hypothesis that the two training 
COUrCGEes are equally crice ts. However, duce to “the 
closeness in values, the choice of rejecting ene ae. 
rejecting the null hypothesiswis stricta elec cicus: 
call. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 
THE MEDIAN? (Y/N). 

Enter Y (If N is entered, the mainframe 
progam ends; or, the Wilcoxon test menu reappears). 


The next prompt is: 
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ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 95. The following is displayed. 

A 95% CONFIDENCE INTERVAL FOR THE MEDIAN OF 
THE POPULATION OF DIFFERENCES IS: 

( -1 < MEDIAN(X-Y) <« 16.5 ). 

The mainframe program ends. The menu- 
driven microcomputer program pauses for input from the 
keyboard by prompting: 

PRESS ENTER WHEN READY. 

Press Enter andthe Wilcoxon test menu 
reappears. 

c. Mann-Whitney Test for Equality of Medians 
(1) Description of Problem 3. A group of Army 
and Navy officers were given the Defense Language 
Aptitude test. From the results, 14 Army and 17 Navy 
officers’ scores were randomly selected. These scores 
are listed below. 

Army (X): 35 350 55°51 28 25 16 63 60 44 20 
42 47 38. - 

Navy (Y): 54 26 41 43 37 34 39 50 46 49 45 
SrmeINS6 SS-42 34. 

(a) Is there sufficient evidence to claim 
that Navy officers score higher on this test than Army 


officers? 
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Gls} ) slenre what range of values = the 

scores between the two groups be expected to differ 90% 
of the time. 

(2) Seluatene To see if Navy officers score 

higher on the exam, we test H@: Mx = My versus 


Ais ix Gein 


(3) Workspace Decision Process. 


(a) Microcomputer: Choose the Mann- 
Whitney Test from the main menu, and the option, Test 
HO: Mx = My versus H1: Mx < My, from the test menu. 


skip to the Program Interaction section below. 
(bo) Mainframe: Enter MANNWHIT at the 
keyboard and receive the prompts: 
DO YOU WISH TO COMPARE THE MEDIANS'- OR 
VARIANCES OF THE POPULATIONS? ENTER: 1 TO COMPARE 


MEDIANS; 2 TO COMPARE VARIANCES. 


Enter 1. The next prompt as: 

THE NULL HYPOTHESIS STATES - THE 
MEDIANS OF XK AND Y ARE EQUAL; Mx = My. WHICH 
ALTERNATIVE DO YOU WISH TO TEST? ENTER: 


Te PORS HT 2 Mec hye. 2 FOR Mx > My; % FOR Mx # My. 
Enters). 
(4) Program Interaction. The prompt is: 


ENTER xX DATA (MORE THAN ONE OBSERVATION IS 
REQUIRED). 
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Enter the X data separated by spaces. The 
next prompt is: 

ENTER Y DATA. 

Enter the Y data. The following is 
displayed. 

THE SUM OF THE X RANKS IS: 224. rie Uy 
STATISTIC EBQUALS: 119. 

THE P-VALUK FOR H®@: Mx =My versus H1: 
Mex < My IS: . 50 TS. 

We donot reject the hypothesis of equal 
population medians and conclude the median of all Army 
scores is equal to the Navy’s. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 


THE SHIFT IN LOCATION (My - Mx)? (Y/N). 


Enter Y (If N is entered, the mainframe 
program ends; Or the Mann-Whitney Test menu 
reappears). The next prompt is: 


ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 95. The following is displayed. 

A 95% CONFIDENCE INTERVAL FOR THE SHIFT IN 
LOCATION BETWEEN POPULATIONS X AND Y IS: 

( -10 < My-Mx < 10 ). 

The mainframe program ends. The microcom- 


puter program pauses for input from the keyboard by 


PrOMpi ing : 
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PRESS ENTER WHEN READY. 
Press Enter and the Mann-Whitney test menu 
reappears. 


d. Mann-Whitney Test for Equality of Variances 


(1) Description of Problem.. Referring to 
problem 3 in. SeGtvonumecial )- Is there sufficient 


evidence to claim that Army scores vary more than Navy 


scores? 
(G2) Sent seine To see if Army scores vary 
more, we test H@: Vx = Vy versus H1i: Vx > Vy. 
(3) Workspace Decision Process. 
(a) Microcomputer: Choose the Mann- 
Whitney Test from the main menu, and the option, Test 


H@: Vx = Vy versus Hi: Vx > Vy, from the test menu and 
receive the premnoc: 

ENTER THE DIFFERENCE OF THE MEANS ' OR 
MEDIANS (Mx - My). 

Because we believe the population 
medians to be approximately equal, We enter @. Skip to 
the Program Interaction section below. 

(db) Mainframe: Enter MANNWHIT at the 
keyboard and receive the prompts: 

DO YOU WISH TO COMPARE THE MEDIANS’ OR 
VARIANCES OF THE POPULATIONS? ENTER: 1 =TO COMPARE 
MEDIANS; 2 TO COMPARE VARIANCES. 


Enter 2. The next prompt is: 
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THE TEST TO COMPARE VARIANCES, REQUIRES 
THE TWO POPULATION MEANS OR MEDIANS TO BE EQUAL. IF 
THEY DIFFER BY A KNOWN AMOUNT, THE DATA CAN BE ADJUSTED 
BErORm APPLYING THE TEST. ENTER THE DIFFERENCE OF 


MEDIANS (Mx - My) OR 900 TO QUIT. 


We enter @. The next prompt is: 
THE Nome SHY rPOTHESTIS STATES - THE 
VARIANCES OF XK AND Y ARE EQUAL; Vx = Vy. WHICH 


AOPERAATIVE DO YOU WISH TO TEST? ENTER: 
1 FOR Hi: Vx < Vy; 2 FOR Vx > Vy; 3 FOR Vx # Vy. 
Enter 2. 
= eatogenelinteract@en. The prompt is: 

ENTER xX DATA (MORE THAN ONE OBSERVATION IS 
REQUIRED). 

Enter the X data separated by spaces. The 
Wet Drompt ws: 

ENTER Y DATA. 

Enter the Y data. The following is. 
displayed. 

TSoeeoUMe OF ihe X RANKS £5: 166. THE U 
STATISTIC EQUALS: 61. 


Toe P=VAGWwE FOR H@: “Vx = Vy versus H1: 


vx > Vyers: mien 2. 

Consider a significance level of Oo 
Since a P-value of .9112 is less than .@5, we reject 
the null hypothesis of equal variances in favor of 
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Vx > Vy and conclude that Army scores do vary more than 
Navy scores. 

The mainframe program ends. The 
microcomputer program pauses for input from the key- 
board bY prempyins: 

PRESS ENTER WHEN READY. 

Press Enter and the Mann-Whitney test menu 
reappears. 

e. Kruskal-Wallis Test 
(1) Deseription of Prebpiem 4. During@= a receas 
Monster Mash involving four Navy SEAL Teams, one of the 
events consisted of the number of pushups aman could 
do in 2minutes. Hight men were chosen randomly from 
each Team. The following scores were recorded. 

SEAL 1: 90 96 102 85 65 77 88 7@. 

SHAL 2: 64 79 9959S 5e75i oo oie 

SRAL 3: 101 66 93 89 7 isa (Guo 

SEAL 4: 72 78 73 81 85 92 94 86. 

Are the different Seal Teams cConsiderea sce 
be equally fit? 

(2) Solution. To see if the Seal Teams are 
equally fit. we test the hypothesis that all the 


population medians are equal. 
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(3) Workspace Decision Process. 

(a) Microcomputer: Choose the Kruskal- 
wer is test from the main menu; and, once the test menu 
is displayed, press Enter. 

(bob) Mainframe: Enter KRUSKAL at the 
keyboard. 

(a) Program Interaction. “Ihe prompt is: 

ENTER THE NUMBER OF POPULATIONS TO BE 
COMPARED (MUST BE GREATER THAN TWO). 

Enter 4. The next prompt is: 

ENTER YOUR FIRST SAMPLE. 

Enter the SEAL 1 data separated by spaces. 
Pe next prompt is: 

ENTER YOUR NEXT SAMPLE. 

Enter the SEAL 2 data. The next prompt is: 

ENTER YOUR NEXT SAMPLE. 

Enter the SEAL 3 data. The next prompt is: 

ENTER YOUR LAST SAMPLE. 

Enter the SEAL 4 data. The following is 
displayed. 

THE H STATISTIC EQUALS: mio D 

THE P-VALUE FOR H@: THE POPULATION MEDIANS 
ARE EQUAL versus H1: AT LEAST TWO POPULATION MEDIANS 


ARE NOT HQUAL IS: See). 
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We do not reject the null hypothesis that 
the population medians are equal and conclude that the 
SEAL Teams are equally fit. 

The mainframe program ends. The microcom- 
puter program pauses for input from the keyboard by 
DIS OUND) canoes 

PRESS ENTER WHEN READY. 

Press Enter andthe Kruskal-Wallis test 
ae reappears. 

tf. Kendaiiis se: 

(1) Description of Mrroblemms- In order to 
determine if cold weather affects target marksmanship, 
Naval Special Warfare recorded small arms marksmanship 
scores and corresponding air temperatures for a period 
of one year. 2®@ men were chosen at random, and their 
scores averaged for different air temperatures. The 
average score for each air temperature is shown below. 

Air temperature (X): 5®@ 55 20 50 65 55 30 
52 40 6@. 

Average scores Cy) 212 20@ 165 165 260 
215 17D WOT seCe2oe. 

Can it be said that colder temperatures 
have an effect on marksmanship scores? Is that effect 


positive or negative? 
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(2) “so imGtion. We test the null hypothesis 
that no association exists between cold temperatures 


and marksmanship. 


(3) Workspace Decision Process. 
(a) Microcomputer: Choose Kendall’s B 
Test from the main menu; and, once the test menu is 
displayed, press Enter. 
(db) Mainframe: Enter KENDALL at the 
keyboard. 
(“Ge reoccatmimeteraccton = The prompt is: 
ENTER X DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 
Enter the X data separated by spaces. The 
next prompt is: 
ENTER Y DATA (NUMBER OF Y ENTRIES MUST 


EQUAL NUMBER OF X ENTRIES). 


Enter the Y data. The following is 
displayed. 

KENDALL’S B EQUALS: 1 Set 

THE P-VALUE FOR H®@: NO ASSOCIATION EXISTS 
versus: H1: DIRECT ASSOCIATION EXISTS IS: -@0045. 


THE P-VALUE FOR THE TWO-SIDED TEST OF 
PYPOrHmols IS: .@009. 

Since the P-value for the one-sided test 
equals .0@0045, we reject the null hypothesis that no 


association exists between temperatures and 
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marksmanship in favor of direct association. We 
conclude that colder temperatures tend to cause lower 
marksmanship scores. 

The mainframe program ends. The microcom- 
puter program pauses for input from the keyboard by 
prompting: 

PRESS ENTER WHEN READY. 

Press Enter and Kendall’s B test menu 
reappears. 


gg. Spearman’s R 


(1) Description of Problem 6. When fitness 
reports are written, officers of the same grade are 


ranked against each other based upon their demonstrated 
level of performance. Last marking period, the 
Commanding and Executive Officers separately ranked 9 


Ensigns as shown below. 


Ensigns 
A BCDEFGH I 
CO (%)e 6 40465 2 8 3 ao 
KO (CY): 5°96 3 2 eo mec 


Does any association exist between the two 


sets of rankings? 


(29a So lumaone We test the null hypothesis 


that no association exists. 
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Po me vornkaeace Decision Process. 
(a) Microcomputer: Choose Spearman’s R 
Test from the main menu; and, once the test menu is 
agreplayedws press Enter. 
(db) Mainframe: Enter SPEARMAN at the 


keyboard. 


GP) Frogram Interaction. The prompt is: 
ENTER xX DATA (MORE THAN TWO OBSERVATIONS 
ARE REQUIRED). 
Enter the XK data separated by spaces. The 
next prompt is: 
ENTER Y DATA (NUMBER OF Y ENTRIES MUST 


EQUAL NUMBER OF X ENTRIES). 


Enter the v7 data. The following is 
displayed. 

SPEARMAN’S R EQUALS: - 200. 

THE P-VALUE FOR H@: NO ASSOCIATION EXISTS 
versus: H1: DIRECT ASSOCIATION EXISTS IS: -@664. 

THE P-VALUE FOR’ THE ‘TWO-SIDED TEST OF 
Peer OlHBoIS IS - .13528. 

Consider a significance level of  .@5. 
Since a P-value of .9664 exceeds .95, we do not reject 


the null hypothesis that no correspondence exists 


between the two sets of rankings. 


el 


The mainframe program ends. The microcom- 
puter program pauses for input frem thes keyboard) jay 
prompting: 

PRESS ENTER WHEN READY. 

Press Hunter and the Spearmeanas@.8ccsG. mon 
reappears. 


h. Nonparametric Simple Linear Regression; Least 
Squares 


(1) Description of senoplem 7. Battery-powered 
Swimmer Proplusion Units are sometimes used to aide 
swimmers during long underwater swims. Recent tests 
have shown that a nearly linear relationship exists 
between water Len pe cere and? Dattery Lite ter egese 
Und tS. The following 17 data points were randomly 
selectedetroemache test ere mnt. 

Water temperature (X) Battery life (Y) 
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(a) Find the fitted regression equation. 


eZ 


(bob) For the following water temperatures, 
bwemetet thew bautery life of the units: 61 52 46 356. 

(c) Can we determine with any certainty 
if the slope of the regression line equals .@5. 

(d) What range of values could be used 


as the slope of the estimated equation line 90% of the 


time? 

(2 )*Setution- To determine the estimated 
peeression equation, we use nonparametric linear 
regression. 

(3) Workspace decision process. 

(a) Microcomputer: Choose Nonparametric 
Simple Linear Regression from the main menu; and, once 


the test menu is displayed, press Enter. 

(bd) Mainframe: Enter NPSLR at the 
keyboard. 

C= Paeeran Interaction. The prompt is: 

ENTER xX DATA (MORE THAN TWO OBSERVATIONS 

ARE REQUIRED). 
| Enter the X data separated by spaces. The 

nest Drompt is: 

ENTER Y DATA (NUMBER OF Y ENTRIES MUST 
EQUAL NUMBER OF X ENTRIES). 

Pai er tae Y data. Pico iol Lowing 1s 


displayed. 


73 


THE LEAST SQUARES ESTIMATED REGRESSION 
BQUATION IS: 

Y = =-1. 26550 2 Vevcco 

The next prompt is: 

DO YOU WISH TO ENTER SOME X VALUES TO GET 
THE PREDICTED Y’S? (Y/N). 

Enter Y¥ (If N is entered, the program skips 
to hypothesis testing for the slope). 

ENTER X VALUES. 

Enter 61 52 46 36. The next prompimae- 


THE PREDICTED Y VALUES ARE: 2.44 Tee 


WOULD YOU LIKE TO RUN SOME MORE X VALUES? 

Enter N. The next prompt is: 

WOULD YOU LIKE TO THST HYPotTeesis Clee e las 
SLOPE OF THE EQUATION? (Y/N). 

Enter Y¥Y (If N is entered, the program skips 
to confidence interval estimation). The next prompt 
is: 

ENTER THE HYPOTHESIZED SLOPE. 


Enter .@5. The following is displayed. 


SPEARMAN’S R EQUALS: a5 Looe 
THE P=-VALUE FOR HO: 385-9005 Versueea |) eee 
05 IS OU a 


THE P-VALUE FOR THE TWO-SIDED TEST OF 


HwPOTHESks 15; 201 3S. 
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Consider a significance level of .@05. Since 
a P-value of.@9079 is less than .95, we reject the null 
Mmypothesis “that B “= 2.05 in favor of B > -@5, and 
conclude the slope of the regression line is greater 
than .@5. The next prompt is: 

WOULD YOU LIKE A CONFIDENCE INTERVAL FOR 


THE SLOPE? (Y/N). 


Enter Y (If N is entered, the mainframe 
program ends; or, the Nonparametric Regression test 
menu reappears). The next prompt is: 


ENTER THE DESIRED CONFIDENCE COEFFICIENT; 
FOR EXAMPLE: ENTER 95, FOR A 95% CONFIDENCE INTERVAL. 

Enter 90. The following is displayed. 

A 90% CONFIDENCE INTERVAL FOR B, THE SLOPE 
OF THE ESTIMATED REGRESSION LINE, IS: 

(Seis o> <5 < 2.07 ). 

If the estimated slope, does not lie within 
thie confidence interval, the following would be 
displayed. 

THE LEAST SQUARES ESTIMATOR OF B_ LIES 
OUTSIDE THE CONFIDENCE INTERVAL. DISCARD THE LEAST 
SQUARES EQUATION AND USE: 

Y= -1.4458 + .960833xX. 

THIS EQUATION IS BASED ON THE MEDIANS’ OF 
THE xX AND Y DATA, AND THE MEDIAN OF THE ‘TWO-POINT 


SLOPES CALCULATED FOR THE CONFIDENCE INTERVAL ON B. 


Ve 


The next prompt is: 
DO YOU WISH TO ENTER SOME X VALUES TO GET 


THE PREDICTED Y’S FROM THE NEW EQUATION? (Y/N). 


To compare results, let us input the 
temperatures in the new equation. Enter Y (If N is 
entered, the mainframe program ends; On, the 
Nonparametric Regression test menu reappears). The 


next prompt is: 

ENTER X VALUES. 

Enter 61 52 46 36. The following is 
displayed. 


THE PREDICTED Y VALUES ARE: 2.265 1. 72 


The next prompt is: 

WOULD YOU LIKE TO RUN SOME MORE X VALUES? 
Enter N. The next prompt is: 

WOULD YOU LIKE TOSEERS Tey PotnkclewON) Bie mes 


SLOPE OF THE EQUATION? (Y/N). 


To compare results once again, we enter Y 
(If N is entered, the mainframe program ends; or, the 
Nonparametric Regression test menu reappears). The 


next prompt alee 
ENTER THE HYPOTHESIZED SLOPE. 
Enter .@5. The following is displayed. 


SPEARMAN’S R EQUALS: So ore 
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Mabe P-VALUE POR H@: B= .05 versus Hi: B > 
we mel: -O@00Q. 

THE P-VALUE FOR THE TWO-SIDED TEST OF 
HYPOTHESIS IS: .@@@0. 

The mainframe program ends. om ibiaet. oC OMm— 
puter program pauses for input from the keyboard by 
Proemp cine: 

PRESS ENTER WHEN READY. 

Press Enter and the Nonparametric 


Meeression test menu reappears. 


of 
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APPENDIX D 


MAIN PROGRAM LISTINGS FOR MICROCOMPUTER WORKSPACE 


V Be AsAA;B3BX3873C:CX GID s DD; Px: DY; DXY;3;S3POS;NEG;XX;JY;N3DEN3NN3NUM;P; 


Vals 35g sSViT Us VAs Z:X3¥ :WW3CHA:E:PV3Q:R° 
Q Carehr by’ COMPUTES THE LeNDALE B SPATISTIC WHICH IS A MEASURE 
a OF rae RTA FION BETWEEN SAMPLES. P-V RE GIVEN FO TESTING ONE 
a AND TWO-SIDED HYPOTHESIS FOR NO ASSOCIATION VERSUS ASSOCIATION. 
a SUBPROGRAMS CALLED 3Y THIS FUNCTION INCLUDE: TIES, TIESK, KENDALP, 
a INTERP, INPUT AND NORMCDF. 


A DISPLAY TEST MENU AND INPUT DATA 
Ni: ee KENQ3J 


MENO MAINOBJ 
Bi: ee 2 


+1 
+14 (Q+1)42 
Be 
ORDER Y IN INCREASING ORDER OF X 
"A+ZCAX] 


Q ORDER X IN INCREASING ORDER 
BX CAK] 


A Ae COMPUTE CURRENT RANKING OF I 
< 
BE oat ane NOW ORDER Y RANKS IN INCREASING ORDER 
<< 

I® CIES EXIST IN EITHER X OR Y RANKED VECTOR USE MID-RANK METHOD 
DD+1 TIES D 
tel PIES 3 

FIND ORIGINAL RANKING OF I WITH TIES RESOLVED 

"YY+DDCC} 


A COMPUTE fone OF DISTINGUISHABLE PAIRS aN 
Sepo 1))+ 
+90 


AA+ 
H POSITIVE ONES COME FROM A RUNS UP CONDITION; NEGATIVE 1 FROM RUNS DOWN 
a ZERO IS SCORED FOR TIES. MULTIPLY THE RESULTS FOR EVERY ELEMENT AND SUM 


“BRS CERCA)» (aayxxy) 


D 


CX+ &x AA} <(AAtYXX 
<_ 

B+ (ZZCAa > (AAVTY)) _ 

CY¥+(YY(AA}< (AAV) )x(71) 

DY+BY+CY 

DXY*«DXxD 

POS (DX¥>0) _ 

NEG+«(DXY<0 )x(71) 

S+S,POS,NE 

a? (AAS (= -1))/51 

ae SUM FINAL VECTOR TO DETERMINE S 

i ree IN THE NUMBER OF TIES IN EACH VECTOR USING THE TIESK FUNCTION 
< 

V+TIESK D 

SU*+/ (210) 

oiV 


SV++/ 
A aes THE B STATISTIC INCLUDING THE CORRECTION FOR TIES 
+ ((NN-SU )x(NN-SV))*0.5 


are a 
+(N>13)/NORM 
* pe REND ALP CALL KENDALP TO CALCULATE THE RIGHT TAIL OF TRE CDF OF B 


a abb NTERP Ae ao nbCuULars P-VALUE BY INTERPOLATION 
PVA Leal INTE cRP 

>(PVA4L= 1)/2 

PVAL+«0.5 


+53 
a CALCULATE P VALUE USING NORMA 
NORM: SHON RAs eC Fat 8 85 tt, Aes 
DEN+(2x ((2xN)+5 ))«0 
Z+NUMtDEN 
Pe a sO Re ORE Ve pprnt OUT DIRECT 
A E 
Peedreu ee Cl ASSOCIATION. 
Chas! INDIRECT 
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A+«'DIRECT' 
<2xPVAL 
$1)/L8 


-~coty + ot + 
=F gtuA 


cee 
try = 
Oyrshe og 


8) 
L KAEROALL''S 8 Bane (4s7) ,OTCNL 
P-VALUE FOR HO: 6 asso OCrATION EXISTS VERSUS' 
t Bl at CEADT" ASSOCIATION EXis7S JS: ear O7CNL 
'THE P P-VALUE FOR rae rio. *s1DED TEsi Of BIPOTHES!IS JIS: 1 tusPy), TCNL 
'DPESS ENTER WHEN AEADY. 
WW<9 
3 Nt 
v 


CONT ee reeawenre a ae e 
WO 00 00 60 CO CO NINN 
NEOMPOROWO IMAM FE 


V AAWL :NUM;DENOM+AiC +9; 5:K:AA:BB;DDiEsP :PiPVAL: Ei SOPR:SR:PSOR;CHA:B 
TH2S FUNCTION COMPUTES Taz XUSKAL Siipiis risr sharrorie a wHICd 2S 


250C‘(a 

3% 9 3 MEASURE OF "HE FQUALITY OF < ZNDEPSNDENT SAMPLES. 

33 a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: TIzS INDEX2LS 
41 a FOTSTN, INTERP AND THE VARIA BLES PMATERW 20, bwareeaa SL Te Dea hxyau 
3 a DMATEWU1, PMATKWUu2, AND PMATXWU3 

71g MENU CHOICES AND ROUTE TO PROPER STATEMET FOR ACTION. 

3} Ni: B+MaNU KRWLOBT 

g +(B=1)/81 

10 MENU MAINOBJ 

ie) ei, 

as 


335) 2h ZNTER THE NUMBER OF POPULATIONS TO BE COMPARED (MUST 3E CREATER THAN T 


~ 


= ($R<3 Jv ((11K)20))/24 


SII eee oer 


pK)> 
R INITIALIZE VECTORS E AND F AND VARIABLE C 
ae ore? 
aeeaet Preset PACILITATES ENTERING THE SAMPLE VECTORS AND STORING THEM 
Caa+ os 


« — -T 


Bie 
RENTER YOUR ',(3CHA),' SAMPLE.! 
=(¢ogD )= 0)/NEXT 
R CONCATENATE SAMPLES AS THEY ARE ENTERED AND STORE THEM IN VECTOR E 


NEXT 3: E+E , D 
R "RECORD THE LENGTHS OF THE SAMPLES AS THEY ARE ENTERED 


NWNNNVNho hth +4 bp 


SAM EWNRPROWODDINAM FEWNHPOWDINMDM FCLWNHPrPOWODIMDMNE WHitonuw nu F&F 


en 
>(C<K)/L1 
a RECORD SIZE OF ALL SAMPLES WHEN COMBINED 


NetZe 
AReerrta SAMPLE SIZES LARCEST TO SMALLEST 


oF LWA) 
Rg CAF] ORDER COMBINED SAMPLE VECTOR TO BE USED BY TIES FUNCTION 
+ 
R CALL INDEXPLS TO INCREMENT INDEXES WHEN TIES OCCUR WITHIN ONE SAMPLE 
AAtF INDEXPLS E 
R CALL TIES TO BREAK TIES BY MIDRANK METHOD 
aa Llés D 
“ 
THIS LOOP CALCULATES THE #2 STATISTIC 


A 
ce MOF RANKS FOR EACH SAMPLE IS U 
IS CALCULATED 
Sha+/BBt (tel AACE ))]J e 
CALCULATE SUM P BANKS SQUARED DIVIDED BY THE INDIVIDUAL SAMPLE SIZE 
SR+(S2*x2)+F(C] 


SOERSSOFR SR 
+(C<K)/L2 
2 SUM ACROSS ALL SAMPLES 

ae ee GATGULATE FNAL 2 § Sir 

: r Rares ene 
"ae (zsopn (sae Cantgetys 2) tae Cees ge rf 
a ZGa C C IES 
A+TIESK E 
NUM+(+/(Ax*3))- (+/A) 
DENOM+Nx ( (Nx2 )-1) 
H+eH+(1- (NUM DENOM) ) 
M GF PD HeSEGAL STATEMENTS ENSURE PROPER PROB. IS ACCESSED 
OX 


dD 


Db 


SroRE EACH CALCULATION 


vd 


oo oe ee oe oe ee eo oe oe ee oe ee ee ee ee ee ee ee ee ee er ee ee ee ee ee ee ee 
AM 
Dy oy + 
5 Cy os 

A te 
rN = 
briny 

u 

pa 
Hw 

VS 

~ 

t 

js 


R 


gw 


STE 
J 
/ 
0 OF(1]>3))/FAPPROX 
F Ceor 


ENON NTNETNTETETNEEEAENEET en 
NOAM AOMOHOOOOWWGHIANMNNUM EEE FEEL EF FWWWwWWwwwwWWwhryb . 


OwWONIMM FWhyOrHOwo 


Fever? 


a2 


C71) Rong ee 3°21 1 1))v(A/(OF= 3 2 1 1))VCA/(OF= 2 3 1 1) v4) Cr eee 








>(K=4 )/Pu3 
}>4)/FAPPROX 
= 21 1))vCA/(COFR= 354, 12)0)0 7207 ecu7 
=u )ACOF( 31 s=1)) 7225 
= 2 A 1 3 ene 3.3 .1))vQa/COF= 3 3 2))VCa/COF= 3 353700723 
= 2))v(OF(1J=4))/P34 
CALL ete aes VARIABLE ACCESS CDF 


3) 
P3i: P<{PMATKW31C0 (N- a 
P33: P+PHATEW33((N-5)3 
P34: P<PMATKW34L(N-6)33] 
Pui:P+PMATRWHYLC(N-5):; 
PU2; PaPMATRW42( (N-5)3 5] 
Pu3:P<«PMATRWU3((N-7):3] 
A CALL INTERP LO CALCULATE P-VALUE BY INTERPOLATION 
PM: PVAL~H INT ERP 
eres 13) OUfPUT 
a CALCULATE P-VALUE USING THE F DIST W/ONE LESS D.F.IN DENOM APPROX. 
ae eae )xH )+(CCK- ae See 
bes Le1-(((K-1),¢C(N-K)-1)) FDISIN F 
On7 PUT: PVAL<'GREATER THAN 223! 
La:'THE BH SUAPISTIC EQUA. (usd) .OTCNL 


THE P-VALUE FOR HO: THE soporarton MEDIANS ARE EZQUAL aes s7s ' sd2CNE 
gi: Af ZEAST THO POPULATION MEDIANS ARE NOT cQUAL 25: |, (4USPYAL), 


J 
3 
2 
j 
orcyr 
‘DRESS ENTER WHEN READY. ' 
1 Aa+G 


ONT TIT 
FPHRPPP PERO WDWWW DOD NDDDD MDM MOWONAII IASI) 
OO00O9O BOO OWDANIOMNEONHONMOODNOHMNFEWNPOWONMAMEWhN 


sailet 
E1: "ERROR: ENTER A SINGLE INTEGER VALUE GREATER THAN 2; TRY AGAIN.' ,OTCNL 
Nit . 


V MANW:N:M:PV2: G:MM:NN:RX;0:NM1;P:NU;PVAL;NM:NUMZ;NUMZ1;DEN;DENC:DE 
NC3 5 ie rea snd 2 isdePaa; -CDF: s INDEX ; sFahe CI: UALPHA;8B;CC:U1;02;PV;NN1:NN 
bvr: bir gni23 roe GG:A 1a PV PV 93 
2 eT NehtON eee wai S P PRANKS: ROCEDURE TO CALCULATE THE MANN- 


A 

A WATPNEY aks sTIC WHICH 3 SEE IN SA THE P-VALUE FOneareE Tes” 
Rn OF LOCATION aND SCALE. Pu C.t. FOR (MY-MX), THE SHIFT IN LOCATION IS 
a ALSO COMPUTED. SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: ILLES; 1 leo 
A ras DEXPLS, VARMW, MANWP, INPUT, CONFMW, NORMCDF, AND NORMPTE. 


N1: TENU MANWHELP 
A MENU U CHOICES AND ROUTE TO PROPER STATEMENTS FOR ACTION. 
N3:D+CHOICEM PAGEDMENU MANWOBJ 
DIFF« <0 
SAREE a. 934et 
>+(D= 
+(D= tS /N4 
MENU MAINOBS 
a3 pENTER THE DIFFERENCE OF THE MEANS OR MEDIANS (MX - MY).! 
>((pDIFF)>1)/E3 
A ENTER DATA VECTOR 
See X DATA (MORE THAN ONE CL OE PVATION IS REQUIRED). 


= Pens =0)/E1 
ER Y DATA.' 


Hed 
"ve See INVOLVE VARIANCES ADJUST X BY THE DIFFERENCE IN MEANS 
CONCATENATE X AND Y SAMPLE VECTORS 


nde M.A AA 


Ae, M 
ee e DETERMINE SIZE OF & AND Y VECTORS’ AND ASSIGN TO NN AND MM 
NN«+o 
ar 
COMPUTE SIZE LIMIT OF LEFT TAIL OF NUL DISTRIBUTION 
WM (NNXMM 542 
eet 


ORDER A AND ASSIGN TO B 


Et ie) ee A 
8 et wreee : TIES FUNCTION TO BREAK TIES USING MIDRANK METHOD 


A we. FALSE CALCULATE TEST FOR VARIANCES 
Pbiics VARMW TO GENERATE RANKS REQUIRED FOR VARIANCE TEST 


Cc 
" GG+VARMW (NN+M 
a CALL TIES TO "RECORD TIES IN THE DATA AND BREAK TIES IN GG 


oe ce ee ee ee ee ee oe i | 
FEE FE EWWWWWWWWUWINNNNNNNNNPPHPPPHPPHPHPBPOONOUEWNH 


EWN POWDNOMWNE WN OWONHMHMNEWNHPOWONOMNE WNP Olssz wa 


— aS ae eeeaee ae e e e 
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LL TT ae 


6462-R RPRRPP PRP PPE .PHRPPPYD DODO DODO DOW ODDO WOO WWOWWOODY)~JYQOYQWI YQ QYURNOIMNAHOMMOONOOHMMnMNIONNIMNIM FE EE Fe 
AAA R RPRPRPRP POO OOO DOKOC OOO AIHW EWNVOPOW DOWN FWNYHROW ONO FWP POW OID) FWHOPOW OND EFWNRP O100O NOM 


POWWO NOME WNYIPOW WO IMI £Upory Owo0 SO F-WNH Oli i w an nt ll a 


a a a a ee ee eee ee ee ee eee Oe ee ee ee ee ee ee ee ee ee oe ee oo oe oe ow oe 


PRRPRPPRBPHHP PPP) 
WWNNMNNNNNNNNE: + 


‘ 
era was ea eee a ee TT OO ON iO EE EE aaa ene ool 


Tre 

Bs:ax eer che nO SRAM _ 
meet! Ne CONVERT TO MANNWEIT U STATISTIC 
Be ( (NNx (NN4+1 ) )#2) 


IP SIZE OF X TIMES SIZE OF Y > 803; GO TO NORMAL APPROX 
7 5 (WMx2)>80)/L2 
NNi+« /WN MM 


aiaie FUNCTION CALCULATES LEP? TAIL CUMULATIVE PROBS. OF U STATISTIC 
P+NN2 MANWP NN1 


)/Z10 
LOGicaL STATEMENT ENSURES ONLY LEFT SIDE OF NULL DIST IS USED 


A 
+(USNM1)/L3 

Q CONVERT U STAT WHEN GREATER THAN LEFT TAIL VALUES 
U+(NNxMM)-U 

a ali ae FRACTIONAL, INTERPOLATE P VALUE 

ae gt 0)/NON 


220350) /?4 
PV<2-((PCU2+1) D4 2) 


oa S oyet-((PCU21+PC02+11)42) 


PS: PU I< (PCU2+1) +PlU2+2 J+ +2 
NON: +(U>0)/CO 
PV<1 


PZ 
gO: DTet Pry 
Eyre. ot (741) 


33221 [CoiSwitt)/zu 

Byop oy 

evI+PV2 

a "TP CONFIDENCE INTERVAL DESIRED GO TO L10 

a COMPUTE 78) yore AL APPROX! “MATION W/CORRECTION FACTOR 
3)-NM 

Senet COMMAND GEMCNN= 1))#12)*0.5 

Z<NUMZ+DEN 

21+NUMZ1+DEN 


NM ) 
DENCL©< (CCWN4Mi=4 x (DEN#2 334 (NSHM=23)=< 


EDESTN 26) s2 


a 
+MM~2) TDISTN TC1)))+2 


D 
(NN+MM-2) TDISTN BO 
2) TUTSIN TE) se 


M)/SECO 
Bese the CDF _Z)+((N 
es 1-(NORMCDF 21)) 


N+MM-2 
+(1-(¢ 

SECOND: PVI+( (NORMCDF 2d+f 
) 


-2) 
((NN 
1-(( 
PV+((1-(NORMCDF 21) NN+MM- 


LY: aye Eva 
ee N5 


3«4 
N5:PVM+«(3,1)p(PVI,PV,PV3) 

AEBE SUM OF THE X RANKS IS: ',(s8RX),'. THE U STATISTIC EQUALS: ',(8U1), 
A 


oe 
1698 Bale MX = MY VERSUS Bier 5 (SLOGIC(D-1;1))," MY IS: ', 
p 
a 


(PV, PV, PV3) 
ron = VY VERSUS 81: VX ',(sLOcIc{(D-53:1]),! VY Is: ', 


1 787 
HEN’ READY. 


ct. 


PE 
(6 & 3PVM 


>N3 
L8:'WOULD Lou EEN A CONFIDENCE INTERVAL FOR THE SHIFT IN LOCATION (MY - 
MX)? (Y/N).',O7CNL 


nw ( 332! 0") /N3 
Beoscceiyrur 35 
wbPdas (290-2 320 00 
WOOL 2 20 ReetAe aor OAs Er OR@GONE. INT. OF GARCLR SAMPLE SIZES 
"> ((WMx2)>80 9/L5 


"CDF P COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 
he 
INDEX POSITION OF VALUE IN CDF s ALPHA 

“INDEX* (47 (CDPSALPHA)) 

nue ae 0)/L6 

ae 

> 
a COMPUTING Re aoa Re USING NORMAL APPROX. W/C.F. 
L353 Poe ee ee 

UALPHA+ DENS (NORMPTH * ALPHA) ao 
A ROUND UALPHA DOWN “NDS INCREMENT BY ONE 


81 


INDEX+| UALPHA+1 
L6 :iPX«NN , INDEX 


CI+IPX CONE MH A 
A ti(2 CONFIDENCE INTERVAL FOR THE SHIFT IN ' 
' LOCAPTON ‘BEDWERY | POPULATIONS X AND Is: ,07C NL 
oe et ae Ss Mi = MX s-",($€7021),° >  (@tene 
‘BRESS ENTER WHEN C orAbe 


& 


+N3 
Sepa SAMPLE CONTAINS LESS THAN TWO ENTRIES; TRY ACAIN.' ,OTCNL 
= 
Eee ee YOU dAVE ENTERED MORE THAN ONE VALUE. TRY ACAIN.' ,OTCNL 


+ 


MEAP Ppp papapapage 
FFF FE EWWW wwWWwu 
FWNPOWODYA Mu kwr 
LJL_KH_ ALAA M8 


V NPLE;N;SUMX:SUMY :XBAR: YBAR; SUNX2 s SUMXY iB: A:WW3XX;BB;U;D;ALPHA;P3CC;CDF; 
LeeA LOPES ;22:52:77:DENOM: INDEX: (are ear Coen DY 

A S5ERIM CONDUCTS NONBIGAMETSIC LENELS REGRESSION Tas gest SQUARES 
Q ESTIMATED REGRESSION LINE IS COMPUTED WITH HYPOTHESIS TESTING AND 
q CONFIDENCE INTERVAL AVAILABLE FOR THE SLOPE B. I? B DOES NOT rie IN THE 
9 aN ALTERNATS 2EGRESSION SINE <5 PROPOSED. SUBPROGRAMS CALLED ARE 
A Senin, RKENDAL?, NORMPTH, tear’ AND CONFLa. 
A 


P< 
DISPLAY MENU AND INPUT DATA. 
N1:E*«MENU NPLAOBI 
+(EB=1 /N2 
new MAINOSZ 
N2: pe ae 2 
Q+1+2F 
Weeenc os ne 


INMP PPP BME WONGUE Wns 


TT Ea 


poe 
Aa ASSIGN THE SIZE OF X (AND XY) ZO N 
6 

9 COMPUTE THE SUM OF X'S AND Y'S 
SOMX«4+/X 
SUMY<+/Y : 

f 3 COMPUTE [HE MEAN OF X AND YZ 

7 XBAR«SUMX+N 
YBAR+SUMY+#N 
Q COMPUTE THE SUM OF THE X'S SQUARED 
SUMX2<+/(X*2) 
Q COMPUTE THE SUM OF X TIMES Y 
SUMXY<+/ (Xx 
a COMPUTE THE SLOPE OF THE ESTIMATED D LEAST SQUARES REGRESSION LINE 
pe ¢ (NXSONXY) - UsUNExSUMY) )4 (SUNS ) = (SUM 
A COMPUTE , HE Y-INLTERCEPT 


See eo 
ee Pete Rises Fe ee eC RGuTEL ae EQUATION IS:',O7CNL 
'DO rou prise ZO "ENDER "Rone x VALUES TO GET THE PREDICTED Y''S? (Y/N). 


WW+) 

>+(WW='N' )/D4 
L2:'ENTER X VALUES.' 

XX<O 


A CALCULATE PREDICTED Y'S 
YY<A+BxXx 


'TRE PREDICTED Y VALUES ARE: (37z) OPCNL 
WOULD YOU LIKE TO RUN SOME HORE VALUES? (Y/N). 


seals "72 
bis 'HOULD YOU LIKE TO TEST HYPOTHESIS ON B, THE SLOPE OF TRE EQUATION? (Y/ 


WW+0 

+(WW='N' )/D3 

Been THE HYPCTEEST ZED SECEE. 
< 


¢ COMPUTE UI'S 
U+Y-(BBxX) 


A<'> t 


a CALL SPMANP TO COMPUTE R AND ASSOCIATED P-VALUES. 

D*X SPMAIP 

=(Btt3>9)/£44 

143 TPV+2xD(2] 

>(Pys1)/b18 
Li9:'SPEARMAN''S R Equals: 1 (43D(1)),O27CNL 

' THE P-VALUE FOR HO: 1 (sBB),! bERSOS 

a : ,(3CHAC1 23), (CsaB Gos (4sDC2] ) OTCNL 

‘THE P-VALUE FOR THE yo ~SIDED TEST OF HY a dorubers ‘Ss , (43 t2xD(2])),07CN 

Q IP USING THE NEW REGRESSION EQUATION BASED ON MEDIANS, EXIT HERE. 


L3:+(FF= 18 
LERESS ENTER WHEN READY.' 


ONO FWNRPOWOSONMAWUEWNHROWON MNMNEWNHOPONMOMDNMMNEWHOPROWDDNIMDNEWNHHMOWWONIMDW FW OLAw ww 


Lm_A Aj 8 «8 A RAs el A eM A ht 


NHOOD WDHOHOOMUGIGNIMNMNUMMEEE FEF EFF EWWWWWWWWWWNNNNNNAD? 


Tn mamnNlM ewe eR Una eee de 


eZ 


CO TEE EAT IIIS Eo 


mo 
sy 


+N1 
a COMPUTE CONFIDENCE INTERVALS ON 8B 
Eee eee YOU LIKE A CONFIDENCE INTERVAL FOR THE-SLOPE? (Y/N). 


WW 
> (WH="N")/N4 
B20: CCCLNEUE -S 


AQ CHANGE Bocce VALUE TO ALPHA 
ALPRA+(i00~ -CC )+ 
Q Eis iG NORMAL APPROX. FOR CONF. INT. OF LARCER SAMPLE SIZES 
> 
Q COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 
cbEert2is a 3 

Li cael 


INDEX POSITION OF VALUE IN CDF s ALPHA 
"INDEX*(+/ (CDPSALPHA)) 
+(INDEX>0 )/Z6 
INDEX+1 


>L6 

Q COMPUTING CONFIDENC 
bs: DENOM#¢ (Nx (Ho 1 9x COROT 
mALPSA©DENOMx (| (NORMPTR Aas 


eer : INDEX] 
aa Cad CONFLa Z 

NN<«1 
SLOPES#14CI 
Bael CON. TALPHA)+#2) 


a 0)/Z20 
393 >SR<r (1+((NN+TALPHA )4#2)) 
+(Sas (sSZ0PES))/521 


SH<cSLOP=S 
‘tink pSbTGinen aSOUELOEUG® TNEEPTAL 7Qhfp THE SLOPE OF | 
BSLOPESTAR?), aus < | (SeerorEstse.),!' ).',G2CNE 
I DEESS ENTER HEN 2EADY. 


Ren als USING NORMAL APPROX. W/C.F. 


*Q 


L 


pas 
eae Sein 
Fatvaat 


tr 


410.110.0120 10.010: 0D: 60 £9. 1D WD 00 00 00 00 00 00 00 00 020 0108 NININININI SSIS II 
OW ONDMN FW OM ONHMNFE ONHOW OND FWNHOW 


2 O)l- 4 AAA A 


ome 
a f- BRCUTSIOE THE C. ve CALCULATE NEW EQUATION BASED ON MEDIANS 
47 ( (BRSLOPES (221 )A (383% OPES UGE I T/NY 
ORD a 
Eek EA 
Y<Y 


fa, IF THE SIZE OF SS IS EVEN OR ODD FOR FINDING MEDIANS 


"9 ¢(2Yi5=0) S03 7S 
MPUTE MEDIAN FOR ODD CASE 
, eainaiaiant 2)] 


Q MPUTE MEDIAN FOR EVEN CASE 
Sa: B+ (SLOPESt Wis ee see) o eetoe 
AQ O THE SAME F THe 2 AND JF VECTORS 
S2:2((21N)= 03/33 

YBAR<+Y pote ane a 

ae N+1)4#2 
S3:YBAR+(YC(N#2)]+700(N Nee 

XBAR+ (XU(W#2 )J4X0 0 ee 2)42)]) 


E NEW thPeRcePr [ae 
Our: A+YBAR-(BxXB 


aay ‘ 
fae LEAST SQUARES ESTIMATOR OF B LIES OUTSIDE THE CONFIDENCE INTERVAL. ! 
, DISCARD THE LEAST SQUARES E 7 AND USE:', Pee 


Y , (3A (3B),1X 
; THIS EQUATION “rg BASED on ree MEDIANS OF THE X AND Y DATA ,AND TH 
MEDIAN OF THE TWO-POINT SLOPES CALCULATED FOR THE CONFIDENCE INTERVAL O 
a “ALLOW USER TO DO SOME ANALYSIS ON NEW EQUATION 
DO YOU WISH TO ENTER SOME X VALUES fO GET PREDICTED Y''S§ ! 
'FROM THE NEW EQUATION? (Y¥/N).' 


FF<0 
+(FPF='!¥!)/22 
FR+«!yY! 

ae 


CARA AM A jk 
iB 
wy 
rr’dy 
WY 
my 
= 
tH 


Pap dpb pa pd PD pa PAP SPA PAPAPDPAPAPAPAPAPAPAPAPAPAPAPAPSPDP DEDEDE dG ap ay 


FEWWWWW WW WWWWNNNNNNNNNNHPHPHPHHHHHEFEOOOOOOCOOO 
4 


POWOSONAM FF WNHHrPROWOWMAMFEWONHPOW ONMMNFWNH OW DO YVNMMNE Ghose 


p> pt pt pt pe pap 


v SIGN; A; jPVAL ; X:MO:N: CDF: ALPHA; 4ST Y:AA;BB;CC;3;DD;PV;PVI:NNN;KPOS;ORD 
ORD 54° kabpn -2.24-QUAsH Be 

225.08 S FUN Fai is BEF ah” ObbT Mipys éTéN TEST TO CALCULATE THE K 

STArTIstte NGaLO AND TARP TD ENCE INTERVAL AS A TEST FOR MEDIANS. 

THE LAS P bprr ION peep DISPLAY A TABLE OF CONFIDENCE INTERVALS OF ORDERED 

STATIS TCS WITH CONFIDENCE COEFFICIENTS. 

SUBFROCRANS CALLED BY THIS FUNCTION INCLUDE: BINOM, NORMCDF, NORMPTE, 


>NU 
ae SIGNHELP 

MENU CHOICES AND_ROUTE FOR PROPER ACTIONS 
0) Nu: C+CHOICES PACEDMENU SIGNQBI 


MMO cOoOnIMin Ee WhoP 
bDDDD 
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WWIII IT EIEIO YE 


FE FFE FE FF FEWWWWWWWWWONNNYNNNNNNNPPPPPPPHBE 
WONTOUN FWNHOWDNDA FWHPOWOYNDMNEWNHOWMYMM FWHP 


ke a 
D 


DOONOMNFEWNHPOWDD WMO EWNPOWDDOWNMHMWNEWNHROWDM NDMNFL WNHPOWOYNAUMEWhrRoO 


Low 8 8 ef Uk lk ll nk rg Ol 


1010: 10:10:10.0 1D 11D WH DM WW ODWDWONNNNNNSNNSNOMD DOMOOM BMONDMWMnMnMMMnninwne 


TT TTP Ero tote 


SAAQY 


A INPUT DATA FOR SINGLE SAMPLE CASE 
L8:AA+l 
se 1 


D<X-MO 
oa 


A 
L9:AA<2 
R<INPUT 2 
« 


PAIRED SAMPLE CASE 


COMPRESS D TO REMOVE ZEROS 
RECORD LENGTH OF A AND ASSIGN TO N 
KEEPING TRACK OF POSITIVE SIGNS 


aa 
£11:A«(D+0)/D 
R 

N+pA 


a 
XPOS<+/(A>0) 
+(N>25)/NCRM 
LVAL]|eZINOM N 
= (CAPOS $>0)/P1 


PV I-11 


+P2 
oui es, 
ee J 


tot 5S een s)- ae W/ GONTINUIZT? CORRECTION 


=s— 


ODF 2 
PYI+1= (NORNGDE Lo 
A PARED SAMPLE TEST GO TO £17 FOR OUTPUT STATEMENT 
6 :Pv3e2x (17 (PVF IO) 


aes 
N5:PVM+(3,1)0(PV,PVI,PV3) 

s COMPUTATIONS ARE BASED ON A SAMPLE SIZE OF : (aN 4! »OTCNL 
‘THE TOTAL NUMBER OF POSITIVE SIGNS IS: nue TCN 


qo7217 
THE Be vaLoe POR Ho: M = (3MO),' VERSUS 81: si cca. ae ook 
oHO), i TS: ', (SsPVM(c-13 433 sbrent 


> 
L17:'THE P-VALUE FOR COMPARING THE MEDIAN OF THE POPULATION OF ' 
‘DIFFERENCES TQ THE RYPOTHESTZED MEDIAN, * OTCN 
'HO: M(X-Y) 310 VERSUS S #1: M(X-¥) |, (sZocrIcl(c-6)3:1]),' ',(aMO),! 
Is: (Sepyctee =5 24348 
Lis: twoutp YOU LIKE A SA DP TD ENCE INTERVAL FOR THE MEDIAN? (Y/N).',O7CNL 


*(BB='¥')/L16 
*+QUANT 
A INPUT SIZE OF CONFIDENCE INTERVAL 
Li6:CC+INPUT 5 


ALPHA+(100- eeunece 
ees 
OMPUPING CONFIDENCE INTERVALS BY EXACT P-VALUE 


* CDP*BINOM N 

Q INDEX POSITION OF CDF FOR ALPHA + 2 

IU ae 

>(B> 

Bel 

>+SKIP 

2 COMPUTING | CONFIDENCE | INTERVALS BY NORMAL APPROX 
NORM1:KALPHA+((0.5x (NN 5) )x(NORMPTH ALPHA) )+(0.5xNNN)-0 


3 ROUND Aa LCHA Sonn 7 VESRES? INTEGER AND INCREMENT 3Y SE 
SeURALC HATS 
ae SINGZE SAMPLE CASE CO TO L7 


a 

Shr i (C=2 004.507 i 
3] CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 
DS ORDD DDT ADD 


¢I+ORDD(BI ,ORDDL(TT- (B-1))] 
Cc) CONFIDENCE INTERVAL FOR THE MEDIAN OF THE ' 
' POPULATION’ OF DIFFERENCES Is:! 
ae t (sCIC1]),' S MEDIANCX-Y) sql e7021)," JipezevE 
A CALCULATE AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 
L7:ORDX+X CAX] 


YY<opORDX 
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PREM QOO000 00 oo 
WN HRO 


b> papa papa pe 
OWOIAMN NKDEWNHKPO WOWONDMWNVE 


oe eo a | Crews eee oer 


b>p A pepeea 
hob 2r PA} 
Lut. Nh 


a sk ws 


WON NONNYNNYNNNPRPHPRPHPHPHPHMIAR OM UMMN FE wpe 


FOO BONMAMWNEWNrPOWOONMMN FWhorOr« 


rwwy EE eee 
Lm Leen A A 


MRRP MPHPHPIiDO~IOuWt FWP 
DO IMM FWMH Oswalt ALK Bw 


a a a ok oe oe ee ee ee oe om | 


I“ORDX(B] ,ORDX((Y¥-(B-1)}1 
$3 (3 é23? : CONETDENGS CE INTERVAL FOR THE MEDIAN OF THE POPULATION IS: ' 


f (SQM MEOTAN <= ', (sC7T(21) 
QUANT: WOULD YOU LIKE CONPTDENC E’ INTERVALS FOR ‘4 BSEettrep buaheie se (Y/N 


2 


- 


WW+ 

TOWW = 20 )/BL 

>i 

ee LE SIZE OF THE SAMPLE.' 

31: 'ENTE EE DESIRED QUANTI“ZE; FOR EXAMPLE: ENTER 20, FOR THE 207TH QUANTILE. 


” 


UA<O 
=((@ga0)v (qua>100))/E4 
NNN QUANC QUA 


t_ uxexx THIS TABLE GIVES CONFIDENCE COEFFICIENTS FOR VARIOUS INTERVALS ; 
'WITE ORDER STATISTICS AS END POINTS FOR THE ',(sQUA),'TH QUANTILE.',OTC 


NL 
Peeoo oNeen wae f£EADY.' 
53+) 


N 
eee CeeeuGaNieere ZameeeMUST Sie seTWeeN 0 AND 200; DAY ACali.' 


CK LA ot Ll 


V gO Ate ee Caan Bases PV 

A PATS FUNCTION COMPUTIS TYE SPZARMAN a Seal lsisG AHICd MEASURES 

a as ee OF CORRESPONDENCE 3ETWEEN 2ANK-NGS OF TWO SAMELES . Tee 
avVacds GLVEN 205 O2STiING ONE AND "WO-SLDED 32 POTRESLS OF es 
a SUEPROGRAMS CALLED gi Tais FUNCTION INCLUDE: TicS, @Pil2esKk, sPEdi 
a INPUT, SPAPROX, INTERP, AND THE VARIABLE PMATSP 
ale DISPLAY MENU a INPUT DATA. 


ee SPMANOES 
+(B=1)/81 
oe wh hope 


Bi: TR+INPUT Z 


<14h 
Ye (gst)sa 
+€ 
A CALL _SPMANP TO CALCULATE TRE STATISTIC AND ASSOCIATED P-VALUES 
A+tX SPMANP Y 
>+(Al1]>0)/51 
Guae! InbrRECr’ 


ie apne! DIRECT! 
L2:PV+2xA[2) 
sieve IES 


PV<+1 
L3:'SPEARMAN''S B EQUALS: (43A(1J),Q7CN 
' THE P-VALUE FOR HO: NO assocrarion’ EY fers VERSUS! 
Hass (sCHA) ASSOCIATION EXISTS IS: (4a 421) Orch 
PEE P-VALUE FOR THE ro. Cua ed TEST OF HYPOTERESIS IS: ! aie a), SeoN 


L 

‘PRESS ENTER WHEN READY.! 
BB<O 

ia 

V 


v WISIG: A:B:D:E:F :PV2 ; $213 ;Z;DEN;NUMZ ;NUMZ1 ; PVAL:X3MO3N; TPLUS;CDF;TALPBA; AL 
uid £i4aiz :BE tC bd; >PV;POS;fPOS;NM:PVI:TPOSi:NNN;C;PVM:PV3;2:Q:N 
A ruts Citi iE V bse ie ben WILCOXON SIGNED RANK TEST TO CALCULATE THE TPLUS 
a STATISTIC oe =VAr Be, AND CONFIDENCE INTERVAL AS A TEST FOR MEDIANS. 
a SUZBPROGRAMS CALLED’ BY THIS FUNCTION INCLUDE: TIES, WILP, NORMCDF, 
3 NQaMe re CONEW, INPUT? AND THE VARIABLs PMATRIX. 


Wi MENU WILZZL> 
MenU CHOICES AND ROUTE FOR PROPER ACTIONS. 
vu esgegrery, PAGEDMENO WILOBI 
+(C=6.7,8,9)/L9 
+(C=1 tha? 
MENO MAINOBI 
> 


A INPUT DATA FOR SINGLE SAMPLE CASE 
L8:AA<l 

eae 1 

Cr 85/L16 

MO<INPUT 3 


85 


D+X-MO 


A 
Cat 
ORDER _ POSIZTIVE SIGNS TO COINCIDE WITEVEROPER FOS ito seen ao 


ee ee ee ee ee ee i on on en ie en een me 
WMP FEEFE LES FEWWWWWWWWWWHONNNYNNYNAN 


0 
A Set 
alg PAIRED SAMPLE CASE 
3) 29:AA<2 
u)  R+<INPUT 2 
5 Qei tk 
6 +14(Q+1)4R = 
q Y<(Q+1)+R = 
8} DD+«xX-¥ 
9 NNN+o0 DD 
0} +(C=9)/L16 
1) MO<INPUT 4 
2 De(X-Y)-MO 
31 a COMPRESS D TO REMOVE ZEROS 
u} 211:A+(D=0)/D a 
: * ae, RECORD LENGTH OF A AND ASSIGN TO N 
+6 
7 Arona KEEPING TRACK OF POSITIVE SIGNS 
«(4 
3 "Ee TIRE TRE ASSOLUTE VALUE OF A; ASSIGN 70 3 AND ORDER 3 
<_ 
i] 
3 "POS*P0Si3( 1a 71 
4 CALL FUNCTION TO BREAK TIES 
5 “Bet TIES B 
6 CALCULATE TPOS BY ADDING ACROSS ALL POSITIVE VALUES OF E 
33 "ppose+/ (POSE) 
3) 7TPOSi<«TPOS p ‘ : 
91 A GIVES SIZz OF LEFT TAIL OF PROBABILITY DISTRIBUTION 
O} NM<(LC(+/1W)4 _ 
11 4 6 70: “STATEMENTS BASED ON LENGTH OF VECTOR 2 
°53) +(N>9)/ZL3 
33) A GENERATE NULL DISTRIBUTION FOR TPLUS 
54} FeWILP N 
55] a IF TPOS FALLS TN LE: PT HALF OF PROB DIST CALCULATE PVALUE AS NORMAL 
56 bi: >(TPOSS (NM-1))/TR 
57 Tae RVISE a OSE THE NEGATIVE T STATISTIC 
(33 NTPOS*(+/.H1)-2 DOS 
r59) a IF TS FRACTIONAL USE BOTH "HE INTEGER ABOVE AND BELOW AS TPLUS 
50) PEC: re Teas he / NON 
61 TPLUS«| TPO Ss 
62} Pi: PV<1-((FCTPLUS]+FC7TPLUS+1] )#2) 
63) PVI«(FCTPLUOS+11+F(TPLUS+2] )#2 
64) +CHECK 
65) NON: *(PP0S>0)/¢0 
66] PV< 
B71 +P2 
68] CO:PV+1- Et <TPOS) 3, 
69) P2:PVI<FE 
70 PPAR ei Te 1) 756 
71] PV2«PV 
72) PV«PVI 
73 PYI*PV2 
> 
75) a COMPUTE NORMAL APPROX. W/CONTINUITY CORRECTION FACTOR 
76) 23: TRAP<(Nx(N+1 
77) NUMZ<(TPOS+0.5 “TRA 
78) NUMZ1<(TPOS-0.5)-T 
79 DEW<( (Wx (V1 35 ((2xH +1424 x0. 5 
80 Z«+NUM 
81 oT MUMES DEN 
821 a COMPUTE STUDENT T APPROXIMATION WITH CONTINUITY CORRECTION FACTOR 
83) NOM+« (7P9s -TRAP) 
84) DENT< iB EIa235 (N- 1) 2s Cf (UM. 0 5 )*2)+(N- 1)))s *0.5 
85) DENT1<(( (Nx (DEN*2))+(N-1)) (NUM+0.5)*2)#(N-1)))*0.5 
86 TC+(N .5)+DENT 
87 Pete NUMGO § )tDENT1 
881 a COMPUTE AVERAGE ge TC AND ZC 
89 z(TPoss( (4/1 )42) )/ SECOND 
gO ee MCDF 23) )+(1-C(N-1) es TC1) 742 
32 PyI« (NORMCDF Z)+((N-1) PDISTN TC) 
> 
»23) SECOND: PV+((1-(NORMCDF 21))+((N-1) TDISTN 1TC1))+42 
gu) PYT~( (NOR ORMCDF 248 +(1-((N- HG TDISTN PON) 22 
(95) 56:2732+2x(1/(2V,2VI)) 
ee 3(2V3s2 \/NS 
3 / 3-1 
ee NS: Pym (3 4 top (Pur PV,PV3) 
99}. 'COMPUTATIONS ARE BASED ON A poAMeLE SIZE OF ',(eN),'.! ,OTCNL 
100 'THE TO A L SUM OF POSITIVE RANKS IS: (sTPO$1) CNL 
102 oe IF PAIRED SAMPLE TEST CO TO o be POR Sain STATEMENT 
103 ‘THE p. VALUE FOR HO: M oMO),' VERSUS 81: M ',(sLOGIC({C-1:1]),! ',(s 
Be w0)8 IS:', (usPyM{(c-1; 1)), ‘arene fekSUS 
105) £17 THE P-VALUE FOR COMPARING THE MEDIAN OF THE POPULATION OF ! 
106 ' DIFFERENCES TO THE BYPOTHESTZED MEDIAN, ',OfCNL 
1073 'HO: M(X-Y ERSUS H1: M(X-Y) !,(sZOGIC(C-531]),' ',(sMO),! 


TS: ae seMec- 5; ae ArcNnz 


/ 


86 


PAPA RAED PDE papa papaya 


FE EWWWWWWWWWWUNNDNNYNNYNYNNEHBHEEEEEAEEOO 
HOW KO IMDWFLWNEOWWOIDMEWNYNHROWOIMNUMNEWNHOWO 


Ce A MU LK 4 ANd AAA nl XA ok 


PAP PSE PARP PAP AP AEAPA PAP dea G ry AE DE OG Dp AP apa Papo pa EO 


Cw we OEE Oe aE eee eee 


FFEF 
Mm Fa) 


od oR ee YOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (Y/N).' 


*(BB=!2')/Li16 e 


Bia: CoeINeUT 5 
ALPHA+(100-c¢ )+20 


A 

+(NNN>93)/LY 
CDF«<WILP NNN 

a (core INDEX POSITION OF CDF FOR ALPHA + 2 
CIF<+ - 


0)/CD 
PALDBL=(+/ (COPEALPEA)) 
* (TalezAzo )/ JUMP 
TALPHA+ 


oJ 
a COMPUTING eee Sau INTERVALS BY NORMAL APPROX. W/C.F. 
Du: TRAP+ (NNNx (NNN+1 ) 
DEW< ( (NWN (NNNS2 )x ((2xNNN )+1) 424 )%0 » 2 

gare, ORS AGEAR It iaaAr-O. > 


va SA 
A ND Ear oae DOWN TO INTEGER VALUE AND INCREMENT BY ONE 


ZOU 
ALPHA+LTALPHA+ 
Q F_ONE SAMPLE CASE GO TO £7 
Joma >(C=2,2,4; 375 
- CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 
L5:CI+TALPHA CONFH 
7 bees pC ONFIDENCE INTERVAL FOR THE MEDIAN OF THE ! 
‘POPULATION 5¢ ’ DIZFERENCES TS:' ,OZCNL 
7 (8 cre i) ee : *MEDIAN(X-Y) s ',(ecr(2]),' )'!,O2CNL 
\ >RESS ENTER WHEN BEG 
Ld 


LO ° NORMAL APPROX. FOR CONF INT OF LARGE SAMPLE SIZE 


23 
~VG 
2 one caLCUL ATZ AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 
tad f ee ail G wv 

Be wecy, nONPIDENCE INTERVAL POR TEE MEDIAN OF THE POPULATION IS:',, 

epee t i): S MED Pa seeCOETL 2). )' (ET CNL 
(PRESS ENTER WeeH ‘pEs 
ago 
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APPENDIX E 


MAIN PROGRAM LISTINGS FOR MAINFRAME COMPUTER WORKSPACE 


a a NT a a BE TN TN OOF 


SNNYINMNOMOMHOHDMHD GO DAnnNnnnNnNAKN eS EEK EEE ES EWWWWWWWWWWNNNNNNMYUNNNEEPEt 


£b2t+$4 A190 00~700in FWP 


FWNPOWD DYNAM EWN POW DONMNMNEWHPOWONMHMNELWNHPOWOWNMNMN LEON POW ONAN FWNHPOWO~IWDUW FWDHYOUsA «A «wn 


ie lA ek 


itt 2 a ai 8 ke ek ek ALA 


Vv RENDALL 3A sAA BiB -BX; XsBY 5s CXiCEsD 2D: :DD i Dx; DY; DXY:S;POS;NEG:XX:YY3:N;DEN;NN;NU 
nhs "bon ction ‘one pin wes’ ree ig RD ASF b parAriszic WHICH IS A MEASURE 
on $ P-VALUES ARE GIVEN FOR TESTING ONE 
AND 7 DOreEsre. FoR” 30 PSSGerArTON VERSUS ASSOCIATION. 
SUBPROGRAMS CALLED 3Y THIS FUNCTION INCLUDE: TIZS, TIESK, KENDALP, 
INTERP INPUT AND NORMCDF. 


Od 
ny 
ho 
4 
WY 
dS 
om 
se 
I 


ppDDD 

| 

= 

Oo 

—~ 

Hy 

| 

See 

2a 

‘ay 
ks 


eee 2 
Fie cQe4 ar 
¥<(Q4+1)+2 
a YCAX) ORDER Y IN INCREASING ORDER OF X 
ke 
_ Poy ORDER X IN INCREASING ORDER 
<< 
Ip Rr COMPUTES CURRENT RANKING OF I 
a= 
Q NOW ORDER Y 2ANKS IN INCZEASING ORDER 
D<A CAA] 


A IF TIES GXIST IN EITHER X OR Y RANKED VECTOR USE MID~-RANK METHOD 
DD<Ter res 2 
AATIGIIECS 2 
R FIND ORICINAL RANKING OF Y WITH TIES RESOLVED 
Re ne 
=0 


a COMPUTE NUMBER OF DISTINGUISHABLE PAIRS 
Sopot 
<p0 


AA«< 

A SOSITIVE ONES COME FROM A RUNS UP CONDITION: NECATIVE 1 FROM RUNS DOWN 
ea Lp a SCORED FOR TIES. MULTIPLY THE RESULIS FOR EVERY ELEMENT AND SUM 
Bho CCA) > Ca - 

CX« xxldal< AAVXX))xC 1) 

ate 

BY<(YYCAAJ>(AAVYY 

Some ever eSoat Tiare. 

DY+BY+CY 

DXY<«+DXxDY 


weg: (Dx r<¢ 
NEG+(DXY<0)x("1) 


SaaS ed /L1 
ie y) SUM FINAL VECTOR TO DETERMINE S 
<< 


S 
A OBTAIN THE NUMBER OF TIES IN EACH VECTOR USING THE TIESK FUNCTION 
O<TIESK B 


ULATE THE B STATISTIC INCLUDING THE CORRECTION FOR TIES 
"pest (NN SU )x(NN-SV))x*0 


A 
> (W>12) NORM 
CALL KENDALP TO CALCULATE THE RIGHT TAIL OF THE CDF OF B 
* peRENDALP N 
A CALL ae eee P-VALUE BY INTERPOLATION 


>(PVAL= 1)/L3 
GeO. 2 


ge €ALCILATE 2? VALUE GSING NQRMAL APPROX. 
NORM: NUM<(2*AT)x( C2xXWN )x0.5) 

DEN< (2x ((2xN45))x0, 5 

Z«NUM+DEN 

PVAL+1-(NORMCDP Z) 

A IF B a POSITIVE PRINT OUT DIRECT ASSOCIATION. 
L3:+(f>0)/L5 

CHA«'!INDIRECT' 


+7 
L5:CHA<'DIRECT' 
L7:PV<2xPVAL 
+(PVS1)/L8 


PV+1 
L8:'KENDALL''S B EQUALS ',(4sT) 
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SININISNI™N 


TET EEE Oe eee 


FEF FWWWWWWWWWWNNNNNNMYINIEAR PHPBB oo “SIOcn FE Gh ps 
WN ODDO ION EW ODONMMN EWN OWONMM FWor Ob 44 Laas as 


iti AN ll 


HAINAOAOMNMNNnMNNMNNMUME EEL EL 
Onin FE 


£ 
PP ONHWFWNPROWDONA WEWYROWOBIMMEWNROWWY 


EE OWN eee 


SININISINSININI NI OOH Mo 


a oe oe ee ee ee ee es ee ee 


,2ae P-VALUE FOR HO: NO ASSOCIATION EXISTS VEESUS' 
HI: ',(€8Ca@),'ASSOCTRTION EXiSTS IS: ',(4sPVAL) 


' TRE P-VALUE FOR THE TWO-SIDED TEST OF HYPOTHESIS IS:',(4sPv) 
- 


V KRUSKAL:NUM:DENOM:A K:AA:BB:DD3E:F:N:0F3P:2VAL:2:SOFR;SRE:;TSOR:CHA 
Q THIS ®UNCTION eanbohee? hee EB KUSKAL-WALLIS 7287 STATISTIC a WHICH IS 

a A MEASURE OF THE EQUALITY OF K INDEPENDENT SAMPLES. 

a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: Treo. TTESE INDEXPLS, 

a FDISTN, INTERP AND THE VARIABLES PMATKW20, PMATKW31, PMAPKW33 

a SuATEW SS, PMATKW41, PMATKWU2, AND PMATKWU3. 

— 

Bi; \ENTER THE NUMBER OF POPULATIONS TO BE COMPARED (MUST BE GREATER THAN I 


K+ 
a K<3)v((11%)=0))/21 


(con )>1)7E2 
A INITIALIZE VECTORS E AND F AND VARIABLE C 
em Fr<e50 


“PHS Pes ENTERING THE SAMPLE VECTORS AND STORING THEM 
«! 
Li: 

t ENTER * your TCECHA),' SAMPLE .' 


De 
ee NERS 


oer peo SAMPLES AS TREY ARE =NTERED AND STORE THEM IN VECTOR & 
Av 7 3E+e, 
a BECOAD THE ZENCTES OF TBE SAMPLES AS THEY ARE ENTERED 


PtP, 0 
CRA+'NEXT! 
ees (ha 7 E24 
S(Cex3/ 21 

= 

ee RECORD SIZE OF ALL SAMPLES WHEN COMBINED 
—— 


aA ORDER SAMPLE SIZES LARGEST TO SMALLEST 


OFP+FCYF] 

ECAF} ORDER COMBINED SAMPLE VECTOR TO BE USED BY TIES FUNCTION 
+ 

A CALE INDEXPLS TO INCREMENT INDEXES WHEN TIES OCCUR WITHIN ONE SAMPLE 

AA«F INDEXPLS E 

A CALL TIES TO BREAK TIES BY MIDRANK METHOD 

a Lies D 
“4 


THIS LOOP CALCULATES THE H STATISTIC 


=) 
L2:C«C+1 
SUM OF RANKS | FOR EACH SAMPLE IS CALCULATED 
OTe PEG uae’ 
A CALCULATE SUM SoplaRe SQUARED DIVIDED BY QYHE INDIVIDUAL SAMPLE SIZE 
SReCe eae ERO 
STORE EACH CALCULATION 
SOFR*«SOFE,SR 
a SK b2 


" PSOR++/SOFR 


H (TSORx (124 (Nx (N+4)))) = (3x (N+1)) 

a RECALCULATE H WITH CORRECTION FOR TIES 
A«T7I Sok E 

NUM+(+/(A*3))-€+/4) 

DENOM€Nx ((N*2 )-1 


D 


D 


SUM ACROSS ALL SAMPLES 


gi EF Shem one Denon) STATEMENTS ENSURE PROPER PROB. 
+ (BELL 22 0 /OUTSUE SE ala 
> rate 
+(K=3)/TP 
>+((A/(OP= 3 3 3 3))v(OP[1]>3))/PAPPROX 
+(A/(OF= 2-1 1 1))/QUTPUT 
(Rau a (OPT 7=3 $47 DE3 
s((a/(OF= 3°72 2°2))v(A/(OFs 2.2 2 1) v(a/(OFs 3 3 4 1))v(A/(OR= 33 2 1)) 
hes a 
>(Z=u)/Puz 
Oe ota (or= 3 1 1)))/0urrUr 
= VCA = 
saat SCALES, say ovearcare 2 3 
= VLA eo ViA = = 
Beene 6 4s Os Sty ee 2))v(A/ (OFS Suede sy) 7 ES 
CALL APPRO RIATE VARIABLE ACCESS CDP 


Pi: P«+PMATKEW310(N-5); : 
> 
P33:P+PMATRW330(N-5)33] 
>PM 
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P24 3 Ppa TNS 
oe 
rab ey ooh 


Pua: P<PMATKWU3C(N-7); 
A CALL aioe LO iiicutare P-VALUE BY INTERPOLATION 
PM: PVAL< I aa 


TER 
* (BVAL= 1)/ OUTPUT 
a CALCULATE P-VALUE USING THE ¥ DIST W/ONE LESS D.F.IN DENOM APPROX. 
te iBet (We K xB tC ae x ((N-4)- ae 
PVAE~1=( (Kal) (Nk eer 


>ia 
OUTPUT: ae ZHAN 23% 
LS: 'THE F STATISTIC EQUALS: ', (48H) 


THE P-VALUE FOR HO: THE POPULATION MEDIANS ARE EQUAL VERSUS ! 
81: AT SEAST TWO POPULATION MEDIANS ARE NOT "EQUAL ITS: 4), (sPvac) 


> 


HBE&HOR: YOU MOST 2NIZR A SiNGis INTECER VALGE Catarina sah ace 


OT EIT ETAT 
HfrREWIDDWDO WO WOWWO MW OWWOWOWMMDOMOBOMMO™~ 


hi 


FW PH Ol A411 A nn 


Lt t wad 


OO BDOOWOAMMWEWNYPrOWONHMHME WOH OW 


NNWHIT :N3M:PV2:2:3:C:G:MM;NN:PX;U;NM1;2;NU:PVAL + NM: NUMZ:NUM241:DENC1 32 
-21:ALPEA: COP; DOS Beet DY - OF Ga Loas | jee icciOts U2>DV:NNLINN2:2VI: DIFF 24:66 
>PYM:PV3:D:9:8:DEN:DENC:TC:TC1: 

THIS FONCTION USES THE SUM OF VO ANKS PROCEDURE [0 CALCULATE THE 
MANN-WHITNEY y Srarrerte USED IN COMPUTING THE P-VALUE ¥OR THE TEST 

0 AND SCALE. THE C.1. FOR M(Y)-M(X), THE SHIP 

LOCATION, IS ALSO COMPUTED. SUBPROGRAMS CALLED INCLUDE: Pres, TIZS2 
pepe S;, VARMW, MANWP, INPUT, CONFMW, NORMCDF, Nop NORMPTH . 

t ? 


'DO YOU WISH TO COMPARE THE MEDIANS OR VARIANCES OF THEE POPULATIONS?! 
ca ENTER: 1 TO COMPARE MEDIANS; 2 TO COMPARE VARIANCES.' 
ale 
iV Erb y/ eee 
> = 
E TEST TO COMPARE VARIANCES REQUIRES THE TWO POPULATION MEANS! 

‘OR MEDIANS TO BE EQUAL. IF THEY DIFFER BY A KN OW WN AMOUNT, ! 
'THE DATA CAN BE ADJUSTED BEFORE APPLYING THE TEST 
B3:! ENTER THE DIFFERENCE OF MEDIANS (M(X) - MCY)) OR 900 TO QUIT). 

= 
>ep ee 
=_>> 
V 


DDDDD 
ry 
sty 

Q 

tiOr 
hy 
w hy 
© 
= 


DIFF=900)/0 
Tae z NOLL BRYPOTHESIS STATES - THE POPULATION VARIANCES ARE EQUAL; V(X) = 


WHICH ALTERNATIVE DO YOU WISH TO TEST?! 

EGE : 1 FOR H1: V(X) < VOY); 2 FOR Hi: V(X) > VCY); 3 FOR Hi: V(X) 
1)A CD42 )AC D243) )/E4 

7 


gyn NOLL BYPOTHESIS STATES - THE MEDIANS OF X AND Y ARE EQUAL; M(X) 


WHICH ALTERNATIVE DO YOU WISH TO TEST?! 


: Wey 1 FOR H1: MCX) < MC(Y); 2 FOR H1: MCX) > MC(Y); 3 FOR #1: MCX) 
((D21 a (D#2.)a(D23))/Bu 


CLL TN DFC LLL Eee 


WNNEE EF EFE EFF FwWOWW WWWW DNYNDY NHNNDD NHNNHPRPRPPPRPHPRRH WOW OnMFEWnNH 


a ENTER DATA VECTORS 
eee ER « DATA (MORE LZHAN ONE OCOBSERVALION IS VREQGULAED). 


+((p9N)=20)/E1 
'ENTER Y DATA.! 


was 


aren Oe ae INVOLVE VARIANCES ADJUST X BY THE DIFFERENCE IN MEANS 
«<N a= 


FWNRPOWDONIMMW FWNHrF OWON ANEW HPOPOWONHMHMNLWNYH Ors» «4121 


eee ee ee ee ee ee ee ee ee eee ee eee eee ee eee ee eee ee eee eee ee eee ee ee ee ee 
= 
ay 


rrr ww oo le 


5 "A vou CONCATENATE X AND Y SAMPLE VECTORS 
+ 
f ee oN DETERMINE SIZE OF X AND Y VECTORS AND ASSIGN TO NN AND MM 
+ 
9 (Mae M 
0 ; fone te SIZE LIMIT OF LEFT TAIL OF NULL DISTRIBUTION 
a "NM+ (NVXMM #2 
2 NM1+ NM 
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bnOimnina 


SOME ON rPROWONHUMW EMO HROWOYNMU FW 


TE AE eo oe 


00 DOBYNS NN DHHMMHAMAHHNMHuNYG 


FOWtROwo 


Ae kA ee Al A A A 


“wae 


B22 bpd kt kt (1 LO WW WL) WW) WZ 00 © 00 CY © 0 0 
DNRNN NNNYNBPEPBPPBPBRBBPHOOOOOOQOOOOWMOUMUEWNHOWMOQHU 


SOME WHPOWDNHAMEWNPOWONMNU EWN Os 


LO EE eee ae 
ht Me tk ke hl 


PAPAPAPAPAPAPAPAPAPAP APD pay apa 
FE FWWWWWWWWW Wh 
NROWONHDMWNEWNH Ow 


ORDER A AND ASSIGN TO B 


BA (AAJ 

AOeNEE s ae eae ae A 

ae 4. o9re ZL IIES FUNCTION TO BREAK TIES USING MIDRANK METHOD 
+ 


Eg" 
A oF a hl Te TEST FOR VARIANCES 
ts 
ae vane eee ZO GENERATE RANKS REQUIRED FOR VARIANCE TEST 
oe 
ee RECORD TIES IN THE DATA AND BREAK TIES IN GG 


CALCULATE SUM OF X RANKS 
35:RX++/(GL(CL1; (1NN)J)}) 
A CONVERT TO MANNWHIT U STATISTIC 
oieo ((NNx (NN+1))42) 


IZE OF X TIMES SIZE OF Y > 80; GO TO NORMAL APPROX 
A gC (NNx2 So80 5 sro 
has! hi My 
+! iV 

ama FUNCTION CALCULATES LEFT TAIL CUMULATIVE PROBS. OF U STATISTIC 
a ven ee 85544 STATEMENT ENSURES ONLY LEFT SIDE OF NULL DIST IS USED 
As asa CONVERT U STAT WHEN CREATER THAN LEFT TAIL VALUES 

ae x -J 
a F U1 IS A FRACTIONAL, INTERPOLATE P VALUE 
£33 2(G10)2 =0)/NON 


P41: PV+1-( (ietua sense] ye2 
NOW: (ro 0)/GO 
>? 
CO:PV+1 ane) 
P2:PVI<+PCCU 
CEECK i» LT ene /Lu 


bye yr 
PYI*PV2 


A "COMPUTE THE NORMAL APPROXIMATION W/CORRECTION FACTOR 
£2. eNO Z+ (070 . 
NUMZ1+«(U- vu 
HONE ono SNM ao. 5 
Z+NUMZ+DEN 
Z1+NUMZ14+DEN 
NUM+ po 7am) 
DENC+ OTe eae NN+MN- hee) yor 
DENC1+ (NN+MM-1)x(DEN*2 ) )+(NN+MM-2)) ¢ 
TC+(NUM-0.5)#DENC 
TC1«(NUM+0. 3pEDENC2 
+(USNM)/SEC 
(NOSHEO? Z)+(( iN? TDIESTNS TC ))+42 
1-(NORMCDF 21))+(1-((NN+MM-2) TDISIN 7C1)))+2 


PVI+«((NORMCDF 2Z)+(1- ( (NN +MM< 2) TDISTN PC )))+2 
- (NORMCDP Bd OE 2) TDTSIN Tei) + 
ibe WORtSHE peVT 


a 
Sun5 rn 


PV 
NS: ae. 1)9(PVI,PV,PV3) 
‘THE SOM’ OP THE X RANKS IS: ',(@RX),'. THE U STATISTIC EQUALS: ',(8U1) 


(AA uae STATEMENT FOR VARIANCE OUTPUT 

ts = 

' THE P-VALUE FOR HO: M(X) = M(Y) VERSUS H1: M(X) ',(@ZLOGIC{D3:1]),'! MCX) 
2S: 1, (4aPVM(D:1)) 


= 


L8 
VAR: PVM+(3 1 dp (PY, pvr, Py3 
'PHE 2-VALUE V(X 


<3: : CE SPUMED ola. 


Sey Olja vieovs Ai. VC. )y", (SLOGICED;11),"* VC) 


~0 
“eRe. 700 IKE 4 A CONFIDENCE INTERVAL FOR THE SHIFT IN LOCATION(M(Y) - 


BE+O 
cipB='Nt)/0 
L10:CC+INPUT 5 
ALPHA+«(100-CC)+ 
ROUTE T NORMAL APPROX. FOR CONF. INT. OF LARGER SAMPLE SIZES 
+((NMx2)>80)/25 
COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 


D D oD 


COF+P 
INDEX eapenaS OF VALUE IN CDF s ALPHA 
SP ea Gente eT aA ) 
+(INDEX>0 )/L6 
INDEX+ 


PPP ppp pppppappppppppppyp 
OORNOOOMOOMNMNMNNHnNnUINNNnELELK EEK 


ta tk ll 


ONOMNEOOrOWONAUFWMORrROWONMMEW 


i a ee ek oe eo | 


Age 
MPUTING CONFIDENCE ee OSING NORMAL APPROX. W/C.F. 
Ls: vaxeaa=COtons* VIS ALPHA ) )+NM- 
ROUND UALPHA DOWN AND INCRENENT BY ONE 
" INDEX+\ UALPHA+ 


L6:IPX<«NN TNDEX 
CI«IPX CONE MH A 

A’! .(3CC),' PERCENT CONFIDENCE INTERVAL FOR THE SHIFT IN ' 
7 BETWEEN POPULATIONS X AND ¥ TS: 


( ', CeCIC1IM,' s MC) = MG@rs ', @Cri2i Se 


= 


0. 
218 ERROR: THE SIZE OF YOUR SAMPLE IS LESS THAN TWO; TRY ACAIN.' 
ao 
3 |ERROR: YOU HAVE ENTERED MORE THAN ONE VALUE. TRY AGAIN.' 
3 
| ERROR : YOU HAVE NOT ENTERED A VALUE OF 1 OR 2; TRY ACAIN.' 


EN: | ERROR: YOU HAVE NOT ENTERED A VALUE OF 1, 2, OR 3; TRY AGAIN.' 
+(AA=2)/N2 
+B6 


7 NPSLA;N;SUMX;SUMY; XBAH; Y3ZAR; SUMX2; eee A; WW s ax bie Us 333 ee 7P3CO 3 Cie 


;TALPHA:NN;CL;SLOPES;2R:SA3 77 3 DENOM: INDEX: F?:X: 
a _'DROGRAM CONDUCTS NONPA iniraic “aiNean BECRESSION Phe Phase SQUARES 
a ESTIMATED REGRESSION LINE IS PUTED WITH HYPOTHESIS 7ZSTING 
a CONFIDENCE INTERVAL AVALLABLE FOR THE SLOPE 3. T? 3 DOES NO? Lis IN 
a THE C.I. AN ALTERNATE REGRESSION LINE IS PROPOSED. SUBPROGRAMS CALLED 
A ARE: SPMANP, KENDALP, NORMPTH, INPUT, AND CONFLR. 
+ 

A INPUT DATA 
R«INBUT 2 

eee 

ge ath eR 
*w c ASSIGN THE SIZE OF X (AND Y) TO N 

+p 

COMPUTE THE SUM OF X'S AND Y'S 
" SUMX++/X 
SUMr++/2 
COMPUTE THE MEAN OF X AND Y 

" XBAR+SUMX+N 
YBAR+SUMY+N 
A COMPUTE THE SUM OF THE X'S SQUARED 
SUMX2++/(X*2) 
A COMPUTE THE SUM OF X TIMES Y 
SUMXY++/(XxY 
a COMPUTE 'B HE SLOPE OF THE ESTIMATED LEAST SQUARES REGRESSION LINE 
B+ ((NxSUMXY )- SOME xBUHT )3 (Os SUMX2) - (SUM 
A COMPUTE 'A', THE Y-INTERCEPT 


ne (BxXBAR ) 
THE LEAST me" ESTIMATED REGRESSION EQUATION IS: ! 
oe y= ' 5.060), 0 + 8 Kooga 


DO YOU WISH TO ENTER SOME X VALUES TO GET THE PREDICTED Y''S? (Y/N). 


> GH Oe! =< 
bz: VENTER x VALUES. ' 
XX<O 


a CALCULATE PREDICTED X'S 
+f 

‘THE "PREDICTED Y VALUES ARE: ',(®YY) 

1 WOULD YOU LIKE TO RUN SOME MORE X VALUES? (Y/N). 

i= 12" )/52 
i ZOU LIKE TO TEST AYPOTHESIS ON 8, THE SLOPE OF THE ZQUATION? (Y/ 
w+ 

+(WW='!N S 

ENTER "THE HYPOTHESIZED SLOPE." 


ae 
COMPUTE UI'S 
“Usz- (BBxX) 


* GALE E_SPNANP TO COMPUTE RHO AND ASSOCIATED P-VALUES. 
=(Dtt}>0)/041 
ig aoe 2xD(2] 
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TEE EEA wo oeoy 
WOM MOO INN INYIYIYINMMHOMHOHH BMHMHuUwH 


MEWWRP Ow 


PPEPPPPPPBPPPPPPPPPPPPPPRPPHRPPPHPHP DWWWOWWWWWW MMW 


WWWW WWW WWNNNNNNNNNNPRPPPRPPRPPEHP POO OOO OOOCOO0ONHMFWNRPOMWO YUM 


CO EEE een ooo re ee a a 


cn te 
£19; ' SPEARMAN''S R EQUALS: ',(4sD(1iJ) 


Tae B= Hae FOR #0: 8 = ',(eBB),' VERSUS H1: B ',(scBA(1 21),¢8BB),' IS: 
' 3D 


t 
TRE P-VALUE FOR THE TWO SIDED TEST OF BYPOTHESIS IS: ', (4sPv) 
Le Te THE NEW REGRESSION EQUATION BASED ON MEDIANS, EXIT ERE. 


Qa 
Lave r= Tt 
a COMPUTE CONFIDENCE INTERVALS ON 8B 

ne YOU LIXs A CONFIDENCE INTERVAL FOR THE SLOPE? (Y/N). 


a Hs" ee) 
£10: c+ INBUT = 

a CHANGE ieee? VALUE TO ALPHA 
ALPRHA+(100- eae 


q ct Bee at” YOSMAL APPROX. FOR CONF. INT. OF LARGER SAMPLE SIZES 
as > a3 
COMPUTING CONFIDENCE INTERVALS oY EXACT P-VALUE 

eee N 

a” p< 7 

INDEX POSITION OF VALUE IN COF s ALPHA 

" INDEX«(+/ (CDF FSALPHA)) 

> (INDEX>0)/L6 

ail 

> 


q COMPUTING CONFIDENCE INTERVALS USING NORMAL APPROX. W/C.F. 
65 :DENOM+( (Nx (Nw1 )xC(2*N ) +5) 428 )%0.5 

eres | oka rs ALPHA) 

=f 

be: LAZeaAee’ 3: INDEX] 

o3:Cl<«x CONFLaA ag 

NN«Lscl 

SLOPES+1+CI 

BB<l CNN ae 

a )/L20 


L20:SRk<i ee ee 
= (S25 (oSCOPES } VE 
Sp aaa 

2a ew ' PERCENT CONFIDENCE INTERVAL FOR 8, ZHE SLOPE OF '! 


‘2 ree repens pen REGRESSION LINE, IS: 
- ( ', (SSLOPES(RR]),' < B < !,(8SLOPES(CSR]),' ).! 


p 
a IF B OUTSIDE THE C.I. CALCULATE NEW EQUATION BASED ON MEDIANS 
a+ ( (BBSLOPESLAR] Ja (BSSLOPESTSHI)I (0 
ORDER xX A 
beinty 
+ 
» ip id 6) oo IF THE SIZE OF SS IS EVEN OR ODD FOR FINDING MEDIANS 


A RE 
nie N= OO7S 

t6 MPOTE MEDIAN FOR ODD CASE 
apnea #2)] 


Q MPUTE MEDIAN FOR EVEN oor 
$1:B+(SLOPESC(NNS2 1 +SLOPEST C(NN#2 94501 
Bs . O THE om FOR THE X AND rECTORS 
YBAR<Y $$ 


j= 
W083) 
ee a N+1)#2 
S3:YBAR«(YC(N+2)1+Y(C((N+2)#2 ]}+2 
area Rs TE Sass 1 ETE LE )#2)] 
a COMPUTE NEW LNPERCEPT tat 
Our: A<YBAR-(BxXBAR) 


2 
E 

'THE LEAST SQUARES ESTIMATOR OF B LIES OUTSIDE THE CONFIDENCE INTERVAL.! 
| DISCARD TRE LEAST SQUARES EQUATION AND USE:' 


Y= !,(3A),' + ',(oB),'X ! 


! 
: THIS EQUATION IS BASED ON _THE MEDIANS OF THE X AND Y DATA AN 
eS el CM tHE 2VC-POINT "SlCr eo CALCULATED FOR THE CONFIDENCE INTS DRV 


ON 
t 
ALLOW _USER TO DO SOMZ_ANALYSIS ON NEW SQUATION 
“1 p90 YOU" wIs# TO ENTER SOME "x VALUES 70 GET PREDICTED Y''S FROM TRE NEW E 
ATION? (Y/N). 

Fe 
Gee" 1" )/L2 
PR+ty! 
+L1 

v 


9) 


POW OYIOMDOEWNMOWONIMDUWEWN POL tA_nww_s an 


Crwrwaw ee ea a aaa rea a aa aoe oO rn oe eo eee ee 
A kk a ll 


FEE EWWOWWWW WWWNMNNDPNNNNNNEHRHMMHPEHE BHW ONDUF WHF 


WNrOWOIMO FW po 


c 
oO 


Li 


FEEE 
NOM & 


a ek ee ee ee ek ee ee ee ee ee ee en ee See ek fee oe * dee | 


“PINDDHD BRHMOHMDMNnMNINIANIAnnIni£ & 
NOHrOWON MNEWNMHMOWONMDMNFWHHOWD 


a a a a a oo oe | 
oe oe ee ee ee ee a ee ee ee ee ee 


00 00 00 00 00 00 CO MJNIMNI~INISI 
MNFWNYrOWO IMDM FW 


qv SIGN: N3A3Ci 33D -D: EVAL Xm :MO: iN icDR: ALPHA CI: -AA;BB:CC:DD;PV;PVI:NNN;KPOS:ORD 
TOF Be RO i ieees THE ORDINARY STGN TEST TO CALCULATE THE K 
sraristic P-VALUE, AND CONFIDENCE INTERVAL AS A TEST FOR MEDIANS. 
TAST APTION WILL DISPLAY 42 TABLE OF CONFIDENCE INTERVALS OF 
AE DERED STATISTICS WITH CONFIDENCE COEFFICIENTS. 
SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: BINOM, NORMCDF, NORMPTH, 
INPUT ,AND QUANC. 


' DID YOU _ENTER THIS PROGRAM FOR THE SOLE PURPOSE OF GENERATING CON 
PIDENCE INTERVALS FOR A SPECIFIED SAMPLi SIZE AND QUANTILE? (Y/N). 


WW<+ 
oe Wiis 'Y')/Bu 
THE NULG RYPORASSiOo S2A725 > PEE ape MEDIAN (M) IS EQUAL T 
TRE HYPOTHESIZED MEDIAN (MO); HO: @M MO 


DODDDD 


! 

9 

WHICH ALTERNATIVE DO YOU WISH TO TEST?! 
ae ENTER: 1 FOR H1: M < MO; 2 FOR Hi: M > MO; 3 FOR H1: M = MO.! 
C+ 


+( (C21 4(CH2)40G-—o ys 
B2:'ENTER: 1 £OR SINGLZ-SAMPLZ PROBLEM; 2 FOR PATRED-SAMPLE PROBLEM.'! 


eos NOTA Z 2/2 
+(AA=2)/L9 
A INPUT DATA FOR SINGLE SAMPLE CASE 
X*+INPUT 1 


NNN<o0X 
MO+INPUT 3 
D+<X-MO 
>L1ii 
Qa 
L3: 5 ane Zz 
$13 (os2)4R 
Y¥<(Qt+i)+# 
DD+xX- 7 
NNN<+oDD 
MO<INPUT u 
De (X=-Y )-MO 
a 
511:A<(D=90)/D 
A 
N+paA 


“Ergss54 439) 


PAIRED SAMPLE CASE 


COMPRESS 92 TO REMOVE ZEROS 
RECORD LENCTH OF A AND ASSIGN TO N 
KEEPING TRACK OF POSITIVE SIGNS 


Z 
Pi: eyo oes es 
Eee J 
=> 


a IFN IS GREATER pean 20 USE NORMAL APPROX W/ CONTINUITY CORRECTION 
NORM: Z1+( (KPO (0,5xN) (9. 5x (140.5) 
Ze ( (KPOStO. 53-00. SS TYLI Nx0.5)) 
PV<NORMCDF Z 
PVI+i-(NORMCDF Z1) 
9 IF PAIRED SAMPLE TEST CO TO L17 FOR OUTPUT STATEMENT 
be: Puse2x (1 /(V, coun 


Pyaes 
NS:DVine(3 Jo (P PV ,PVI,PV3) 
i CONPUTAPTONS ARE BASED ON A SAMPLE SIZE OF: ',(#N) 


'THE TOTAL NUMBER OF POSITIVE SIGNS IS: ',(®KPOS) 

+(AA=2)/L17 

'THE P=-VALUE FOR HO: M = ',(3MO),' VERSUS Hi: M ',(sLOGIC(C;1]),' ', (smo) 
»' 2S3', (4oPVM(C31)) 

{ 

+£18 
Li7:' THE P-VALUE FOR COMPARING THE MEDIAN OF THE POPULATION OF ' 

' DIFFERENCES TO THE HYPOTHESIZED MEDIAN, 


HO: M(X-7) = _',(3M0),' VERSUS di: MCX=-7) ' \G@eeOGTCiC: ii.) ' Ferro), *, 2s 
(SBaPYM(C31] ) 


i 
gS Mn YOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (Y/N).! 
€ 

Re 

+QUANT 
A INPUT SIZE OF CONFIDENCE INTERVAL 
L16:CC«INPUT 

ALPHA+(100-CC)+200 

gg eso) Nae s 

COMPUTING CONFIDENCE INTERVALS BY EXACT P-VALUE 

"“CDP+BINOM NNN 
INDEX POSITION OF CDF FOR ALPHA + 2 


"Be+/ (CDESALPEA ) 
+(B>0)/SKIP 
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aS 
S SKIP 
QA COMPUTING CONFIDENCE INTERVALS BY NORMAL APPRO 
NORM1: KALPHA<((0.5x(NNNx0.5 ) )x(NORMPTH IEPRAY +0. SxNNN )-0.5 
ats aOUND SALZEA DOWN TO NEAREST INTEGER AND INCREMENT BY ONE 
pea 
ber Om 1)/L7 IF SINGLE SAMPLE CASE GO TO L7 
RQ ? CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 
D3: ee 


7¥<00RDD 
CT eO2DDEB »ORDDE (Z¥~(B-1))] 

tA 1) (38sec PERCENT CONFIDENCE INTERVAL FOR THE MEDIAN OF THE ! 
i POBULATION OF DIFFERENCES IS: 


Cot (46701) a Mera x-Y) Ss ',(ecrC2]),' )' 


! 
+QUANT 


a a oe ee 
=a 


PRPPPRPRPPHPHPDONWNOWNHWWWWMaWwoO 


f CALCULATE AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 
Zio ORDe 
GArorperg} »ORDX( (2Y-(B-1) )J 
A PEACENT CONFIDENCE INTERVAL FOR THE MEDIAN OF THE POPULATIO 


iv 
1 


{( ',(s@701]),' s MEDIAN s ',(SCIE2)),' )° 


Bu: 'ONTER DESIRED SAMPLE SIZE (SINCLE INTECER VALIE).' 


20a 
ora staco vie UA>100))/E1 
NNN ZUANC @ 404 


\ aw THIS TABLE GIVES CONFIDENCE COEFFICIENTS FOR VARIOUS INTERVALS W 
| ORDER STATISTICS AS THE END POINTS POR THE ',(8QUA),'TH QUANTILE.' 


+0 
Ei: \ERROR: THE QUANTILE VALUE MUST LIE BETWEEN 0 AND 100; TRY ACAIN.! 
>+Bt1 ‘ 
E2: | ERROR: YOU HAVE NOT ENTERED A VALUE OF 1 OR 2; TRY ACAIN.' 

»B2 

BQ: | ERROR: YOU HAVE NOT ENTERED A VALUE OF 1, 2, OR 3; TRY ACAIN.'! 


+53 
Eu:'ERROR: YOU HAVE NOT ENTERED A SINGLE, INTEGER VALUE; TRY AGAIN.' 
+BY 

v 


FLEWWWWOWWWWWWNYND NNIODND DPNFFHNAAEAR PRR OOODODOOOOOCOMO~INMNEWNHNHPOWMO™~ 
POWONMHAMNEOANFROWOs] MMF WP POWOWUNE DNYNHPO NMOWDYNMAMN ECON OA «a s+ at ne 


OTe aaa ee oe ae rwew erwl Ooo wows oho ou} rere ws CO OE Oa 
Hprprp 
aa ae ae eae Ss Sea SS Se (mami ALAN ae ee ee ee ee ee ee | (Jt 


PRPPPPPHP PPPs Ppp sy tp Pepe peepee 


R;CHA; 


Ky SPEARMAN; x; Pee 95 
THIS inept N iohpires THE SPEARMAN R STATISTIC WHICH MEASURES 


a THE DEGREE’ OF CORRESPONDENCE BETWEEN RANKINGS OF TWO SAMPLES. THE P- 
a VALUE IS GIVEN FOR TESTING ONE AND TWO-SIDED HYPOTHESIS OF ASSOC. 
a SUBPROGRAMS CALLED BY THIS FUNCTION INCLUDE: TIES, TIESK, SPEARP, 
a INPUT, SPAPROX, INTERP, AND THE VARIABLE PMATSP. 
R<INBUT 2 
+1+(Q+1)4+R 
Y<(Q+1)+F 


a CA ae TO CALCULATE THE STATISTIC AND ASSOCIATED P-VALUES 
Ae SeMANE 2 
+092 5-0 )/ Li 


a PypI ASC. | 
L1:CHA<'!DIRECT! 
£2:PV<2xA(2] 
>(PYs1)/L3 
L31'SPEARMAN''S R EQUALS ',(4sAC1iJ) 


' THE P-VALUE FOR HO: NO ASSOCIATION EXISTS VERSUS' 
; ; H1: ',C8CHA),' ASSOCIATION EXISTS IS:', (4saAl2]) 


THE P-VALUE FOR THE TWO-SIDED TEST OF HYPOTHESIS IS:', (4sPv) 
V 


MNEWNrPOW OID MF Ob) ' Otsa_sa_sa_ www ws 


ee 8 SS OE — | 


EAT TOT ESI 
NNNYNNP PPP PHAR WOONIONONE OPP 


oa 


QUANT: WOULD YOU LIKE CONFIDENCE INTERVALS FOR A SPECIFIED QUANTILE? (Y/N 


oii! SEs DESIRED QUANTILE; FOR EXAMPLE: ENTER 20, FOR THE 20TH QUANTILE. 


ONE ANE III 
NMNNNNNRPRPRBRRPRPPRBMRPRDON OMFEWPDOH 


FWP OM ONMAMEWDH OW ONOMMEWpoHrOtxr ra tania wias 


SNININNONMNAMNMNMNNMOMONMNNNANANNOMEEER PE EE EF CWWOWWWW WWW Wh 


— sk A Uk An en i A ll 


SAN EWNMRPOWDONMNW FUOINHPOWDODNAM EWN HOW ONMNM FWNHOWONHMMEWHHOWONIMM 


rewy www YT wy eee ee Ieee 


OO OOOO MOO JI~IJ~I~J 


TPLU 
1 


v WILCOX: A:B:D:E3F: PY? ; $24 5 Z:DEN:NUMZ:NUMZ13;PVAL:X:MO;:N; PELUS s CDP: PALPBA; iA 
LPHA;H:CI:¥:Aa: BB 3C DD: PV:POS;fPOS;NM:PVI:7EOSi:NNN;C;2VM:2V3;328 
RAP: 1 hom Ch iw pe pene wre 
THIS’ PROGRAM ise ize CATLCOXON SIGNED RANK TEST TO CALCULATE THE T 
STATISTIC, P- VALUE, AND CONFIDENCE INTERVAL AS A TEST FOR MEDIANS. 
SUBPROGRAMS CALLED’ BY THIS FUNCTION INCLUDE: TIES, WILP, NORMCDF, 
NORMPTH, CONFW, AND INPUT’. 


THE NULL RYPOTHESIS STATES _- THE eee MEDIAN (M) IS EQUAL T 
TRE HYPOTHESIZED MEDIAN (MO); HO: M MO 


DDDD 


! 
! 
Q 
; , WHICH ALTERNATIVE DO YOU WISH# TO TEST?! 
aoe ENTER: 1 FOR Bice < MOs 2 FOR #i: M > MO; 3 FOR Hi: M = MOQ.! 
— 
= (Yce4 ya (c#2)a (G23 ))/E3 
B2:!ENTER: 1 FOR SINGLE-SAMPLE PROBLEM; 2 FOR PAIRED-SAMPLE PROBLEM. ' 
— 
eee eee 
>(AA=2)/ 
A INPUT DATA FOR SINGLE SAMPLE CASE 
XA*+INPUT 1 
NNN<0X 
MO+INEPUT 3 
=X -MO 
>L11 


A PAIRED SAMPLE CASE 
LQ: era 2 


<4 
+14 (Q+1)+R 


Y<(Qei )+2 

ies 

NNN+oD 

MOMINPOT 

De (X-%)-MO 
a COMPRESS D TO REMOVE ZEROS 
211:A+«(D=0)/D 
om A RECORD LENGTH OF A AND ASSIGN TON 

+p 
Q KEEPING TRACK OF POSITIVE SIGNS 

POS+(A>0) 
A ; TAKE THE ABSOLUTE VALUE OF A; ASSICN TO 38 AND ORDER 3 

» 

BeBtAB) 

REQRDER POSITIVE SIGNS TO COINCIDE WITH PROPER POSITIONS IN B 
"POS*POSt ({4)] 
CALL FUNCTION TO BREAK TIES 
“Est gree B 


CALCULATE TPOS BY ADDING ACROSS ALL POSITIVE VALUES OF E 
"rpoges/(BOSxE E) 
LECoieTPOS 


A Ca § SIZk OF LEFT TAIL OF PROBABILITY DISTRIBUTION 
NM#(L ((+/9N 32) +1 
AON iG O TO STATEMENTS BASED ON LENCTH OF VECTOR E 
oe WILE ¥ GENERATE NULL DISTRIBUTION FOR TPLUS 
a IF TPOS FALLS ae Paap HALF OF PROB DIST CALCULATE PVALUE AS NORMAL 
» (PROSE M- -1))/ZR 
A THE phy tee USE TRE NEGATIVE T STATISTIC 
eee eal -TPOS 


LCE OS lS ea ie USE BOTH THE INTEGER ABOVE AND BELOW AS TPLUS 
Precis (GI2POS)= 0)/NON 


+L 7PO 
>PV+1-( (FC 7TPLUS)+FP(2PLUS+1) )4+2) 
PYT« (PCTPLUS+1]+¥(TPLUS+2] )+2 


+CHEC 
Dea 


+P2 
GO:PV+1- meee a 
P2:PVI«F((T 

Chick: >(7POS13 (Hi =1))/L6 


Pee 
PVLery2 


+G6 
9 COMPUTE NORMAL APPEOX. W/CONTINUITY CORRECTION FACTOR 
Ea :TRAPH(Nx (N+1 ) )£u 
HOWE (TPOS +0. 5) STRAP 
Be Le TPOS- leet 
DEN ( (Wx (W013 ( 2xN} 41) 424) 40. 5 
Z<NUMZ+tD 
Zi<NUMZISDEW T APPROXIMATION WITH CONTI x ECTION F 
A NOITY CORREC. AGTOR 
vile (B08- TRAP) 
DENT< oar yee oa Eon Cae 0.5)*2)4+(N-1) sao 
DENT1<(( (Nx (DEN*%2 ) )4#(N-1)) (NUM+0.5 )*2)#(N-1 )*0.5 
TC+(NUM-0.5)+DENT 
TCLs(NUM+O. SST BENTs 
A COMPUTE AVERAGE OF 1G AND ZC 
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[ Bee ee 7s Rea ND 
eerste seme De 21 ))+(1~( (N=1) FDISTN TC1)))+2 
aa (NORMCDF 2)+((N-1) TDISITN TC))+2 ,' 
SECOND: PV<((1~(NORMCDF as we TDIGeu82 C2 ))+2 
PVI<( (NORMCDP ag ITDISTN 7C)))+2 


Eyes 
NS; ByH«(3, 1)o(PVI PV. Ea) 
‘COMPUT TIONS ARE 3ASZD°ON A SAMPLE SIZE OF: ', (aN) 


“TRE TOTAL SUM OF POSITIVE RANKS IS: ',(8TPOS1) 
IF PAIRED SAMPLE TEST CO TO £17 FOR OUTPUT STATEMENT 


beter tAbAR RPP OWNDDOWWWMWNwWOo 
rPOOOO9O BDOOOV CWO 1WMMNFWNHHOWC 


OwWO NHN CFWOHOLw asain a 


>(AA=2)/217 
'THE P-VALUE FOR HO: M = ',(8MO),'! VERSUS Hi: M ',(s8LOGIC(C;1]),' ', (MO 
5) £S:', (OSPYMCLC;1)) 
1 t 
£13 
Pee oe VLU FOR COMPARIVG aE MEDPAN OF THE POPULATION OF ' 
DIFFERENCES LO THE EYPOTHESIZED MEDIAN, ' 
SOs O(k=-7) Coy. Vemma fils A-£) ';(SL0GlCCC3:23),' ',(3M0O),', £ 
:', (USPUMEC; 3133 


et YOU LIKE A CONFIDENCE INTERVAL FOR THE MEDIAN? (Y/N).'! 
<< 
BB='Y¥')/Z16 


oa 


tobe } 
tro 
fo 
hy) 
mC y 
4 + 


(209-c2)42 9 


CvNN> oe ag 8. ” NORMAL APPROX. 27OR CONF INT OF SARGE SAMPLE SIZE 
aes 15 w= 
COMPUTING CONFIDENCE INTERVALS 3Y EXACT P-VALUE 


“CDP+HILP Sai 
INDEX POSIZION OF CDF FOR ALPEA + 2 
“TALPHAs (+/ (CDPsALPHA)) 
+(TALPAA>0 )/ JUMP 
TALPHA<1 
~JUM? 


SSS eae Oe eee el 


Dp ti 


Spd pd pr pdp apap ag dy rapa prpapa 
NHB EAEAP EAR Rp 


ah 4 AAA gM eg ay ee ee ee a a ee | 
mM 


3 COMPUTING CONFIDENCE aoe ao a NORMA APPROX. W/C.F. 
CPAEEESC ASIA SOEs ABB aU + 
TALPHA+(( (DENOMZx(NORMPTH ALPHA) + (NNN (NNW? 4) )- 2)74) 
g TALPHA DOWN TO INTEGER VALUE AND INCREMENT BY ONE 
" PALPHA+LTALPHA 


A ** TP ONE SAMPLE CASE GO TO L7 
JUMP: +(AA=1)/L7 
A CALCULATE AND PRINT OUT CONF. INT. FOR PAIRED SAMPLE CASE 
LS:CI«ZALPRA CONEW DD 
UA toc)! PERCENT CONFIDENCE INTERVAL FOR THE MEDIAN OF THE * 
' POPULATION OF DIFFERENCES IS:' 


' Ge oec ri 1)).' = Mepmiy«x-}) Ss ',(ecr(2)),' }* 


rot 
+0 
CALCULATE AND PRINT OUT CONF. INT. FOR ONE SAMPLE CASE 


Q 
L7:CI«ZALPHA CONFW X 
Wala » (@ 8CC),' PERCENT CONFIDENCE INTERVAL FOR THE MEDIAN OF THE POPULATIO 


' ( ',0@cI(C1]),' s MEDIAN s ',(€ecI(2]),' ):! 


+0 
E2: \ERROR: ZOU HAVE NOT ENTERED A VALUE OF 1 OR 2; TRY AGAIN.' 


ria 

EQ: \ERROR: YOU HAVE NOT ENTERED A VALUE OF 1, 2, OR 3; TRY ACAIN.' 
+B3 

V 


FUNPOWONMUW FWHNROWDNMTMEWNHHROW OWN FWNHOWOYNOMWNEWNPE 


CLL LL LDN ORR Oe eee eee ew ww wow we 
PAPDPAPAPAPAPAPAPAPA PAPAPAPAPAPAPAPDPaPaPapapapapapapag ap rpay 
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APPENDIX F 


LISTINGS OF SUBPROGRAMS BASIC TO BOTH WORKSPACES 


VY CBN+«BINOM N X3;K:CDP 
a THIS FUNCTION. iss A’ SUBPROGRAM OF THE SIGN TEST (SIGN). IT CALCULATES 
a tHE CDF OF THE SrNOHAL WHEN PROBABILITY = .5. N=SAMPLE SIZE. 


inf Wh, 


K<0,1N 
CBN NCE IN) ee 


, V U<+ZL CBIN R 

Z THIS FUNCTION IS A Seams OF CONFIDENCE INTERVAL GENERATOR FOR 
a ae Bone QUANTILE. If RETURNS THE VALUE OF THE BINOMIAL CDF AT R, 

u N,P=L WHERE N= SAMPE E SIZE AND P=PROBABILITY. 

ae 
On 


Fe De1i4001] DEDC1] x (LC2] ~% te/e1t 4201] x C1-002] )*« L501] -~14+114+001] ) CR+ 


Cc.aA_oow 


Vo CON*XX Bo Vis SB AGS «AMA en Sea 
PHILS 2£UN 


1) A CTION IS A Shipbocnan OF WONPARAMETRIC LINEAR REGRESSION 
a5 4 (YE PIR). If rALCo ILA ES. THE THO-POINT SLOPE POR ZACH PAIR OF POINTS 
Bn (Xe AND AND Xz=X2. EL. OF OBES ARE ORDERED AND 
u] 9 USED ?O FIND THE “Sa pTBENEED INTEROAL POR 8, THE SL Ove OF ESTIMATED 
3 a EQUATION. XX= X DATA SAMPLE AND Yi= Y BAPA SAMPLE 
7) 9a RECORD THE SIZE OF XX AND INITIALIZE VARIABLES 
g BB+pXX 
10] Raso- ° 
13 q THIS LO LOOP COMPRESSES XX AND YY DOWN TO WHERE THE XX < ALL OTHER XX'S 
<J4 
13 AS ECAAT xx) 
44} XR«A/XX 
15 TRAY 
16 Seen 
no me =0)/L3 
ao 
19 g THIS L LOOP CALCULATES THE SLOPE OF EACH PAIR OF PAIRED DATA. 
24 $2(71 CAA - ~YR(C] )+(XX(AA]-XR(C]) 
22] Sg§<SS,5 
393) +(e BY/L4 
2u) £23: Cea its B)/L 
25 gone (oss), (33581) 


V CONFM+AA Sa eeate BB;A $ liee ee 
PHT Ss orUN 


cet 
AaINNER AND QUTER LOOPS CALCULATE ALL POSSIBLE DIFFERENCES; EVERY Y 
a ELEMENT MINUS EVERY X ELEMENT. H VECTOR STORES THESE DIFFERENCES. 


L2:G<i 
L1:E+DCF]-CCC] 
H<H,E 


1] a ny is’ 5 ‘CiBb ROCRA Cie OF THE MANN-WHITNEY TEST (MANW). 
25 a If COMPUTES CONFIDENCE INTERVAL ENPOINTS FOR THETA, THE SHIFT IN 
33 «a LOCATION BETWEEN X AND Y. AA=INDEX POSITION OF C.I. ENPOINT. 

4) a BB=COMBINED DATA SAMPLES. 

5 Rs pee SIZE OF X VECTOR TO A; INDEX POSITION FOR CONF INT TO B. 
wie << 

a rad C3] 

3) ASSIGN X VECTOR 70 C: Z VECTOR 70 D 

3) c+a+aB 

10 D<A+BB 

113 A REORDER X AND Y VECTOR VALUES TO ASCENDING ORDER 

13) 55088) 

3 D<+D 

su) a INITIALIZE VECTOR H AND VARIABLE FP 

5 B<p0 

16 

17 

18 

19 

e549 

[21 


Sa ae eee ae aS ee ae 
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Crna ss eT Ow OL oe 


NNNNMDH MH HMRPR MMH PW O UMN EWP 
FWNRPOWWNOMF WHOM Ola.4 sas 


LL Nae a a a a a aoe ae a aoe dt 


FEE FFE EE EWWWWWIWIWWNNYNNYNNNNNHMMH HR PRR wWO~WUMU FEW} 
OnNMHOMNECWNOHrPOW OWN FIN HOW ONMMNEWNHrPOWO~ DWMM F WNHH Ot 4444 _2_a 


a ee ee eS ee eee ee ee ee ee ee eee eee ae ee a ee ee ee ee eee 


=(és9c))/L4 


ee ee 
REORDER 8 VECTOR VALUES TO ASCENDING ORDER 
“aed CAa] 


INDEX APPROPRIATE _CONF INT VALUES FROM # 
"GONFH+ELE3 pa (oe (a1) i 


VY CONF<A CONFW 3:C:D:E:G:H8 

THIS FUNCTION "IS A ibepboeRAe OF THE WILCOXON SIGNED RANK TEST 
(WISIC). IT PROVIDES CONFIDENCE INTERVAL END POINTS BASED ON THE 
AVERAGES OF ALL PAIRS OF DATA SUCH THAT ALL XI $s XK. A= INDEX 
POSITION OF C. I. END POINT AND B= INPUT DATA. 


START H ACCUMULATION VECTOR OFF WITH ORGINAL VECTOR VALUES 


D DDDD 
a 


a<20 
Bate c 
De 


"Qu2srDE LOOP INCREMENTS D AND RESETS FE TO D 
°D<D 


a 

“elo 

Q rabeeee LOOP CENERATES NEXT SET OF AVERAGES AND CONCATENATES TO ORIGINAL 
b1 


Epp hcune 


(Bele) re nTsNUE INNER LOOP UNTIL EF EQUALS THE SIZE OF C 
"CON NeSyoE OUTER LOOP UNTIL D EQUALS THE SIZE OF C LESS ONE 


= Cc 
yeh eC) -2 0/5 
ORDER FINAL ACCUMULATED VECTOR HB 


"+H CAR) 
Q INDEX el INT VALUES OUT OF @ 
CONESELAL, ,2C Gee )-(A-1) 


a P«DP FPDISTN X:A:M:N:RM:RN:LM:LN:SM:SN:M2:N2 

THIS A ONCTTON Mi SUBPROGRAM OF THE ipo sxac- WALLIS TEST (KRWL) AND 
a THE STUDENT T DIST. (TD DISTN). IT CALCULATES APPROX. CUMULATIVE PROBS. 
g OF X USING THE F DISTRIBUTION W/ DF=(M,N) DEGREES OF FREEDOM. 


A TAKES VECTOR ARG. 

At ,MxX+N+XxM<«1i+DF ,N«1+DF 

+Lx1(M>2)vN>2 

a ae ae 87. FOUR SPECIAL CASES. . 
+2 


ore a a neo" 5) +042 
oo P+Ax0.5 a= 
24: P<1-(1-4)*0.5 
£22: P<A 
+0) 


a BEGIN THE GENERAL CASES. INITIALIZE THE QUANTITIES. 
L:RM<«(1|M+2),0/RN 2 
“Site (ah Ox <BMSL Me M20 5xM- 2 
xBNSLN«|N2<0.5xN- 
a TREAT S27 ASE N-EVEN PrRSe, 


SWeies/#\(1=A)e. x14(7~140.5xM)41N2 
+>MEVEN 


Vena ees 
Ne -2 
A TREAT THE PORTION OF SN THAT DOESN'T CHANGE W/ M ODD OR EVEN. 
>(2+1+0LC )x10=N 
SN+SNx1++/x\(1- Bayo X142XM2414+2x12IN 
A fhe. 22E M-2VEN SUBCASE 
~( (M22), 2=2{M)/MEVEN ,34+170LC 
SN<+SNxx/13E2%1M2 
+>MEVEN 
a TREAT THE M-ODD SUBCASE 
SN<SNX2<401 
>+ENDX11=M 
SN<SNxXx/144 142x10.54M2 
MEVEN:*MODDx11=2|M 
+ENDxX12=M 
SEL at 7x\AS. X1-+2x1M2 
+>END 


QR NEXT TREAT THE SPECIAL CASES ve ODD-éM. 
MODD: SM< (OxM=1 )+((Ax%x0.5 )#042)x3=M 
>ENDx1(1=M)V3=M 
SM« (A405 )x(t0t2 )x1++/x\Ae0.x1-t14+2x10M 
Bee tS CONe ANY 2) + (4xRNxRN- RM) +( (RMXBNX 10A*0.5 )404#8 )-2xRNxSMx(1-A)*0.5 
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SN TTT TET NTE 
FEWWWWWWWWWWN NNNNMNVNNYNPEPRPPRPPPBPBREOONODMNEWNH 
FPOWDNMAMNEWNHROW DNHNMNEWNRPOWONMNMN FWNH Otsasavae wn a 


ee 


OL EE OOo 


FE FWWWWWWWWWWINVNNNNNYNNPRPPPRPPPRPPROONOONEWNER 


V XX«N INDEXPLS 3: faa 3B; ci DD; axhi NN; NC 3H 
THIS FUNCTION red GRAM o peiskShe tiaitrs (KRWL), AND MANN- 
WHITNEY (MANW). one oe Eee MATR orx OF RANKS OF THE ORIGINAL DATA 
WHERE THE Ppoh Row 6% THE MATRIX IS THE FIRST SAMPLE. IF TIES OCCUR 
IN THE RANKING WITHIN A SAMPLE, SUCCEEDING TIZS ARE INCREMENTED 3 
OWE. N= SIZES OF THE SAMPLZS BEING PASSED IN B. B= ALL DATA SAMPLES 


RECORD SIZE OF N TO DETERMINE NUMBER OF SAMPLES IN 8B 
LOCATE LAaGEST SAMPLE SIZE 
BB<[f/N 
DD«BCAB] ORDER 8 SMALLZ@ST TO LARGEST 
Pa 
XX«(AA, Pane UP MATRIX OF SIZE REQUIRED TO STORE SAMPLE RANKING. 
< 
OO Nea TENATE 0 AND N 


:N 
"FIND CUMULATIVE SUMS OF SAMPLE SIZES. 
eas 


<0 
a THIS | Aaa INDEXES OUT ORGINAL SAMPLE VALUES FOR FURTHER CALCULATIONS. 


AA+pN 


Se: DD D D PDYOYPNYD 
= 


9 atincleetsintect } ORIGINAL SAMPLE VALUES IN 8 
xX1+¢ ze 

ve po Pa CORD POSITIONS OF ELEMENTS OF X VECTOR IN 8B AND ASSIGN TO C 
CULL 


De 
R THIS LOOP DETECTS TIED Ee Cee AND INCREMENTS THE INDEXING OF 


EAGA See TIEDPEOST ZEON 
A TEST DTH ELEMENT OF C AGAINST BEST OF C FOR TIES 
8 


eae F SQUAL TO THE APPROPRIATE SIZE (DTH SIZE) ZERO VECTOR 

“U0 

Be F AND VECTOR OF 0'S AND 1'S (1'S APPEAR WHEN TIE OCCURED) 
© % 


A ADD RESULTANT F VECTOR TO C 
CHC+F 
ee 

CONTINUE POR ENTIRE C VECTOR 
"= (Ds (pe C))/ 52 


XX(CC31<+C, (BB-NICC] )po 
=(cechA ri Cee 


V IN*INPUT A -Y 
a THIS VENCTTON bbks MOST OF THE INPUT PROMPTING AND ERROR CHECKING. 
a A=THE TYPE OF PROMPT DESIRED. 

a IT ts A SUBPROGRAM OF SICN, WISIG, MANW, KEN, SPMAN, AND NPLR. 
+Cid G2 62, bh .L5 bad 

bai tENr ER THE bara (MORE THAN TWO OBSERVATIONS ARE REQUIRED). '! 


4h Se IN)s =0)/E1 


N)=2)/E1 
oe 'ENTER X DATA (MORE THAN TWO OBSERVATIONS ARE REQUIRED). 
“hoeinegyee 
yen Nowe R Y DATA (NUMBER OF Y ENTRIES MUST EQUAL NUMBER OF X ENTRIES). 
nee 
¥)=0)/E2 
{fp OR z(oL) 3 7E2 
Ide (pk, ¥ 67 


L3:' ENTER THE HYPOTHESIZED MEDIAN.' 

ee 

eee THE HYPOTHESIZED MEDIAN FOR THE DIFFERENCES OF THE PAIRED DATA.' 

seth eyes ee 

53: ENTER THE DESTRED CONFIDENGS COl2 feo: 

1 £9 FOR EXAMPLE: ENTER 95, FOR A 95 PERCENT CONFIDENCE INTERVAL. ' 
= ERM EB}eSIIR00 =" 


PN IJ20 755 
Ei: ERROR: THE SIZE OF YOUR SAMPLE IS LESS THAN THREE; TRY AGAIN.'! 
=(A= =1)/L1 


ze : aa SAMPLE SIZES ARE NOT EQUAL; WANT TO TRY AGAIN? (Y/N).! 


woh! )/ 52 
ENTER RICHT ARROW + TO QUIT.' 


100 


tf A 


CT Er EEE ee eee 


NN NNMNNYDY HAM RHR RPHRHMHWOO QYMDU FE Wp 
SONNE WPOrH OW ONDINE WOH OL «a t_a_a_ naan 


(M44 B_A nd A ek 


VP OWOYQVD FWRPOKHPOW WO INDMEWPVORPOWODNIDAMNEWRWOHO Lassa el 


Le AAA 


LF FWWWWWWWWWWIIhhVWNDNNNNPHPHPHHPPHPPHHE DONMOME Wh 


BO 

ER: \ERROR: THE AIPOTHESIZED MEDIAN MUST BE A SINGLE VALUE; TRY AGAIN.! 
*+(A=3)/Z3 

>D4 
Ev: ERROR: THIS VALUE MUST LIE BETWEEN O AND 100; Z8Y AGAIN.! 


25 
ES: ' ERROR: TaIS VALUE MUST BE AN INTEGER; TRY AGAIN.' 


>Z5 
v 


Y INTER+A INTERP B:C:D:E:F er 
THIS FUNCTION IS &A sob eRoGhan OF THE KRUSKAL-WALLIS (KEWL) , 

KENDALL'S B (KEN), AND SPEARMAN'S R (SPMAN AND SPMAN1). WHEN 

C AS THE LEFT ARGUMENT IT CALZULATES THE 

A TATED P-VALUE BY INTERPOLATION OF THE TABLE VALUES OF STATS. 

AND CDF'S WEICZ ARE PASSED AS THE RIGHT SBGUMEN?. 


op oEPARATE THE CDF TABLE AND STATS. INTO SINGLE VARIABLES. 
Ge 2 °3i3 


bppypp 
daly 
Ts 
—A 
—Y 
by 
Le) 
far 
xy 
Ge, 
hy 
by 
—| 
Hy 
Sn 
ts 
Hy 
Hy 
ts 


s° 


PIND WHERE A FIRST EXCEEDS OR EQUALS ONE OF THE TABLE VALUES. 


R 
coc (ABP 
Bde eee SEA7- DOES NOT EQUAL ANY TABLE VALUES SET INTER = “1. 
-> 
8 _o._ INDEX LOCATION OF FIRST OCCURENCE OF MATCH 
01 
A TF INDEXED POSITION EQUALS ONE INDEX P-VALUE OUT OF GG. 


>(D>1)/51 
eee Cowen 


A OTHERWISE CONDUCT INTERPOLATION TO GET PROPER P-VALUE. 
LIVEePP(D-15-FFCD) 


F(D-1)]-A 
G+GG{D1-GCcCD-1) 
PxG)t2 
INTER+GOCD-13+PL 
>0 
L2: :INTER<"1 


ag KENP+RENDALP N:A:83C:D:E:NN:X3XX:AA3 F3T:G;:P 

THIS FUNCTION IS A SuBbR OcRAN OF’ KEMbaiT B (KEN) AND NON- 
A PARAMETRIC LINEAR REGRESSION (NPI? ). IT PAECULATES THE CUMULATIVE a 
8 DISTRIBUTION OF B FOR A st ON ee Fle Ny. 


xX* 11 
“ye coi FREQUENCIES FOR X FOR SAMPLE SIZE NN. 
“(LCN pee eS re OF RIGHT PROB. TAIL VECTOR 
A COOTER LOOP INCREMENTS THROUGH THE N SAMPLE SIZES TILL THE DESIRED ONE IS 
GENERATED 
£1:D<p0 
C+0 
F+o s 
NNeN 
Be ( (NX CNW 1))#2)+1 
g INNER LOOP GENERATES NN+1 FREQUENCIES FROM THE VECIOR OF NN FREQUENCIES 
C+C+XC(B-A)] 
Mash 
&e 
F D 2954 eeu NN AND INDEXES OF X STILL REMAIN GO ZO Ly 
“Sago RNA (CB-A 
"acca -A)SPITED sez" Pont NTE TO INCREMENT THRU OLD FREQS 
ne WHEN LEFT HALF OF NN+1 VECTOR IS COMPLETE GO TO £5 


a als BOOP aLiLOWS ONLI NN TERMS TO 3E USED TO GENERATE NEW VECTOR 


- * 
3 ie has 


CectetiB- ~A)]-X(AA]J 
D+D 


C 
unCCE- ~1)#2)+3 )> (B= A))/LY 
a NVERT VECTOR D AND ASSIGN TO E 
L5:£+«eD 


AR FOR NN OF APPROPRIATE SIZE EITRER DROP FIRST VALUE OFF E OR NOT 
» (NN=3, oom eto, tia, 24, 159,18,19 3/703 


E 
Pee NEW VECTOR OF FREQS X BY CONCATENATING D WITH E 
Pe CWN EN) /L1 COMET OE ONTIL SIZE OF SAMPLE N IS REACHED 
A GENERATE VECTOR OF CORRESPONDING P STATS OF PROPER SIZE G 
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CO TN eT Eo roe Io 


WWNNYNVWNNDYNNPPBRERPHPHPHEPWOOYUOMEWNHE 
POWONIOAMA EWN OM ODO YOM EWP Ola» 4a Aa 


Ciel 


orm 


PRE RODNIOUW EWE 


WOR Ot A841) 


m 
ay 
on 
4 


PRP RPPPPRPPHPWOOADin KFWNR 
WO IM FWHM OWUsA-aAs 4 4 


Perv ww YP Vw ooo ek lc lem 


af Sle 1(G=-1) 


CALCULATE B STATS | PROM THE P VECTOR 
“rel (axe 3 (Nx (N-1)))- 
TARE ONLY CG ENTRIES PROM THE FREQUENCY VECTOR 
EBLe Ts Wx W=1) 5 +2 


ae a 
B pees TO CDF VALUES AND OUTPUT B STATS W/APPRO. CDF VALUES 
"genps(3, Ne pr CC4\KX% )4#C1N)) TPL 


V MAN«N MANWP 
THIS Pune ton pees i CiebRoe RAM Ob Oobe MANN-WHITNEY TEST CHAN). 


A 
A IT GENERATES THE CUMULATIVE DIST. FOR TREOU STATISTIC ey gee 
A LARGER SAMPLE; M = SIZE OF OTHER SAMPLE. 


ACOMPUTE ag OF TERMS TO BE INCLUDED IN LEFT TAIL DISTRIBUTION LESS 1 


pean 2 
350 > VECZIOR SQUAL TO 1 CONCATENATED WITE MM Zanes 


" P+1 ,HMp0 
9 PDP EQUAL TO THE MINIMUM OF N+M OR MM 
Pel /C(N+M) MM) 
AG LM MM SET Q EQUAL TO THE MINIMUM OF M OR MM 
~ 7 
A GO TO ee DENOM IF MM IS LESS THAN N+1i(SIZE OF X+t) 
>+DENOMx1 (MM<N 
Boy Eo eels CENERATE FIRST BLOCK OF RECURSIVE RESULTS USING NUM LOOP 
L2:9<MM+1 


Rta eee eer USED IN GENERATION OF FI&ST SLOCK OF BECURSIVE RESULTS 
Bleenr ey 

A ere NEW DECREMENTE ED VALUE TO U AND TESTS 228" % TEIS Were 

>I1~x Le ee 


+E Pel einer 
a GENERATE FINAL RECURSIVE RESULTS USING DENOM LOOP 
Bengusses 


D 


a PEP EOE ot Oe USED IN GENERATION OF FINAL RECURSIVE RESULTS 


dw oe 


-$] 
Bee oe 


>DYUXy See ) 
VERT FREQUENCY TABLE TO CDF VALUES FOR FINAL OUTPUT 
"wane 4\ 238 CHE Gait) 


SS 
ts 
ty 
C 
Ss 
h 
Ibe 
Sh 
= 
oO 
pee 


V Z<NORMPTH P: 
a IMPLEMENTS AL 


a FOR A VECTO NDING NORMA 
FOR GREE: 


79134 4.85014127135 2.32121276858 
«(]@0.5-]S)*0.5)e¢.% 1 2)+.xD) 


ee NON So ONO” 
*¥ FR Pe 


"OWD war Ow il w 
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OtsA_A_#_A_A_w of A 


Hi O~IMU Flat e 
| 


Mei rrabPrrearPPereOmr-IgaMneWNH, 
WOM FWMHOU«A A_«A Mew 


mt at 


rw eo owe eo 


PRE WOYNHMNE WOH 


NONNNPRERPRPRRRER EA a OonOMFeWDrH 


ewer ea eee ee eee 


ete) els 
eee OR QF MORE P VALUES ARE QUT QF RANCE.' 


UV P<PERM N:X3:Y32 

a TEIS BUNCHTON on CALLED ee SPEARP. IT GENERATES ALL POSSIBLE 
R eee oN OZ N RANKS. N=SAMPLE SIZE. 5 

~OX Tee = Of 

Z*+PERM N-1 

P<1X<0 
LiseOxyrN<X<X+1 

Ye (wCiN eX)\Z 


¥C:XJ)<«N 
P<((Xx!N-1),N)p(,P),,27 
+51 


V M<N QUANC 
A THIS ronertow'¢ Gives" a CHOICE OF NONPARAMETRIC SON EE a INTERVALS 
a FOR THE QTF TILE OF “A CONTINUOUS POPULATION. =SAMPLE SIZE, 


aA Q=QUANTILE. re CALES THE SUBPROCRAM CBIN. 
Q Cente eS pace? AD The OCfDERS@ogAlloLic NEAREST ESTIMATE. 


19.57 
Q euacit 3 S2ANDERD DEVIATIONS (APPZOX.) 

JL 9.541.35x(2x1-G#1900 )x0.5 

(wi ty X 2 a 
a RUN OUT X 20TH WAYS. 

Ke 1411+. /(T(2) W-2(1) ) 
a CONFIDENCE COEFFICIENTS 

L~-/8(2.90X)0(N,C.01xQ) CBIN( (701) +K-1),(2T[2]-X)) 
M=!O2DE? STATISTICS | COEFFICIENTS! 
aA PART OF FORMATTING OUTPUT 

Mew, {0.23 Oo pee ae . 

Wye? 6 sO QC? bh oCle (CL 623 cK), 1+K-1)) 

re ac Ge t 

Hen, te] Ue. Tams OC Cok), 1 )oL 


V She SPAPROX X 

THIS FUNCTION ig A SUBPROGRAM OF SPEARMAN'S ‘: (SPMANP ) 
IT APPROXIMATES TRE CUMMULATIVE PROB FOR R WHEN PASSED THE SAMPLE 
SIZE IN THE LEFT ARGUMENT AND THE ABSOLUTE VALUE OF R IN THE RICHT 
ARGUMENT. SUBPROCRAMS OF THIS FUNCTION INCLUDE: TDISIN 


penEGOEare gue CONTINUITY CORRECTION 
<«6+Nx 
a TRANSFORM TERE STATISTIC R INTO ONE THAT CAN BE USED WITH THE STUDENT 
Q xe Dist aoay N 

+«(X x((N=-2 )41-(X-Y )*2)*0.5 
Q MOREL THe TDIsST FUNCTION ro CALCULATE THE P-VALUE 
SPN«1 - (N-2) ZDISTN X 


db DDDD 


a SPEAR+SPEARP N3C1:34; °N3:D1;D2:LIM3:R:CDF 
THIS FUNCTION IS'A 24S Cie ReR Th OF SPEARMAN'S R (SPMAND ) . Iv 

A CALCULATES THE EX2CT CUMULATIVE DIST. FOR R FOR THE SAMPLE SIZE 

a PASSED AS THE RIGHT ARGUMENT. BECAUSE OF PEE Lance COMPUTER MEMORY 

a REQUIREMENTS, N IS LIMITED TO SIX ON THE PC AND 7 ON THE MAINFRAME. 

a su PROGRAMS CALLED BY THIS FUNCTION INCLUDE: PERN 


A+Be 
A TWITIALIZE VARIABLES, VECTORS, AND MATRICES. 


C1<)1<+10 

Me(,V)00 

Cae. * (N-1) 
a TETS BOOP GENERATES AN NxN WERAY OF THE POSSIBLE VALUES OF DIFFERENCES 
A BETWEEN ANY TWO PAIRED RANKS BETWEEN SAMPLES. 


L1:8B<«8+1 
MCB;J+C-A 
A<At+1 
>(B 


<N)/L1 
a NOW CALCULATE THE SQUARES OF ALL POSSIBLE DIFFERENCES. 
a aay ie TO LIST ALL POSSIBLE PERMUTATIONS OF N NUMBERS 
A CALCULATE pace ea OF FINAL VECTOR OF & STATS., 


LIM+<1+(((Nx3 )-N 412 
A INITIALIZE VALUES BEFORE INDEXING OUT COMBINATIONS OF ALL POSSIBLE 
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CT 


FEF FE EWWWWWWWWWWHPPNYAA 
FONRPOWDDONOMNEWNHOrROWOYNHMWEL 


Cat tl 


ll lh 


NO FEF ELEE EEE EWWWWWWWWWWNHONYNNYNNNNNNPRPRPRPRPPRHPPRWOYUOU FWP 
WNP O0ODIDMWEWNPOWMODNAMNEWNHKRPOWDONAMWEWHHPROWDONOMW FE WthoHOlLAs saw aw a 


LAM 


Lae 


CO LT NTA Oo eo 


aso ye VALUES. 


E<(N 
rf, LOOP CALCULATES ALL POSSIBLE COMBINATIONS OF THE SQUARED VALUES. 
ae sJ<MCA:DC3A) 
> )/L2 


A<N 
3 a ae ALL ROWS FOR EACH COLUMN TO SUM UP SQUARES COMBINATIONS. 


9 SDD UP NUMBER OF DUPLICATED SUMS OF D-SQUARED VALUES AND COMP3ZE 

a YVECZOR DOWN TO UNIQUE VALUES. : oe 

L3:C1+Cl, $4 (D2) 

Di«D1 i (L/D2)=D2)) 

p2«(((/D2 )2D2)/D2 

+((9C1)<2IM)/1L3 
TRANSFORM SUM OF «SQUARES VALUES TO SPEARMAN'S R STATISTIC. 
<1-(6xC1)2(Nx CCN 


R C1) 
A bec titpas a VALUES ASSOCIATED WITH THE & STATISTIC. 

re (+ moat 

*OZ2M THO 20 Me RIX FOR OUT2UT OF & STATS ANotGo, VYALers- 
sae 2 » (oC1l) eta ene) 


V SPM<«X Seg ae C3;D3DD;013;02:N;DENOMR AX3 IL NS NOMA SE PV ee ae nO: 
UsV *ARBO 3 Od ee 2 
Dols FoNcr ton IS A SUBPROGRAM OF NONPAR LINEAR RECRESSTON (NPLR) 


A 
a AND Boe 2 (SPMAN). IT COMPUTES THE SPEARMAN #8 weal t STI C 
9 AND ASSOCIaTED P-VALUES. PRE LerT ARGUMENT PHAT IS 2ASSED iS THE X 
A SAMPLE: TRE 2&CGAT AZCUMEN™ is. THR] SAN eGs. 
a’ SUBPPOGRAMS OF THIS FUNCTION INCLUDE: ?TZS, TIzZSK, SPEARP, SPAPEOX, 
A INTERP, AND THE VARIASLE PMATS?. 
ps) ORDER Y IN INCREASING ORDER OF X 

V+YCAX] 
Q : ORDER X IN INCREASING ORDER 

0«+X (AX) 
ie AAT COMPUTE CURRENT RANKING OF Y 

<< 

NOW ORDER Y RANKS IN INCREASING ORDER 

“DIV CAV ] 


IF TIES EXIST IN EITHER X OR Y RANKED VECTOR USE MID-RANK METHOD 
" pp«1 TIES Di 
XX<+1 TIES UV 
a Sea ORIGINAL RANKING OF Y WITH TIES RESOLVED 
a 
A RECORD SIZE OF INPUT VECTOR 


N«oX 
A nCAECULATE DIFFERENCES BETWEEN RANKS OF X AND Y VECTORS 
a DETERHINE THE SUN OF SQUARES OF THE DIFFERENCES 
ee THE NUMBER OF TIES IN EACH VECTOR USING THE TIESK FUNCTION 


12 
12 
STIC INCLUDING THE CORRECTION FOR TIES 
ae SV)) 
5)xC(NS-(12xSV))*0.5) 


Anes 
-> 

nee SPEARP TO CALCULATE THE RICHT TAIL OF THE CDF OF B 
"PeSPEARP i 


L42>(w>10)/L9 
PePMATSPE(N 
a GRANGE speP ‘5B? P 70 AN MxN MATRIX 


i] 
cab INTERP TO CALCULATE P-VALUE BY INTERPOLATION 
Fo: PYAL<+ARPHO INTERP 
~(PYVAZL= 13/L7 
2VAL<O. 5 


Sone P_ VALUE USING STUDENT T APPROX. 
3:PVAL<N Seat ARHO 
é: SoPMe (REO) 


iB 
AAG 


a 


104 


Veeck TOTSTN XxX: 


sOxr1sSa/V 


Bb lmVI/1p4140.5-(~)/P 


a a THIS FUNCTION is A SUBP2OGRAM OF TRE CUMULATIVE PROBABILITY GENERATOR 
i R FOR SPEARMAN'S R (SPEARP). If CALCULATES THE CDF AT X USING 

3 a TEE STUDENT'S T pre? ae K DEGREES OF FREEDOM. 

: a 2aes EoNCTION CALLS ON THE 'F!'! DISTRIBUTION FUNCTION (FPDISTN). 

6 P+0.5x(1,K) PDISIN Xx2 

7 >8x10=Vv . 

3 PCLV/19X1+«0.54+V/P 

r10 


tu 


V TI«+3B TIES B -I:N32 °>L3PP:N ; 
TRIS Function ee: suebadghan OF kenbar's B REN SPEARMAN'S 
GSPMARP). ol ia IS (KRNL), MANN-WHITNEY (MAN 5 


Peat 
CE2CAING FOR TizS Jy INCREMENTING TERU TRE VECTOR 
23: C#(%v9 =a (7) 
eh COUNT NOMBER OF DIES; DeenNo TIES GG To L2 


C A 
ae a 2 
lz gn AND WILCOXON mie, 25t G). cueees TAR 2IGAL ARC. VECTOR FOR TIES AND 
C Rg CHANGES THE TIED pasirrons OF THE LEFT ArG. BY THE MIDRANK METHOD. 
f {fees 
YF NO VECZTOR OF BSANKS IS PASSED: GENERATE ONE 
"+ ((098B)= 05 7L6 
BB+iN 
bé: T+90 

L+No0 

C 


e/C 
sie 0)/L2 
A RECORD WHERE TIZS STARTED AND HOW MANY RANKS INVOLVED 


a 
fos 7 gt hGRE EMENT NEXT I 3Y THE NUMBER OF TIES ENCOUNTERED PLUS 1 
ale 4 ges THAN SIZE OF ORIGINAL VecTOR GO TO £3 AND START AGAIN AT NEW T 
Loe 


" PI+BB 

R 
+(Y=0)/0 
Z<0 


am LOCATE THE INDEXED POSITIONS OF FIED RANKS 
LS: ee ee Apri Cir CisZi47l2s2 )- 
FIND THE MIDR K VALUE OF THESE RANKS 
Rs (4 /TIUBE 1534 12 Z) 
A. SE UP VeCcTOR WITH ZEROS AND ONES; ONES WHERE TIE RANKS INVOLVED 
= 


ieee 0 
M<(TI= *PICPPCKI1) 
MM<MM+M 


eee C2 hel ee 
RQ SET OP A pheror WITH ZEROS WHERE TIED RANKS OCCUR 


a TRANSFORM ONES OF MM VECTOR TO MIDRANK VALUE 
+ 
oe ONES OF L VECTOR TO REMAINING UNCHANGED RANKING VECTOR 


RA FILL ron MIDRANK VALUES 
pee 
— 
A 723 772 SAME FOR ANY OTHER TIES INVOLVED BOT WITH NEWLY COMPUTED TI 
> 
V 


Le 


ASSIGN THE RANKS OF THE LEFT ARG. TO TI 
IF NO TIES FOUND QUIT 


FF EFF FWWWWWWWWWWHONYINNYNPNDNHHHE HHH HRBRHBWO QUIN EWP 
£ 


FEE 
POWWOYUNMNMNEANIPrMOWWOIVMMNE ONHOPROWWD YOU FWHM} OWWO~10)tnN FE WN BROU AL. tan 


Crewe wwe we Ww wh ae Oh he wo a a ho ha a a a a os ree el 
Ann lk 


One & 


V ESS (ee aoe ere iC3D3L 


fi] a 2 sie bPocRaM QF XENDALL'S 2 (XEN), SPEARMAN'S 
[et a QZ (Suan) IND (SEMAN). AD SSUSRAL-WAL COS CREWE) om CRECES CUE 
"32004 2IGa? 22GUMENT FOR 7725 GND RECORDS TEE NUMBER OF  OCCUREN 3° oF Fac# 
2] 08 Tie awd ree TOTAL NUMBER OF TIES IN THZ VECTOR. 
6 venga BSSTGN ORDERED VECTOR TO B AND INITIALIZE VALUES 

aad °) 
8 B<AACAAAJ 
9 TIE+p0 
10) +1 


1@5 


Qa CHECKING FOR TIES BY INCREMENTING THRU THE VECTOR 
L3:C«(7T+B )=B(7) 
a COUNT NUMBER OF TIES; IF NO TIES GO TO L2 


+/C 
Sane 0)/Z2 
TrEeTle., (D+1) 
q NCREMENT NEXT T BY THE NUMBER OF TIES ENCOUNTERED PLUS 1 
L2:7+«T+(D+1 


ects Less SHAN SIZE OF ORIGINAL VECTOR GO TO L3 AND START AGAIN AT NEW T 


a eS ee ee ee ee ae es 


NRRPPPHPPpPyY 
OWONMM EWP 


V VAR<VARMW _8;3C;3D;E;D1;E1 
Ds ae ruhebton® IS_ 4 SUBPROGRAM OF THE MANN-WHITNEY TEST (MANW). 


al A 
2} 9a If GENERATES THE RANKING SCHEME USED IN CALCULATING THE 
3] 9 DIFFERENCES IN SCALE (1 ASSIGNED SMALLEST, 2 ASSIGNED LARGEST. 3 
4] A ASSLGNED WEXT GaRer SECOND SMALLEST ETC. 3Y THOS TILL PROPER 
si a SanbEe SIZES TS SESC72D). SAMPLE BT2e Te DASSip iN THe RIGHT Ac. 
ts] Eris00 
— 
3 mel (B22) PLN? 7E00R OF HIDPOINT OF VECTOR AWD ASSIGN TO C 
“= = 
10 batbp: GENERATE RANKING VALUES LEFT HALF FIRST 
ey 
lee EE, 
13 > (CBE) = C)/L3 
{U3 DeD+3 
tS} F<2,D 
16] =(Cp2)<C )/ 22 
173 NOW GENERATE RIGHT ZALF 
15) L3-9ie. 
as Pee ees 
20) £1<F1,D1 
21 *( (923 )=C)/Z5 
22) L620 Tep 1% 
23) Zi~71,D1 
au} > ((aEL)EC)/E5 
255° Uleoies 
26 Bee aes 
27] +((pz1)<C)/L6 
28 isis 2B 20 0/2 VECTOR IS-ODD VALUE CONCATENATE MIDDLE RANK IN BETWEEN HALPS 
30) A EP STZE 78 EVEN CONCATENATE LEFT HALF WITH THE REVERSE OF THE RIGHT 
31) VAR+E, (eE1 
32] +0 
33] L7:VARCE,B,(eE1) 


VY WIL+WILP NN:N3A:P:T:NN:W:PP:NM 
A THIS PUNCELON roca Ciepe eR oM OF THE WILCOXON SIGNED RANK TEST 
a (WISIC). IT GENE RATES THE CUMULATIVE DIST. FOR THE TEST STATISTIC 
a THE GENERATOR USES A RECURSIVE FORMULA. NN=SAMPLE SIZE. 

ee 


peeea SET P EQUAL TO PROB. DIST. WHEN N EQUALS 2. 
cd 

L3:NeN+1 
A<i 


A SET T VECTOR TO PROPER SIZE OF ZEROS. 
T+(+/(1N))p0 


aa (ASN)/L1 IF ASN USE TRUNCATED RELATION TO COMPUTE OCCURRENCES. 
‘> 
"+ (a> >N)/L2 IF A>N USE FULL FORMULA TO COMPUTE OCCURRENCES. 


Pa t A Sea y) IS NEGATIVE TRUNCATE FORMULA TO AVOID A NEGATIVE INDEX. 


A+cA+ 


>I 

R ie Gee ae LARGER THAN THE LENGTH OF P GO TO L6&. 

w2t ela eS j 

ie ee poe POSITIVE: THE RECURSIVE FORMULA CAN S3E USED. 
PCs \ LY epee a 


>L7 

A a5) LARGER THAN THE LENGTH OF P TRUNCATE FUMSTION AGAIN. 
bg: PLAN PU (A-W)) 

fe) a; <(+/¢ Seek INDEX HAS NOT EXCEEDED N(N+1)/2 CO AGAIN. 

> y 

A convern f ANG +p qe CONCATENATE 1 FOR USE IN NEXT ITERATION OR OUTPUT. 
ooo +((+/(1N))o7),1 


P<Ni4+ 
kia Xp ys coanw) 
CHECK IF LENGTH OF INPUT VECTOR EXCEEDS NUMBER OF N'S GENERATED. 
"2 (NNN /L3 


SET IE TET COATT og oe 
WWWWWWMMNDDN ON MOPHMHMPRPPHHMHDOOYNOMiNFWhHe 


MEWANHrPOWOONMMWE WHOM ONMAM FONH OLAass2 at ay 


= ae ee oe ee ee ee ee ee ee ee ee ae eae ee ee ae ae a | 
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APPENDIX G 


LISTINGS OF PROGRAMS USED TO GENERATE C.D.F. COMPARISON 
fea 


v RENTEST sW *ALPHA:B:A:C;TAU:P3NUM:DEN:2;2ZZ3ERRZZ3:M:Z2C:D3AA;PP;NUMC;:I:3H3P 
sFS:S3J:ERRZZC: bei ki 


ATEIS PROGRAM GENERATES TABLES OF C.D.F. COMPARISONS FOR KENDALL'S 3 
SEZ SAMPLE SIZE AND ALPHA VALUES. 


oN=1 a 

nESHa< Seems 2 CO mCo ome moeos O..5 3.28 

2 THIS LOOP INCHEMENTS SAMPLE SIZE. 
Pier P<p0 


M<o0 
KK<p0 


ae ee Bete COMULATIVE DIST AND ASSOCIATED STATS. 
Lek 
Q THIS LOOP CALCULATES ALPHA VALUES AND APPORXIMATIONS. 


e0 TE NORMAL APPROXIMATION. 
ch ee U)x( (Nx (N-1))*0.5) 
DEN«(2x((2xN )+5))*0.5 
Z<+NUM+DEN 
cored » (1 -WORMEDE @) 


COMPUTE NORMAL ne eo WITH CONTINUITY CORRECTION. 
"NOMC«3 x (PAU AA )x ((Nx (N-1))*0.5) 
Z2C«NUMC N 
a a ZC ) 


>(B< 
COMPUTE ERROR DIFFERENCES. 
" BRRZZ+PP- ZZ 
BREZZC*PP-Z2Z2C 
ee PRINT OUT TABLE OF VALUES. 
ata es 


J+195 
Se'TEs? STAT. VALUE ' 
—, ,PROBCB 2 BJ); FOR SAMPLE SIZE EQUAL TO ',(2 0 BN),! 


SESE 0O4V0(18p197),7,619(8p197),/7] 
8 


1n8. 
< 3.9999 el) OFMT(1 7 pKX) 
Q 


Ci 
Fe! ' 267 5 ePP{C]) 


»(C<7)/L9 

L7: hei, rt§ PS 
D+D+ 
ay N2,N3 25D] 
> 

N1iss :S<! yee xAcT PF. 


$53 
N2:S<'ERROR; NORMAL : 
oo 

> 


WNHPOWDNODUWNEWNHMOWONHDMWNEWNPOWOOND NEWNPOWDONMUMNEWHPOW ONDINE WHhPHMOWO~1NUM EWN Ot-aa waa ans 


8A en A nk Oe eae aaa eS aS Se ae SSeS aS See ee ee Se ee ee ee ee ee ee ee eee nA Ad 
Pret tlh 
ee 
tyYD >; 
he 
ty xt 
C7 Maere 
ote ve 
ae 
=e 
uN 
CC 
x ty 


TTT TNT TE NT OI AIA Loa 
SNNNDODOOONAOHNOOHUIINIMNIMWANMNNMNNE EEL FFE FFE EWWWWWWWWWWYONYNNYNMWNNDYNE HbA ERP 01MM FRR 


VOn 


Diez 
iis 


NNN 
ONMWE 


] N2:S*'ERROR: NORM. W/CC } 
2 RELIC 


Lu:M 
+ (N<14)/E1 


V RWTEST:A >PP:PC;NN: ALPHA ; N:K:P:P:H;PVALUE:;PVAL;PP1;PF;PVF:P2VF1;D;=28R 
HB; Tbe? EERE: SCDP:CC:KK3Llid:FF3S:M 


A TEIS PROGRAM CENERATES TABLES OF C.D.F. COMPARISONS FOR THE KRUSKAL- 
a WALLIS TEST. 


A« YuUY 
CC+2 
Bea 0.07 0.02750 .203 50.205 (0208 Osis sO re 


b1:4«<NN 

8 eee eee TO GENERATE EXACT DISTRIBUTION FOR SELECTED SAMPLE SIZE. 
ood 

N<+/A 

CC<CCt+1 


M+<+p0 

PP<o0 
KK<«o0 
PC+90 


FiunpP 


OWO~ JTOONF Whoops OuUn1 A_p_A_aA aA 


9 
ChFeet2:] 
D4: B<«B+1 
a LOCATE POSITION OF EXACT CDF VALUE < ALPHA. 
ea eee 
ae )/L3 


geeePUai ol CORRESPONDING TEST STATISTIC VALUE. 


Q 

re ° 

L3 e aoe 
a 


RR+ 
BECORD EXACT VALUE OF CDF. 
BVALUE<P 2 j Di 

PP+PpP Pp E 

A COMPUTE SonRESrOND IHC P-VALUE USING CHISQ APPROX. 


TT a 
Lok A Mk A 


AL 
Chet CORRESPONDING Pe VALUE USING F APPROXIMATION 
Fe ((N-K)xH)+ ((K-1)*%((NW-1)-#)) 
PyPet (ke 1) (WSR) PDISTN PF 
PF+PF,(1-PUF 
A conph re C. REESPONDIN P-VALUE USING F APPROX W/ 1 LESS D.F. IN DENOM 
PVF1+((K-1 N-K)-1) EDIT een P 
PFi+PF1i,(1- bey5 
>(B<7)/Lu 


< 
ERRH<PP-PC 

ERRF«PP-PF 

ERRFi«PP-PFi 

Ome OUT TABLE OF VALUES. 


FEE FWWWWWWWWWWhNDNNYNYNNDNHPHRPHPBRBRHBH BPH DOOM 


WHPODDNOONE WOH OW ONO F wh 


ov 
LAABA AA 
D 
eas 
sare 
acy + 
“~~ 
ae 
le 


«tao! 


S+ 1 TESL Sa ede : 
M<éM, BCH 2 : DAVC46I, 14; FOR A GROUP OF 3 SAMPLES CONSISTING. OF 4¥& 


, 4. AND ! eces OBS. 
ee 


N1:S*«!ALbHA’ 
KR+PP 
Je'+!? 


A IAT IS rire ieee 
FEEL 
NOUN 


“Ug 

ug 

50 

51 

52 

53 

5u 

oe) £2: few 03 (a7 11) 7), cup Conan me 
<= p Oop 

095 US ee za 

58 Fepo” 

59 H+160S 

60 F+«F, 

61 C+0 

62! +(Dz0)/ZL10 

Odea 1 owed 

54) FFs! Lon 7 u oKRLC]) 

orm Guar ad tated ak 

66 +(C<7)/L11 

67 +512 

68] £210:C+C+1 

69 Fe! |!,(7 5 &KK(C]) 

76 F+F ,F 

7A >+(C<7)/L10 

72} 212:M+«M,(1) F 

73 ee 

74 I< 

75 SCN ,W2 N3,N8,L5)(€D] 

76) VALUE 

174 

78) 
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+L8 
N2:S*!EZRROR; CHISQUARE' 
KX+ERRH 
+18 

N2:S*'2RROR; F DIST '! 
KESERRF 

Nu:S«!ERROR; F W/71 DF! 
KRSERRF 1 

= 


L3:% 
72s 


CO ee 
£0 00 00 00 00 00 00 00 00 00 CON) 


OwWoo~ IM FWHM Ow 


(A, A_A_A_A_ A A 


v EWTESTSM A:R2;PC3;NN:ALPHA:N:K:Y:W:H:sYY:PVAL;PF1:8H:MM:M:X:ERRC:ERRF1;F; 
PVsBsKK:DsPVFisQ3FF:U;sC:Iid:sS 


Bega snOChAN GENZTAATES C.D.%. COMPARISONS FOR #RE XROUSKAL-WALLIS TEST 
a BASED _ON 30000 HANDOM VECTORS ne UO RANKINGS WHICH SIMULATE COMPARING 
a 39 SAMPLES OF 8 OBSERVATIONS EACH 


a GENERATZ A SANDOM VECTOR OF RANKINGS. 


td EACH SAMPLE. 


M<o00 

Bor-ae~ 0.02 0.02 0.03 0.05 0.08 0.13 0.18 
K+p0 

PF1<90 

PC<+p0 

L3 BB+ di 

A DETERMINE DESIRED H VALUE CORRESPONDING TO ALPHA VALUE. 
ee ete) 


KK«K 
Q COMPUTE ae el P-VALUE USING CHISQ APPROX. 


NPOWDONANHE WOROWONMNMWEWHMOWONODUFE WN OLA-aAAr tn 


FEEFE 
SOOne 


Lit 6 oe A tut A _ tii ft nd A On A Ue en A tk Ae hk 
Dp 
™ my 
AANA 
ants $ 2) 
mtr +l 
WO SE 
On + 
FAW 0 
be toh, 
baht 
ty * PO 
nN 


PC*PC (1-PVAL 
u8} an COMPUTE HORBESPONDING P-VALUE USING F APPROX W/ 1 LESS D.F. IN DENOM 
ug Be Gy Cats SR =4 
50 Ee q ws K)-1 FDISIN F 
51 Pie -PYF 
52 47 Se 
53 ERRCCALPHA =PC 
5u) ERRF1«ALPHA-PFi 
55 PRINT OUT TABLE OF VALUES. 
36, I+0.1 
S81) S<'T°UST STA”. VALUE! ; 
S9%° Mem, BeODeta> ' C27iusi,.'1: BASED ON ',(5 9 oW),' GENERATED a@‘''S FOR $ 
SAMPLES OF 8 OBS. EACH. I 


Je+ie-! 


Da 
Le: ae ie!) ae omic" ) yu 


Fe 90 
H<160S 
F<F ,# 

C+0 
=(0z0)/440 
£11:C< 

FF<! a7 3 seKK{(c])) 


TTT TTT I ET IIS ume 
FEE EWWWWWWWWWWHNONYNDDPNNHPH HEHE bt WwW00NMG EWP 
GQ) 


SNHNAMNMOMHMMH HOM 
POW ONIHAMWEWNHO 


F<«<F 
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COGG7 7/011 
Li 


ee ee 
9 WOON SY SYQy 


NPOWDYAMNM FWwp) 
+ 
o~ 
Cc 
A 
x) 
w 
™~ 
ts 
pa 
OQ 


N2:S«'ERROR; CHISQUARE'! 
KK<«ERRC 


>Z8 
N3:S<'ERROR; F W/ 1 DF! 
KK+ERRF1 


Be 


de 990° 30000 1 


atvtAcrl 


+ (Asi )/L1 


L3 


ET AEA aioe ooo 


DN FWHFOWMONMNMN Fw 


CHAWLA 


LO 10:10: 10 t0 LOW © 00 00 0000 00 00 


Z3;2RR2I2 


;M:220;D: 23 NUM 
ee tC ae ~ benGaE Behone 2 ERR ee R°T:Z 


m Gus 
V MANT eee eee 1 NN 3 MMs ene ‘SiR BRR eae 


= mh. ine 


 EBArCZos Kk 


> 

eo Ee 
m 

—_ ; 


: 
5 
a! 


P. 


a TEIS PROGRAM GENERATES TABLES OF C.D.%. COMPARISONS FOR THE MANN- 
A WHITNEY TEST. 


oe ; SET SAMPLE SIZE AND ALPHA VALUES. 


ADDER 0.01 0.02 0.03 0.05 9.08 0.12 0.7 
— s a oake 
7 PETS LOOP CNCRSMENTS SAMPLE ST 
: La: eee Pe 


P<NN MANWP MM 
Po pepe PHIS LOOP CALCULATES ALPHA VALUES AND APPORXIMATIONS. 
7Be 


4047 (PSC) 
A+«(+/(PSC)) 


Ke RK.U 
eae Pray 

COMPUTE NORMAL APPROXIMATION. 
" NUM«U-~ CON XMM 
DEW = ( (NWxMiIx (NN SHM+4))#12)%0 5 
Z+NUM+DEN 
eee NORMCDF Z 

COMPUTE NORMAL APPROXIMATION WITH CONTINUITY CORRECTION. 
"NUMC*(U+0- 5 )-(€ (NNxXMM )42) 
ZC+NUMC+DEN 
eco RaRMCDP ZC 

COMPUTE STUDENT T Ae Dee 
"DEWOMe ( (WN Hit ie (Z*2))+ (NN+MM~2))x0 
1 2x2: DENO 
ee aeane N+MM~2) OT Sele ee 
’ a ccneHtE SLUDENT « APPROXIMATION WITH CONTINUE 22 CCR 2 Clo. 
J d <= ( +0 
; DEWNOMC+( (( (NN+MM=4)x(DENK2))¢ (NN+MM-~2))~-( ((NUM+0.5)*2)%(NN+MM-2)))*0.5 
TC+<NUMT+DENOMC 
rrc<2tC ((WN+MM-2) TDISTN Te) 


>(B<7)/L2 

R COMPUTE ERROR DIFFERENCES. 

ERRZZ+PP-22 

ERRZ2ZC+PP-2Z2C 
ERRTT+PP-TT 
ERRTC+PP-TTC 

ERRAVE+PP- (622420 42) 
ERRTCZC+PP-((Z2Z2C+7TC )#2) 

A ok PRINT OUT TABLE OF VALUES. 
€ 2 


9 J*195 


— 


2 
3 

yu 

5 

6 

7 

8 

9 

: 

: COMPUTE CUMULATIVE DIST AND ASSOCIATED STATS. 
u 

5 

6 

Z 

8 

S 

0 

2 

2 


ft 


‘a 
‘ 
a an 


_S~ — STAT. VALUE 
PROBCO. = Gl; FOR SAMPLE SIZES N EQUAL TO ',(2 0 SNN),' AND 


M<éM 
ijt wn (2 0 8MM),7 
£3; 
i 


67 5 sPP(C)) 
3 FS 


oe 
>(N1,N2,N3,N4 NS ,N6,N7,L5) (DI 
Niise! 'pxAcr’c.b.F. 
» 


oo 
a: aero NORMAL , 


oa: SAG oe 


N32: :5=!22R0R: NORM. wW/CC' 


PER+EZ222 


4226 of oe c 


Nu: 1 S«'ERROP; LeDisT : 
PP+ERRTT 


aL3 
Ns: Se LERROR; © W/CC : 
oe PPSERETC te 


ye 1 S41 ERROR: ave T/T ' 
PPCEREAVGE 


> 
NS SOUERAOR: @AVE C7 2G" 
eee 

> 


On F WH Otlvta_A_a_n An on nL An nN AR_ALA R) ol 

t 
“J 
= 
4 
= 
~~» 

es 
f}- 
to 


ass. AAA 


TT IE IAI ITI oe 
MPRRPRH HF 1000101010 10.00: 10 10 0 WD 0 000000 000 OOO NNN NNN OMMOMOMa od 
OOVKO OOO MO INMNEWNPOWDNIOWN EWM HONMDIMWELWNHOMDONIHMNEWN RO 


Lucu 
gee 9/1 


V ee gs 26 OD EERD DOM 220 DEEL HSPs PS 3S 30 SBRRZZCA2C 


1 
: A THIS PROGRAM CENERATES TABLES OF C.D.F. COMPARISONS FOR THE SIGN TEST. 
2 re SET SAMPLE SIZE AND ALPHA VALUES. 
5 ALPHAS O% oe 0. 0.06 ae) 2 eee 2 
7 Pers. Loop INCREMENTS SaMeEE’ Size. 
8 Li:PP+p0 
9 M+o00 
-10 KK+p0 
alge ZZ2+p0 
12 2oC+p0 
“ti 3 N+N+1 
14 D+0 
1S B<+0 
ae ae pono pose CUMULATIVE DIST AND ASSOCIATED STATS. 
€ 
78) a THIS LOOP CALCULATES ALPHA VALUES AND APPORXIMATIONS. 
19} £22:38+«53+1 
f20) S#@LPFAC8) 
b218 ete / (Psc, ) 
Ze imam) 
"235 ZF<+22,2 
24 PP+PP, ,PCA] 
250) oA COMPUTE NORMAL APPROXIMATION. 
26 2+ (K-(0.5xN) )#(0.5x(Nx0.5)) 
27 35052 CNORMC Grr 2) 
28:\0A ae NORMAL APPROXIMATION WITH CONTINUITY CORRECTION. 
293 Z2C+((K+0.5)-(0 a *0.5)) 
30 22C*22C CRORMCDF DC 
Ba). 2(B<7)7L2 
32 COMPUTE ERROR DIFFERENCES. 
33 " ERRZ7«PP- 22 
3u pee cere ZLC 
35 PRINT our TABLE OF VALUES. 
36 pe 


111 


Eee oes 


INMNMNONONNEEF EF EEE FEE Www 
INEWNRPOWDO-"INMVEWNHHO wWwor~) 


—MowL dh 


mw) ow . 2 ee 


ASNNMHAIMNHONMNMNIMMHA nnn 
NOPOWO OO {AN FWP OWoO~Id 


CALA A AN A Ne 


MFFE FEF EEE FEWWWWWIWWIWWWNHNYNNNYNNANANPRPPPHPPHRHEHHOOUMUOFWNDE 
OOONAMN FEWNHPOWDNIHAAEWONPROWDONMDMNFEWNHPOWOMNMDMNFWNP OUA «a2 


LLL TNE oor on 
Rill ll lll nl 


J*#195 
Se!'!TEST STAT. VALUE ! 
M<M,'! ,PROBCK S K£]; FOR SAMPLE SIZE EQUAL TO ',(2 0 8N),! 


OAVC(¢189197),7,61p(89197),/7] 


29 >' OFMT(1 7 pkKK) 


+1 
(7 5 sPPEC] ) 


>(C ae 
Li: Det (eG FS 


oa 


ath N2,N3,L5)CD 
n1:S<!ExXAc?T Bo pte 
J+198 


+L3 
N2:Se'ERROR; NORMAL us 
PP<ERRZZ 


N3:S¢'ZRROR; NORM. W/Co 
P<ERRZZC 


Lu3sM 
= (N<26)/b4 
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